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DEVELOPMENT  OF  COMPUTER  SOFTWARE 
FOR  THE  ANALYSIS  AND  DESIGN  OF 
MODERN  CONTROL  SYSTEMS 

CHAPTER  1 

PRELIMINARIES 

1.1  Introduction  to  the  Control  Systems  Software  Package 
The  classical  techniques  of  graphical  analysis  and  design  of  modern 
control  systems  involve  the  generation  of  sketches  by  hand.  In  all  but 
the  most  simple  systems  this  can  be  a  tedious  and  difficult  process.  The 
purpose  of  this  research  was  to  develop  computer  software  that  automates 
this  task.  The  software  programs  developed  include  the  Root  Locus,  Bode 
Plot  and  all  necessary  support  routines  needed  to  fully  automate  the 
development  of  these  graphs.  This  software  is  called  collectively  the 
Control  Systems  Software  Package  (CSSP). 

1.2  Computer  System  Requirements  for  CS.SE 
The  CSSP  programs  are  written  for  use  on  IBM  and  compatible  personal 
computers  with  a  graphics  card  and  monitor  installed.  On  machines  with 
an  Enhanced  Graphics  Adaptor  (EGA),  or  higher  resolution  monitors,  the 
graphs  are  produced  in  color.  Computers  with  a  Color  Graphics  Adaptor 
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(CGA)  are  limited  to  black  and  white  plots  due  to  the  low  resolution  of 
the  monitor  when  plotting  in  color. 

A  printer  capable  of  supporting  a  graphics  screen  dump  is  required 
if  a  printout  is  necessary,  but  the  printer  isn't  specifically  needed  to 
operate  the  program.  A  numeric  coprocessor  (8087,  80287  or  80387)  is  not 
required,  but  is  highly  recommended.  This  coprocessor  increases 
computational  speed  considerably.  All  hardware  specific  settings  are 
totally  automatic  and  are  invisible  to  the  user. 

1.3  The  C  Computer  Language 

All  of  CSSP’s  source  code  is  written  in  the  C  language  and  is 
included  for  future  study  and  possible  refinement.  Since  the  C  language 
is  very  potable  and  is  not  tied  to  any  particular  hardware  or  operating 
system,  very  little  effort  would  be  required  to  rewrite  this  software  to 
operate  on  other  types  computer  systems .  This  is  true  for  all  but  the 
video  specific  functions,  such  as  the  windowing  routines. 

The  CSSP  source  code,  provided  in  the  Appendix,  was  written  for  and 
compiled  with  the  Micorsoft  C  Compiler  version  5.0.  The  windowing 
software  used  to  develop  the  screens  and  input  windows  was  "The  Window 
BOSS"  Revision  :  08.15.88  written  by  Phil  Mongelluzzo  of  Star  Guidance 
Consulting,  Inc..  All  the  necessary  source  code  and  libraries  needed  to 
reproduce  CSSP  are  provided  on  a  5  1/4  inch  diskette,  which  is  available 
at  the  Aerospace  and  Mechanical  Engineering  office,  University  of 
Oklahoma . 

This  source  code  is  provided  for  educational  purposes  only,  and  is 
not  to  be  used  in  any  type  of  commercial  application  without  the  prior 
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written  permission  of  the  author.  Copyright  (c)  1988-1989  by  Carl  F. 
Adams,  all  rights  reserved. 

1.4  Thesis  Overview 

The  remainder  of  the  thesis  is  devoted  to  the  discussion  of  the 
classical  techniques  of  analysis  and  design  using  the  Root  Locus  and  Bode 
plots.  Along  with  these  techniques  will  be  a  discussion  of  the  equivalent 
computer-aided  techniques  used  in  the  CSSP  program.  The  following  is  a 
short  description  of  the  remaining  chapters. 

Chapter  2:  The  Root  Locus  Method.  This  chapter  discusses  the  Root  Locus 
method  of  describing  the  transient  response  of  a  control  system  by 
determining  the  location  of  the  roots  of  the  characteristic  equation  as 
the  open-loop  gain  is  varied  from  zero  to  infinity.  This  chapter  also 
introduces  a  set  of  ten  rules  used  to  help  in  the  plotting  of  root  ? ocus 
by  hand. 

ChapL.gr  3.; _ Computer-Aided  Root  Locus  Plotting.  This  chapter  introduces 

the  two  most  common  methods  of  computer-aided  root  locus  plotting.  They 
are  the  numerical  polynomial -factoring  and  angle  testing  methods.  It 
talks  about  why  the  polynomial -factoring  method  is  the  simplest  algorithm 
to  program,  but  the  price  of  simplicity  is  a  longer  computer  run  time. 
It  also  explains  why  the  angle  testing  method  is  the  harder  of  the  two 
methods  to  program,  but  is  usually  the  more  efficient  method  for  computer 
run  time . 

Chapter  4:  Root  Locus  Program.  This  chapter  discusses  the  root  locus 
program  used  in  the  CSSP.  It  uses  the  angle  testing  method  for  plotting. 
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Chapter  5:  Bode  Plots.  This  chapter  serves  as  an  introduction  to  Bode 
Plots.  This  is  a  graphical  technique  that  plots  a  control  system's  gain 
amplitude  and  phase  angle  response  curves  against  the  input  frequency. 

Chapter  6:  Computer-Aided  Bode  PLots .  This  chapter  describes  the  Bode 
plot  algorithm  used  in  the  CSSP. 

Chapter  7:  Summary.  This  is  the  thesis  summary. 


CHAPTER  2 


THE  ROOT  LOCUS  METHOD 

2.1  Introduction  to  the  Root  Locus  Method 
The  root  locus  method  is  a  powerful  graphical  tool  used  in  the 
analysis  and  design  of  linear  control  systems.  This  method  determines  the 
location  of  the  roots  of  the  characteristic  equation  as  the  open-loop  gain 
is  varied  from  zero  to  infinity.  By  knowing  the  locations  of  these  roots 
it  is  possible  to  not  only  know  if  a  system  is  stable  but  also  the  degree 
of  its  stability.  This  method  is  another  way  of  describing  the  transient 
response  of  the  system.  If  a  system  is  unstable  or  has  undesirable 

response  characteristics  the  root  locus  method  provides  possible 
qualitative  ways  to  improve  it. 

2.2  The  Rool  Locus  Method 

Consider  the  simple  single-input  s ingle- output ,  closed-loop  system 
shown  in  Figure  2.1. 
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The  transfer  function  of  this  system  is 

C(s)  KG(s) 


T(s) 


R(s 


1  -1-  KG(s)H(s) 


(2.1) 


and  the  characteristic  equation  of  this  system  is  obtained  by  setting  the 
denominator  equal  to  zero.  Thus  the  general  form  of  the  characteristic 
equation  is  given  by 

1  +  KG(s)H(s)  -  0  or  KG(s)H(s)  =  -1  (2.2) 

or,  in  polar  form, 

KG(s)H(s)  -  1/180°  (2.3) 
where  K  is  the  gain  of  the  system.  The  characteristic  equation  contains 
all  the  information  about  a  system's  frequency  response,  time  response, 
and  stability.  In  other  words,  any  value  of  s  that  causes  the  open-loop 
transfer  function  G(s)H(s)  to  have  a  gain  of  unity  and  a  phase  shift  of 
180°,  and  by  definition  a  point  on  the  root  locus,  will  be  a  pole  of  the 
closed-loop  system.  This  is  a  very  useful  fact  in  that  evaluation  of  the 
usually  simpler  open-loop  transfer  function  tells  how  the  poles  of  the 
closed-loop  system  behave.  Therefore,  by  evaluating  G(s)H(s)  in  the 
complex  s  domain  and  finding  the  values  of  s  that  make  equation  (2.3)  true 
will  determine  how  the  closed-loop  poles  are  influenced  by  changes  in  any 
of  the  system  parameters. 


To  determine  if  a  given  value  of  s  is  a  point  on  the  root  locus  for 
some  value  of  K  between  zero  and  +*>,  it  is  only  necessary  to  determine 
whether  or  not  the  angle  of  KZ(s)/P(s)-KG(s)H(s)  is  an  odd  multiple  of 
180°  and  the  magnitude  must  be  equal  to  unity.  This  is  known  as  the  angle 
criterion  and  the  magnitude  criterion  respectively  and  are  shown  in 


equations  (2.4)  and  (2.5a). 


=  y  (angles  of  the  zeros)  -^(angles  of  the  poles) 
=  (2k -  1 )  1  80  .  k=l  2.... 


(2.4) 


K 


=  1 


(2.5a) 


which  gives 


product  of  the  pole  distances  _  )P(s) 
product  of  the  zero  distances  iZ(s) 


(2.5b) 


2.3  Rules  for  Plotting,  a  Rool  Locus 
To  help  in  the  plotting  of  root  loci  a  set  of  conventional  rules 
have  been  developed.  These  rules  have  been  thoroughly  discussed  in  [1]. 
Here,  a  brief  description  of  these  rules  is  presented  and  applied  to  two 
systems  with  the  following  open-loop  functions, 


example  1: 


_ K _ 

s(s  -r  5)(s  -r  7) 


example  2: 


_ K(s  6) _ 

s(s  +  4)(s2  t  4s  -  8) 


Root  locus  plots  of  the  above  open-loop  functions  are  shown  in  Figures 


Stable  Gain  Range:  8<X<428 

K 

Unstable  Gain  Range:  420<X<« 

(sHs*5)(s+7) 

Figure  2.2  Root  locus  plot,  exaaple  1 


Figure  2.3  Root  locus  plot,  exaaple  2 


RULE  1:  There  is  a  single-valued  branch  of  the  root  locus  for  each  of  the 
characteristic  equation  roots  and  the  total  number  of  branches  is  equal 
to  the  number  of  open-loop  poles. 

Example  1:  There  are  no  zeros  and  three  poles,  therefore;  there  is  a 
total  of  three  locus  branches . 

Example  2:  There  is  one  zero  and  four  poles,  therefore;  there  is  a  total 
of  four  locus  branches. 

RULE  2:  Each  branch  of  the  root  locus  starts  at  a  pole,  where  K=0 ,  and 
ends  at  a  zero,  where  K  -  +a>.  If  the  number  of  poles  exceeds  the  number 
of  zeros,  there  will  be  zeros  at  infinity  equal  in  number  to  the  excess. 
Excess  zeros  similarly  mean  poles  at  infinity. 

Example  1:  There  are  three  more  poles  than  zeros,  which  means  all  three 
of  the  branches  go  to  infinity. 

Example  2:  There  are  three  more  poles  than  zeros,  which  means  three  of 
the  four  branches  go  to  infinity. 

RULE  3:  Along  the  real  axis  the  locus  includes  all  points  to  the  left  of 

an  odd  number  of  real  poles  and  zeros:  no  distinction  is  made  between 

poles  and  zeros,  and  complex  poles  and  zeros  are  neglected . 

Example  1:  There  will  be  points  on  the  real  axis  from  0>o>-5  and  from 
- 7>o>-® . 

Example  2:  There  will  be  points  on  the  real  axis  from  0>a>-4  and  from 
-6>o>-« . 

RULE  4:  If  the  number  of  poles  np  exceeds  the  number  of  zeros  n,,  then  as 
K  approaches  infinity ,  (np  -  n,)  branches  will  become  asymptotic  to 

straight  lines  intersecting  the  real  axis  at  angles  given  by 
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If  n,  exceeds  np, 
above . 

Example  1: 

Example  2: 


_  (2^  - 

1 )  1  BO1 

(np  - 

n2) 

as  K  approaches 

(2k  - 

1  )  1 80 

(3  - 

0) 

(2k  - 

1)180 

(2.6) 


k  =  0,  1,  2... 

zero,  (nt  -  np) 


(4  -  1) 


(2k  *  1)60.  k  -  0.  1.  2 


I2k  -  1)60‘.  k  -  0 .  1 . 


In  both  examples  there  will  be  three  asymptotes,  one  for  each  infinite 
zero,  intersecting  the  real  axis  at  angles  of  60°,  180°  and  300°. 


RULE  5:  The  asymptoter  of  Rule  4  will  intersect  the  real  axis  at  a  point, 
called  the  center  of  gravity,  given  by 


eg 


V 


poles 


-  V 


zeros 


(np-  nz) 


(2.8) 


Example  1: 


eg  = 


0-5-7 
(3  -  0) 


-4 


Example  2: 


eg  = 


+  2j  -  2 


2J 


(4 


1) 


RULE  6a:  A  breakaway  or  arrival  point  st  on  a  real  axis  can  be  found  from 
the  equality 
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(2.9) 

_ 1 _ 1 _ ^  2  (sb  -  n) 

s'b  -  P,  ~  st,  4  z!  ^  <sb  ~  a)2  +  w2 

-  !  _  1  _  2  (sb  -  a) 

;Sb  *  Pr;  "  Sb  ~  Z  r '  T  (-%  -  «)2  -  «2 

where  the  terms  on  the  left  side  of  the  equation  come  from  the  poles  -p, 
and  zeros  -zx  to  the  left  of  the  breakaway  point;  the  terms  on  the  right- 
hand  side  come  from  the  poles  -p,  and  zeros  -z,  to  the  right.  Positive 
signs  are  associated  with  the  zeros  and  negative  signs  with  the  poles. 

RULE  6b:  A  breakaway  point  occurs  when  K  is  at  a  relative  maximum;  an 
arrival  point  occurs  when  K  is  at  a  relative  minimum. 

The  breakaway  point  can  be  found  by  solving  the  characteristic  equation 
(2.2)  for  K  and  then  taking  the  derivative  with  respect  to  s  and  setting 
it  equal  to  zero.  Then  solve  for  the  resulting  roots. 

Example  1:  The  roots  of  (2.2)  are  S--6.1,  -1.9.  Since  the  point  S--6.1 
is  not  on  the  root  locus,  see  Rule  3,  the  breakaway  point  must 
occur  at  S--1.9. 


Example  2:  The  two  non-complex  roots  of  (2.2)  are  -3.1,  and  -7.3.  The 
breakaway  point  is  approximately  -3.1  and  the  breakin  point 
approximately  -7.3. 


RULE  7:  Branches  of  the  root  locus  are  symmetrical  with  respect  to  the 

real  axis  since  all  complex  roots  appear  in  conjugate  pairs. 

This  implies  all  root  locus  points  that  lie  above  the  real  axis  will  have 
matching  mirrored  points  below  the  axis. 


RULE  8:  Two  branches  of  a  root  locus  breakaway  from  or  arrive  at  the  real 
axis  at  angles  of  ±90°. 

For  three  or  more  branches  use  Rule  9  below  to  calculate  the  departure  or 
arrival  angles. 


Example  1:  The  breakaway  point  at  S--1.9  breaks  from  the  real  axis  at 
±90°. 

Example  2:  The  breakaway  point  at  S--3.1  and  the  breakin  point  at  S--7.3 
break  from  the  real  axis  at  ±90° . 

RULE  9:  The  angle  criterion  can  be  used  to  find  the  angle  of  departure 
from  a  complex  pole  and  the  angle  of  arrival  to  a  complex  zero  by 
selecting  a  trial  point  so  close  to  the  complex  pole  or  zero  that  the 
angles  from  the  other  poles  and  zeros  are  unaffected . 

Example  1:  There  are  no  complex  poles  or  zeros  for  this  example; 
therefore,  this  rule  is  not  applicable. 

Example  2:  Figure  2-4  and  equation  2-10  illustrate  how  the  departure  angle 
of  -63.4°,  for  the  pole  at  the  point  (-2+2j),  is  derived. 

-aid-a, -aj-aj+Qi.  -  180°  (2.10) 

RULE  10:  If  np-n^2,  then  the  sum  of  the  roots  of  the  characteristic 

equation  is  constant  and  equal  to  the  sum  of  the  poles. 

This  rule  is  limited  in  use,  but  might  prove  to  be  useful  under  certain 
circumstances.  It  isn't  needed  for  examples  1  or  2 .  Refer  to  [1]  for  an 
example  using  this  rule. 


Figure  2.3:  Departure  Angle 


CHAPTER  3 


COMPUTER-AIDED  ROOT  LOCUS  PLOTTING 

U. _ Introduction  to  Computer-Aided  Root  Locus  Plotting 

Using  the  ten  rules  explained  in  Chapter  2  it  is  possible  to 
construct  root  locus  plots  with  a  reasonable  amount  of  accuracy;  however, 
since  this  method  involves  hand  construction  and  the  use  of  a  calculator, 
it  can  become  confusing  and  tedious.  To  improve  accuracy  and  to  make 
better  use  of  a  system  designer's  time,  a  computer-aided  approach  to  root 
locus  plotting  was  developed.  This  chapter  introduces  the  development  of 
the  numerical  polynomial -factoring  and  angle  testing  methods  used  in 
computer-aided  root  locus  plotting.  This  chapter  also  explains  why  angle 
testing,  which  is  used  in  the  root  locus  program  discussed  in  chapter  4, 
is  usually  the  optimum  method  to  use  for  all  but  the  simplest  systems. 

3.2  Numerical  Polynomial  Factoring 
One  method  of  computer-aided  root  locus  plotting  is  to  repeatedly 
factor  the  equation 

1  +  G(s)H(s)  -  0  (3.1) 

as  K  is  varied  in  step  increments  from  zero  to  infinity.  Figure  3.1  shows 
a  simplified  programming  flow  chart  for  a  program  using  the  numerical 
polynomial -factoring  method. 
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(  STOP  ) 


Figure  3-1  Polynomial -Factoring  Flow  Chart 


A  root  of  equation  (3.1)  is  any  value  of  s-x+jy  that  makes 


D(s)  -  l+KG(s)H(s)  -  0 


(3.2) 


true  for  any  fixed  value  of  K.  This  will  only  be  true  if  and  only  if  both 
the  real  and  imaginary  parts  of  D(s)  are  equal  to  zero. 


-16- 


That  is 

Re[D(s-x+jy) ]  -  0 
Im[D( s-x+jy) ]  -  0 


(3.3) 


or  alternatively  if 


| D(x+ jy ) |  -  0  (3.4) 
The  programing  needed  to  solve  (3.3)  or  (3.4)  is  simple,  but  using  this 
method  to  locate  the  polynomial  roots  is  a  time  consuming  approach.  This 
is  a  two-dimensional  problem  in  x  and  y  and  repeatedly  factoring  this  type 
of  polynomial  can  be  a  lengthy  process.  This  method  is  especially 
inefficient  when  the  gain  factor,  K,  becomes  very  large.  It's  for  this 
reason  that  the  angle  testing  method  was  developed. 


_ Method  of  Ap£ie  Testing" 

Points  s-x+ jy  on  a  root  locus  are  those  for  which  the  angle 

/GH(s-x+jy)  -  180°  (3.5) 

or 

Im[ D( s-x+jy) ]  -  0  (3.6) 

The  angle  testing  method  involves  using  (3.5)  to  find  a  locus  departure 
angle  starting  from  one  of  the  open-loop  poles.  Finding  the  open-loop 
poles  and  zeros  requires  only  one  pass  of  the  polynomial  root  finding 
routine,  when  K  -  0,  instead  of  the  multiple  passes  of  the 
polynomial -factoring  method. 

The  departure  angle  is  used  as  the  starting  point  for  a  process  that 
iterates  the  angle  of  G(s)H(s)  until  the  next  point  is  sufficiently  close 
to  satisfy  (3.6)  to  within  some  given  tolerance.  This  means  that  the 
two-dimensional  problem  of  x  and  y  is  reduced  to  a  one -dimensional  angle 
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problem.  This  process  continues  along  one  locus  until  the  plot  limits  are 
exceeded.  Then,  starting  at  a  new  open-loop  pole,  the  whole  process  is 
started  over  until  all  loci  are  complete. 

The  program  development  needed  to  solve  (3.5)  and  (3.6)  is  much  more 
involved  than  that  of  the  polynomial  factoring  method;  however,  the 
savings  in  computation  time  between  the  two  methods  can  more  than  make  up 
for  this  difference.  This  is  true  for  all  but  the  simplest  problems. 
Figure  3.2  shows  a  simplified  programming  flow  chart  for  a  program  using 
the  angle  testing  method. 
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Figure  3-2  Angle  Criterion  Flow  Chart 


CHAPTER  4 


ROOT  LOCUS  PROGRAM 


4.1  Introduction  to  Root  Locus  Program 
This  chapter  discusses  the  root  locus  plotting  program  which  was 
developed  to  aid  the  control  system  designer  in  the  design  and  analysis 
of  control  systems.  The  algorithm  described  in  this  chapter  follows  the 
work  of  J.A.  Borrie  [2],  Borrie’s  work  was  based  on  the  original  work 
done  by  P.  Atkinson,  and  V.S.  Dalvi  [5],  This  program  uses  the  angle 
testing  method  as  explained  in  chapter  3.  This  computer  program  is 
written  in  the  C  language  and  all  source  code  is  provided  in  Appendix  A. 

4.2  Algorithm  Development 

To  start,  it  is  helpful  to  write  the  system  open-loop  transfer 
function  in  the  form 

G(s)  -  G(or-ji;) 

M 

KFIrra-O- J  U'--z.)] 

- - -  (4.2) 

|  |  j  (  (7  £7pk)  *  J  ( *** ~ 
k  1 

where 

K  -  closed-loop  gain 
M  -  number  of  complex  zeros 
N  -  number  of  complex  poles 
(at i+jw„)  -  complex  zeros 
(opk+jwpk)  -  complex  poles 
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Equation  (4.2)  can  be  written  in  the  form 

M  N 

K II  [  ( rr  a.  )  -  j  ( i;  -  i;,, )]  if  F(  n  -  a )  -  i  ( u  - 
i  1  k=  1 

G(o~j  cc )  =  - — — - - — - - - —  - — — 

I  lKcr-^r- 

k-  1 

Simplify  (4.3)  by  writing  it  in  the  form 


(4.3) 


=  K[X-jY]  <4-4> 

Where  X  is  the  real  part  and  Y  is  the  complex  part  of  (4.3).  To  evaluate 

(4.4),  start  with  the  first  pole  (oPi+ jwPi ) , 


Xj-jYj 


(rr-ap:  )-  j(i— wp,) 

2  2 
(  O  Op  i  )  “  ( -*-  p  ]  ) 


(4.5) 


Finish  evaluating  the  remaining  poles  using  the  recursive  formulas 


(o  OpGY^ _■  pk) ^  ip _ j 

Xk  ~  2  2  ~~ 

(  O  Opp)  (  Cc  “  u^'pk) 

(  ^  ~  )^K_-  —  (  ~  i'pk)Xk_  j 

^  k  ~  —  y  " 

(O- Opp)  “  'Upp) 

k  =  2.3,...N 

Evaluate  the  zeros  of  (4.4)  using  the  formulas 

X ,  . s  =  (o  -  oz,  )X  ,,s- ! —  ( a;  —  ccZi 

Y, *k  =  (o-ozl)V,.N  ;  -(i.'-uZI)X,.K-, 

i=  1,2.  ...M 


(4.6) 


(4.7) 


(4.8) 

(4.9) 


This  gives 

G(s)  =  K[XM.N+jYMtN]  (4.10) 
Notice  that  equations  (4.6)  and  (4.7)  fail  when  the  point  (o+ ja>)  is  near 
a  pole.  This  is  overcome  by  having  the  program  check  for,  and  avoid 
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evaluating,  G(s)  at  these  points. 

Figure  4.1  shows  a  computer  flow  chart  of  the  root  locus  algorithm. 
In  the  interests  of  simplicity,  some  of  the  less  significant  steps  are  not 
shown.  For  further  details  see  Appendix  C. 

BLOCK  1  Input  the  zeros  and  poles  of  the  transfer  function. 

BLOCK  2  Display  the  poles  and  zeros.  The  pole  locations  are  marked  with 
an  x  and  the  zeros  are  marked  with  a  circle. 

BLOCK  3  Set  the  loop  counter,  locus,  to  1  and  establish  the  allowable 
locus  deviation  to  0.00001.  The  error  value  was  established  through  a 
trial  and  error  process  and  is  a  compromise  between  locus  accuracy  and 
computation  time. 

BLOCK  4  Check  the  loop  counter  for  being  less  than  or  equal  to  the  total 
number  of  system  loci. 

BLOCK  5  As  explained  in  Chapter  2,  each  locus  starts  at  a  pole  and  ends 

at  either  a  zero  or  asymptotically  at  infinity.  Use  the  angle  criterion 

to  calculate  the  departure  angle  at  the  start  of  the  present  pole.  The 

angle  criterion  is  best  demonstrated  by  example.  Consider  a  system  with 

the  following  transfer  function 

G(s)  = - -  (4.11) 

(s  +  p,)(s*p2)  (s  +  p3)(s-p4) 

As  shown  in  Chapter  2,  the  angle  criterion  is  given  by 

Y,  o,  =£  (angles  of  the  zeros)  -^(angles  of  the  poles)  ^ 

=  (2k  +  l)n 
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Figure  4.1  Root  Locus  Flow  Chart 
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where 


is  chosen  sc  that 


1.-2,  ±  3, 


-77  -  V  —  77 


(4.13) 


(4.14) 


where  the  angles  of  zeros  and  poles  are  shown  in  Figure  4.2.  This  gives 


(\A  =  a3  —  a,  -  n,  -  a4  -  (2k  -  1 )  7T 


(4.15) 


BLOCK  6  To  find  the  point  S,  near  the  locus,  take  an  incremental  step, 
with  a  length  of  As,  from  the  pole  S,  in  the  direction  of  a„.  This  can  be 
seen  in  figure  4.3.  As  is  computed  to  be  approximately  two  percent  of  the 
largest  open-loop  pole  or  zero  coordinate. 

S'c  -  -  As  (co:>  od  -  jsinr\d)  (4.16) 
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Figure  4.2  Angle  criterion  construction 
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BLOCK  7  In  order  for  a  point  s—a± ju>  to  be  on  a  root  locus  it  must  satisfy 


the  equation 


or  using  (4.10) 


or 


1  -  KG(s)  -  0 

1  *  K(X-jV)  -  0 

Y  =  0 


(4.17) 

(4.18) 

(4.19) 

(4.20) 


The  approach  to  using  equations  (4.19)  and  (4.20)  is  to  find  values  of 
s-c± ja>  that  satisfy  (4.19)  and  use  these  in  (4.20)  to  calculate  the 
corresponding  value  of  K.  Since  it  is  unreasonable  to  calculate  values 
that  will  make  (4.19)  exactly  equal  to  zero,  a  small  error  value  can  be 


established  to  minimize  Y. 
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|  Y |  <  error 


(4.21) 


BLOCK  8  If  the  point  SQ  is  far  enough  away  from  the  root  locus  that  (4.21) 
is  not  true,  a  search  at  right  angles  to  the  line  (Sq-SP)  is  done  to  find 
a  point  S„  sufficiently  close  to  satisfy  (4.21).  This  is  shown  in  figure 
4.4. 


>'R  -  -  Lr  l  -  *innd  -jcosad)  (4.22) 
In  order  to  find  the  value  for  L,,  G(s)=K[X+jY]  is  evaluated  at  two  points 
along  the  line  (S„-Sq).  The  two  points  are  determined  by  using  (4.22)  with 
two  different  values  for  L.  They  are 

L,»Asl0‘ 

L?-2AslO  ‘ 

These  two  trial  points  are  applied  to  equations  (4.5)  to  (4.9)  and  the 
resulting  values  Y(L,)  and  Y(L2)  are  interpolated  together  using  (4.23)  to 
minimize  Y. 


L,  -  L;  - 


I.,  ~  L; 
Y(Li)  -  Y(L2) 


Y(Li) 


(4.23) 
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This  process  is  repeated  until  Y(L,)  satisfies  (4.21)  and  at  which  point 
S„  is  assumed  to  be  found. 

With  the  point  S»  found  it  is  possible  to  use  equation  (4.20)  to 
find  the  corresponding  gain,  K.  Furthermore,  if  the  slope  of  the  locus 
at  S«  can  be  determined,  the  sequence  can  be  repeated  until  the  whole  locus 
is  found. 

The  slope  at  S«  can  be  found  from: 


where 


m 


S  ••  M 

V  ~  Pr:  _  V  —  ^Zl 

fl*  =  k~:  rl 

da  v  a~ap^  _  V  a 

1  V'k  i"7"'  r* 


(2.24) 


TpK  i  {<7  (7pk  )'-(*-  -pkf 

and 

ry-  -  (cr  au)2  -  (x~xz  f 

Using  the  two  points  Sr  and  S«,  a  new  point  St-Ot+Jwt  can  be  calculated 
using  the  (usually  correct)  assumption  that  the  slope  of  the  curve  is 
constant  over  a  small  interval  on  the  locus. 


<7t  -  crs 


— -  ~ 


1  ’  m 


1  -  nV 


111 1  {^r  wp) 


1  *  m 


1 

1  -m2 


(-"R  “■» 


(2.25) 

(2.26) 


Equations  (2.25)  and  (2.26)  are  not  used  for  values  of  m  <  0.01  and 
m  >  100,  when  horizontal  and  vertical  slopes  are  assumed  respectively. 


BLOCKS  9 .  10  Breakpoints  in  the  locus  are  treated  in  this  algorithm  as 
the  intersection  of  two  loci,  and  if  the  wrong  path  is  followed  after  such 
an  intersection,  the  calculated  value  of  K  decreases,  instead  of 
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increasing.  If  a  search  step  produces  such  a  decrease,  it  is  regarded  as 
unsuccessful  and  a  new  step  is  taken  at  right  angles  to  the  unsuccessful 
one.  Since  breakpoints  almost  always  occur  at  right  angles,  this  strategy 
is  usually  successful. 

BLOCK  1 1  A  locus  terminates  either  at  a  zero  or  at  infinity  on  an 
asymptote .  Infinity  is  assumed  when  any  of  the  calculated  coordinates 
exceeds  the  screen  coordinates.  An  arrow  is  placed  at  the  point  of  exit, 
pointing  in  the  direction  of  the  exit  slope,  to  indicate  infinity. 

BLOCK  12  The  locus  is  displayed  as  a  set  of  data  points  joined  together 


by  straight  line  segments.  Line  segments  lying  on  the  left-hand  plane, 
in  the  stable  region,  are  displayed  in  green,  while  segments  lying  on  the 
right-hand  plane,  in  the  unstable  region,  are  displayed  in  red. 


CHAPTER  5 


BODE  PLOTS 

JLI _ Introduction  to  Bode_Plots 

This  chapter  serves  as  an  introduction  to  Bode  plots  which  are  used 
in  the  analysis  and  design  of  control  systems.  This  is  a  graphical 
technique  that  plots  a  control  system's  gain  amplitude  and  phase  angle 
response  curves  against  the  input  frequency.  It  is  customary  to  plot  the 
gain  in  decibels  and  the  phase  angle  in  degrees  against  the  common 
logarithm  of  the  input  frequency.  This  is  the  form  first  introduced  by 
H.W.  Bode.  This  discussion  is  not  intended  to  serve  as  a  complete  guide 
to  Bode  plots,  but  rather  an  overview  to  help  the  reader  to  better 
understand  the  computer  techniques  discussed  in  the  following  chapter. 
For  further  details  refer  to  reference  [3]. 

5.2  Bode  Plots 

Consider  the  open-loop  transfer  function  of  a  single-input  single¬ 
output  system  expressed  in  complex  frequency  domain, 

G(s)|^  -  G(j«)  -  A(w)+jB(w)  -  |G(j»)|  /G< jw)  (5.1) 

where  G(jw)  is  a  complex  variable,  while  A(w)  and  B(w)  are  real  variables. 
To  convert  G(jw)  to  polar  form,  that  is  in  magnitude  and  phase  form,  the 
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following  equations  are  used: 

| G(  jw)  |  -  (A2(u>)+B2(u>)  )1/2  (5.2) 

=  /G(  ju>)  -  tan'1  (B(w)/A(w) )  (5.3) 

A  Bode  plot  of  this  system  consists  of  the  gain  and  phase  curves,  defined 
as  20  log | G(  jo?)  |  and  respectively,  versus  the  angular  input 

frequency  w. 

Example :  Consider  a  system  with  the  following  transfer  function 

G(s  )*=!/(  1+rs)  (5.4) 

Gain:  The  logarithmic  gain  is 

20  log | G  |  -  20  log(  l/(  l+(wr  )2)1/2 ) 

—  10  log(l+(wr )2)  (5.5) 

In  order  to  plot  the  gain  versus  input  frequency,  use  the  following  forms 
for  different  ranges  of  the  frequency.  For  w  <  1/r ,  is  approximated  by 

20  log j G |  — - 10  log(l)  -  0  db,  u.  <  1/r  (5.6) 

For  high  frequencies,  w  ►  1/r,  the  gain  (5.5)  becomes 

20  log  j  G |  -  -20  log(wr),  u  ►  1/r  (5.7) 

while  at  w  -  1/r,  the  form  is 

20  log|G|  -  -10  log(2)  -  -3.01  db  (5.8) 

It  is  clear  from  (5.6)  that  the  gain  is  a  constant  0  db  for  low 
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frequencies,  then  it  drops  to  -3.01  db  at  the  corner  or  break  frequency 
Wc  -  1/r  and  finally  for  high  frequencies,  the  plot  drops  linearly  at  a 
slope  of  -20  db/decade ,  that  is  at  a  frequency  range  coi  <  w  <  oo2  such  that 
ui  =  10w, .  Therefore,  because  of  this  linearity,  it  is  very  convenient  to 
plot  the  phase  20  log | G |  versus  log  co.  This  Bode  magnitude  plot  is  shown 
in  figure  5.1. 

Phase :  The  phase  angle  for  (5.4)  is  given  by 

i =  -tan’cor  (5.9) 

which  has  the  following  asymptotic  and  exact  behaviors  for  a  wide  range 
of  frequencies. 

<i(co)  -  0°,  u>  <  1/r  (5.10) 

*(«)  -  -45°,  co  -  1/r  (5.11) 

<Ku)  -  -90°,  w  <  1/r  (5.12) 

The  phase  plot  is  shown  in  figure  5.2. 
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Figure  5.2  Bode  phase  plot 

For  rational  transfer  functions,  it  is  only  necessary  to  be  able  to  plot 
gain  and  phase  curves  of  the  following  types  of  terms: 

1)  Constant  gain  K. 

2)  Poles  and  zeros  at  the  origin  of  the  complex  plane. 

3)  Real  axis  poles  and  zeros. 

4)  Complex  conjugate  pairs  of  poles  and  zeros. 

A  rational  function  can  be  factored  into  terras  of  above  types,  and  the 
individual  gain  and  phase  curves  plotted.  The  complete  gain  curve  is  the 
sum  of  the  individual  component  gain  curves,  and  the  complete  phase  curve 
is  the  sum  of  the  individual  phase  curves.  A  discussion  of  the  four  types 
of  transfer  function  terms  follows. 

1. Constant  Gain  K.  The  logarithmic  gain  is  20  log  K  while  the  angle  is 
0°.  The  magnitude  plots  are  horizontal  lines.  This  is  shown  in  figure 


5.3. 
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2.  Poles  or  Zeros  at  the  origin  of  the  complex  plane  (ju)*'.  A  pole  or 
zero  at  the  origin  have  a  logarithmic  gain  as 

20  1  og j  ( ju> ) * 1 1  =  ±  20  logo  db  (5.13) 
which  is  a  straight  line  on  a  semi -log  paper  with  slope  of 
±20  db/ decade  and  a  horizontal  crossing  at  w  -  1 .  The  phase  angle 
for  this  term  is  ^(w)  -  ±  90°.  Examples  of  poles  and  zeros  at  the 
origin  can  be  seen  in  figures  5.4  and  5.5. 


-20 


•  40 - 


10  100 
Frequency  u:  (rod  sec) 


1000 


Figure  5.3  Bode  plot  with  constant  gain  K-10 


Figure  5.4  Bode  plot  with  one  and  two  poles  at  the  origin 
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Figure  5.5  Bode  plot  with  one  and  two  zeros  at  the  origin 


3.  Poles  or  Zeros  on  Real  Axis  (1+jco/w,  )*  These  terms  have 

log -magnitude , 

20  log |  (l+jw/w,)”|  -  ±10  log(l+ci>/w, )?)  (5.14) 

The  asymptotic  behavior  of  this  plot  begins  at 
±10  logCw/aj,)7  -=  ±20  log(w/w, )  or  a  straight  line  with  a  ±20  db/ decade 
slope.  The  two  asymptotic  lines  (i.e.  0  db  and  ±20  db/dec . )  cross 
each  other  at  the  point  u  -  u,  or  at  the  corner  frequency  (see  Fig. 
5.1).  However,  the  actual  value  of  the  logarithmic  gain  at  u>  -  w, 
is  ±10  log(2)  -  ±3  db  as  demonstrated  with  (5.4)  example.  The  phase 
angle  -  ±  tan'1  (w/u, )  which  begins  at  0°  for  low  frequencies 

(w  <  w,).  reaches  ^(w,)  -  ±tan’(l)  -  ±  45°  and  approaches 
-  ±tan'(®)  -  ±90°,  (See  Fig.  5.2).  Examples  can  be  seen  in 
figures  5.6  and  5.7. 
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Figure  5.6  Bode  plot  with  poles  on  the  real  &xis 


Figure  5.7  Bode  plot  with  zeros  on  the  real  axis 


4.  Complex  Conjugate  Poles  or  Zeros.  [  l+2£(  jw/w2)  +  (  ju/u})2  ]*  ’ .  Let  the 
ratio  w/w2  be  represented  by  quantity  u  and  evaluate  the  logarithmic 
gain  as, 

±20  log|G(jw)|-  ±  10  log((l-w*)*+4£V)  (5.15) 

while  the  phase  angle  is 

*(«,£)  -  tan-,(2fw/(l-w*))  (5.16) 

Once  again  the  asymptotic  behaviors  of  the  above  plots  will  be 
investigated.  When  v  <  1,  the  magnitude  is, 
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±10  log(l)  -  0  db  (5.17) 

and  the  phase  angle  approaches  0°  .  On  the  other  end  of  the 
frequency  scale,  i.e.  for  u  1 ,  the  magnitude  is 

±10  1  og ( i/4 )  -  ±40  log(i/)  (5.18) 

which  results  in  a  straight  line  with  a  slope  of  ±40  db/decade.  The 
phase  angle  4>(u>)  approaches  0°  asymptotically  as  v  <  1  and  reaches 
±180  as  v  >  1 . 

The  frequency  wr  at  which  the  maximum  magnitude  occurs  is 
called  the  resonant  frequency .  Note  that  as  the  damping  ratio  £ 
approaches  zero,  the  resonant  frequency  w,  approaches  the  corner 
frequency  w2.  The  resonant  frequency  is  obtained  by  taking  the 
derivative  of  the  magnitude  of  (5.15)  with  respect  to  u  and  setting 
it  equal  to  zero.  The  resulting  equation  is 

i/-l+2e?  -  0  (5.19) 

or 

«r  -  (1-2£J)VJ,  £<0.707  (5.20) 

The  maximum  value  of  the  magnitude  itself  is 


Mr  -  |G(«r)|  -  (2^(1  -£’)•’,  £<0.707 


(5.21) 
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Figure  5.8  Bode  plot  with  a  complex  pole 


CHAPTER  6 


COMPUTER-AIDED  BODE  PLOTS 

6,1  Introduction  to  Computer-Aided  Bode  Plots 
This  chapter  describes  the  Bode  plot  algorithm  as  shown  in  the 
computer  flow  chart  in  Figure  6.1.  In  the  interest  of  simplicity,  some 
of  the  less  significant  steps  are  not  shown.  For  more  detail,  refer  to 
the  program  listing  in  Appendix  D. 

6.2  Bode  Program  Description 

BLOCK  1  Input  the  open-loop  transfer  function  information,  the  number  of 
poles  and  zeros  and  their  corresponding  complex  roots. 

BLOCK  2  Input  the  open-loop  transfer  function  gain  K  and  the  lower  and 
upper  plotting  frequency  range  values,  omegal  <  w  <  omegau.  If  the  lower 
frequency  value  omegal  is  not  a  power  of  10,  i.e.  not  0.001,  0.01...,  then 
it's  converted  to  the  next  lower  power  of  10.  Likewise,  the  upper 
frequency  is  converted  to  the  next  higher  power  of  ten.  As  an  example, 
if  the  frequency  range  values  were  input  as  omegal-0.05  and  omegau-300, 
then  their  corresponding  converted  values  would  be  initial-0.001  and 
final— 1000.  This  conversion  is  done  in  order  that  the  abscissa  on  the 
output  plot  has  a  starting  and  ending  points  with  powers  of  10. 
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start  ■ 


j  Input  N,  M,  j 
i  poles  and  zeros, 'j' 


ID 


Input  K.  omegnl 
and  omegau  2 

, - — 

_ i _ 

Compute  iterations 

J  =  0  3. 


Figure  6.1  Bode  plot  flow  chart 
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BLOCK  3  Compute  the  number  of  frequency  iterations.  This  is  a  function 
of  the  frequency  range  and  the  number  of  available  horizontal  output 
screen  pixels.  This  insures  that  the  minimum  number  of  iterations  are 
used  that  will  produce  the  most  accurate  plot  with  the  least  amount  of 
computer  run  time.  Next,  set  the  frequency  loop  counter,  j,  equal  to 
zero . 

BLOCK  4  For  each  iteration  of  the  frequency  loop  a  new  frequency  value, 
w,  must  be  computed,  u  is  given  by  Equation  6.1 

(j  -  omegal  +  jAw  (6.1) 

where  A lo  is  the  frequency  step  size. 

Also,  for  each  loop  pass  through  the  frequency  loop,  the  running 
totals  of  the  magnitude  and  phase  and  the  pole  loop  counter,  i,  are  set 
equal  to  zero. 

BLOCK  5  For  each  iteration  of  the  frequency  loop,  each  of  the  M  factored 
poles  are  evaluated  to  adjust  the  magnitude  and  phase  running  totals. 
Block  5  controls  the  poles  loop  and  sends  the  program  flow  to  the  zero 
loop  once  each  of  the  poles  are  evaluated. 

As  described  in  chapter  5,  the  transfer  function  can  be  factored 
into  four  types  of  terms,  they  are 

1)  Constant  gain  K. 

2)  Poles  and  zeros  at  the  origin  of  the  complex  plane. 

3)  Real  axis  poles  and  zeros. 

4)  Complex  poles  and  zeros. 
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BLOCK  6  Checks  for  poles  at  the  origin  and  if  true  sends  the  flow  to 
block  7. 

BLOCK  7  Pole  at  the  origin.  Adjusts  the  magnitude  and  phase  according 
to  the  following  equations. 

magnitude  -  magnitude  -20  log,0(w)  (6.2) 

phase  -  phase  -  n/2  (6.3) 

BLOCK  8  Checks  for  real  axis  poles  and  if  true  sends  flow  to  block  9. 

BLOCK  9  Real  axis  pole.  Adjusts  the  magnitude  and  phase  according  to 
the  following  equations. 

magnitude  -  magnitude  -  20  log10( l+w/wp)  (6.4) 

phase  -  phase  -  tan’(u>/Wp)  (6.5) 

BLOCK  10  Complex  pole.  If  blocks  6  and  8  are  both  false  the  flow  comes 
to  block  10  and  the  magnitude  and  phase  are  adjusted  as  follows. 

magnitude  -  magnitude  - 

20  log10((l-(w/w,)?)7  +4e?(w/Wp)')’/7  (6.6) 

phase  -  phase  -  tan '( 2^ (w/wp)/(l -  (w/wp)’) )  (6.7) 

Once  all  of  the  poles  are  evaluated,  the  loop  counter  i  is  set  equal 
to  1  and  the  flow  goes  to  the  zero  loop  control  block  11. 
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BLOCK  11  For  each  iteration  of  the  frequency  loop  each  of  the  N  factored 
zeros  are  evaluated  to  adjust  the  magnitude  and  phase  running  totals. 
Block  11  controls  the  zero  loop  and  sends  the  program  flow  to  the  constant 
K  block  once  each  of  the  zeros  are  evaluated. 

BLOCK  12  Checks  for  zeros  at  the  origin  and  if  true  sends  the  flow  to 
block  13. 

BLOCK  1 3  Zero  at  the  origin.  Adjusts  the  magnitude  and  phase  according 
to  the  following  equations. 

magnitude  =  magnitude  +20  log,0(w)  (6.8) 

phase  -  phase  +  n/2  (6.9) 

BLOCK  14  Checks  for  real  axis  zeros  and  if  true  sends  flow  to  block  15. 

BLOCK  15  Real  axis  zero.  Adjusts  the  magnitude  and  phase  according  to 
the  following  equations. 

magnitude  -  magnitude  +  20  log10(l+u>/wp)  (6.10) 

phase  -  phase  +  tan1  (u>/wp)  (6.11) 

BLOCK  16  Complex  zero.  If  blocks  12  and  14  are  both  false  the  flow  comes 
to  block  16  and  the  magnitude  and  phase  are  adjusted  as  follows. 

magnitude  -  magnitude  + 

20  log10(  (!-(«/<*)')'  +4 *'(«/«,)')’'’ 


(6.12) 


-42- 


phase  -  phase  +  tan'1  (2£(w/wp)/(  1  -  (w/w,)2) )  (6.13) 

Once  all  of  the  poles  are  evaluated,  the  flow  goes  to  the  constant 
control  block  17. 

BLOCK  17  Once  all  of  the  poles  and  zeros  of  the  transfer  function  are 
evaluated  the  magnitude  is  adjusted  for  the  constant  gain  term  K.  Only 
the  magnitude  is  adjusted  since  the  phase  change  for  a  constant  is  zero. 

magnitude  -  magnitude  +  20  log10(K)  (6.14) 

BLOCK  18  The  current  frequency,  magnitude  and  phase  data  is  stored  once 
all  of  the  terms  of  the  transfer  function  are  evaluated. 

BLOCK  19  This  block  controls  the  frequency  iteration  loop  and  sends  flow 
to  the  output  block  once  the  transfer  function  is  evaluated  for  each  of 
the  step  frequencies.  If  not  finished,  the  iteration  loop  counter  is 
increase  by  one  and  the  loop  continues. 

BLOCK  20  Flow  is  sent  here  once  all  of  the  data  points  that  describe  the 
magnitude  and  phase  curves  are  computed.  Both  curves  are  plotted  together 
sharing  the  same  frequency  abscissa.  The  left-hand  vertical  scale 
represents  the  magnitude  as  measured  in  decibels  and  is  shown  in  magenta. 
The  right-hand  vertical  scale  measures  the  phase  angle  in  degrees  and  is 
shown  in  green.  The  magnitude  curve  is  plotted  with  a  solid  magenta  curve 
and  the  phase  is  plotted  as  a  dotted  green  curve.  The  above  color 
information  only  applies  to  plots  run  on  a  machine  with  an  IBM  compatible 
Enhanced  Graphics  Adapter  (EGA).  Plots  run  on  a  machine  with  an  IBM 
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compatible  Color  Graphics  Adapter  (CGA)  are  only  shown  in  black  and  white 
due  to  the  low  resolution  of  the  CGA  when  plotting  in  color.  The 
magnitude  is  plotted  with  a  solid  white  curve  while  the  phase  is  plotted 
with  a  dotted  white  curve. 

Hard  copy  prints  of  the  plots  are  provided  by  way  of  the  print 
screen  function.  If  a  printout  is  desired,  press  the  P  key  and  the  output 
is  automatically  sent  to  the  printer.  A  example  output  screen  is  shown 
in  figure  6.2. 


CHAPTER  7 


SUMMARY 


7.1  Thesis  Summary 

^  With  the  wide  spread  availability  of  the  personal  computers  to  the 
control  systems  engineer,  is  would  only  make  sense  to  supplement  the  time- 
honored  techniques  of  using  root  locus  and  bode  plots  in  the  analysis  and 
design  of  control  systems  with  one  that  involves  computer-aided  support. 
Using  the  computer  to  do  the  mundane  task  of  generating  the  graphs  allows 
the  designer  to  spend  more  time  analyzing  the  system.  This  could  also 
allow  a  student  in  a  controls  course  the  opportunity  to  learn  more  of  the 
theory  behind  control  systems  analysis  instead  of  spending  all  their  time 
learning  the  tedious  process  of  generating  the  graphs  by  hand.  It  was  the 
purpose  of  this  thesis  to  develope  computer  algorithms  to  automate  this 
task  in  a  easy  to  use  format.  This  software  is  called  the  Control  Systems 
Software  Package  (CSSP). 


-44- 


REFERENCE 


1.  Hale,  Francis  J.:  Introduction  to  Control  System  Analysis  and 

Design.  Prentice-Hall  Inc.,  NJ,  1973. 

2.  Borrie,  John  A.:  Modern  Control  Systems:  A  Manual  of  Design 

Methods .  Prentice -Hall  International  (UK)  Ltd,  1986. 

3.  Jamshidi ,  M.  and  Malek-Zavarei ,  M.  :  Linear  Control  Systems:  A 

Computer-Aided  Approach.  First  Edition.  Pergamon  Books  Ltd. ,  NY, 
1986. 

4.  McDonald,  A.C.  and  Lowe,  H. :  Feedback  and  Control  Systems.  Reston 
Publishing  Company,  Inc.,  Reston,  Virginia,  1981. 

5.  Atkinson,  P.  and  Dalvi ,  V.S.:  An  Improved  Algorithm  for  the 

Automatic _ Determination _ 2f _ Ropts _ Loci  .  Radio  and  Electronic 

Engineer,  41,  1971,  365. 

6.  Hostetter,  G.H.  ,  Savant,  C.J.  Jr.  and  Stefani,  R.T. :  Design  of 

Feedback  Control  Systems .  CBS  College  Publishing,  NY,  1982. 

7.  Gerald,  C.F.  and  Wheatley,  P.O.:  Applied  Numerical  Analysis.  Third 
Edition.  Addison -Wesley  Publishing  Company,  Reading  Massachusetts, 
1985. 


-45- 


APPENDIX  A 


CSSP  MAKE  SOURCE  CODE 
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IHHHHmiHHHHHm***  i'll, 111  I  IHHhmMHHHHHHHHhHHHHh #  ff 

#  This  is  the  make  description  file  used  in  the  development  of  CSSP.  This  # 

#  make  file  was  used  with  the  Microsoft  Make  utility.  In  CSSP  development,  # 

#  the  Make  utility  automatically  updated  the  executable  file  whenever  any  # 

#  of  the  source  or  object  files  were  altered.  t 

#  # 

#  To  use  the  Make  utility  and  this  file,  type  at  the  DOS  prompt  # 

#  * 

#  MAKE  FILENAME  # 

#  # 

#  press  ENTER.  # 

#  # 


cssp . obj :  cssp . c 

cl  ML  /DMSCV4  cssp.c  /c 

locus. obj:  locus . c 

cl  ML  /DMSCV4  locus,  c  /c 

bode. obj:  bode.c 

cl  ML  /DMSCV4  bode.c  /c 

modified. obj :  modified. c 

cl  ML  /DMSCV 4  modified,  c  /c 

f actored . ob j :  factored. c 

cl  ML  /DMSCV4  factored,  c  /c 

coeff.obj:  coeff.c 

cl  ML  /DMSCV4  coeff.c  /c 

video. obj:  video. c 

cl  /AL  /DMSCV 4  video. c  /c 

root. obj:  root.c 

cl  /AL  /DMSCV4  root.c  /c 

printer. obj:  printer. c 

cl  /AL  /DMSCV4  printer. c  /c 

cssp.exe:  cssp. obj  locus. obj  bode. obj  modified. obj  f actored . obj\ 

coeff.obj  video. obj  root. obj  printer. obj 

# 

link  /NOE  /ST: 10000  cssp+locus+bode+modif ied+factored+coeff+\ 
video+root+printer , , , lwin+llibce ; 


appendix  b 


CONTROL  SYSTEMS  SOFTWARE  PACKAGE  SORCE  CODE 
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/*************************  START  OF  CSSP  ROUTINE  **************************/ 


/*  V 
/*  This  is  the  main  program  driver  for  CSSP.  Control  to  each  of  the  */ 
/*  separate  routines,  such  as  the  Bode  plot  routine,  starts  from  and  */ 
/*  and  returns  to  this  program.  */ 
/*  */ 
/*  This  files  also  includes  the  source  for  the  following  funcions:  */ 
/*  */ 
/*  void  tranfunct):  Transfer  function  calling  routine  */ 
/*  sign(double ) :  Check  sign  of  number  */ 
/*  */ 


/************************  START  GLOBAL  VARIABLES  **************************/ 
^include  "variables .h”  /*  Include  variable  list  */ 


struct  mitem  { 
int  r; 
int  c ; 
char  *t; 
int  rv; 

}; 


/*  menu  item  template  */ 
/*  row  V 
/*  col  */ 
/*  text  */ 
/*  return  value  */ 


struct  pmenu  { 

WINDOWPTR  wpsave ; 
int  winopn; 
int  lndx; 
int  fm; 
int  lm ; 

struct  mitem  scrn[25}; 

}  ; 


/*  popup  menu  structure  */ 
/*  a  place  for  the  window  handle  */ 
/*  window  open  flag  */ 

/*  last  index  */ 

/*  first  menu  item  index  V 

/*  last  menu  item  index  */ 

/*  a  bunch  of  menu  items  */ 


int  signtdouble  );  /*  Check  sign  of  number  */ 
/**************************  qjd  GLOBAL  VARIABLES  *«***•********«***********/ 
/*•*******•*************  START  OF  MAIN  CSSP  ROUTINE  *****»****»*»**********/ 

main( ) 


{ 


WINDOWPTR  wl,  w2 ; 

/*  a  few  windows 

*/ 

WINDOWPTR  qpopupt); 

/*  function  returns  WP 

V 

int  i ; 

/*  scratch  integers 

*/ 

int  watrib.batrib; 

/*  scratch  atributes 

»/ 

int  rv; 

/*  for  popup 

V 

int  tranfunc ( ) ; 

/*  Transfer  function  routine 

*/ 

static  struct  pmenu  ml  “  ( 

/*  Main  menu 

*/ 

00,  FALSE,  00, 

06,  9,  { 

01,  02,  "  Main  Menu”,  0, 

02,  02,  0, 

03  ,  00,  - — - 

04,  00,  "  Press  the  desired  menu  number  or  position  the",  0, 

05,  00,  "  highlight  bar  with  the  cursor  keys  and  press  enter.",  0, 

06,  00,  "— - — - 

8,  13,  "[1]  Input  Transfer  Function  Information",  1, 

10,  13,  "(21  Root  Locus  Plot",  2, 

12,  13,  "(31  Bode  Plot",  3, 

14,  13,  "(4)  Quit  CSSP  and  go  to  DOS",  4, 

99,  99,  "",99  } 

); 


set_print_screen( ) ;  /*  Load  print  screen  program  */ 

_aetvideomode(_TE(TC80) ;  /*  Set  video  mode  to  text  */ 

wn_dmode  ( FLASH ) ;  /*  Set  window  speed  to  fast  V 

fort ;  ;  ) 

( 

/* 

*  Set  window  attributes: 

* 

*  border  -  blue/whita  box 

*  window  -  white  background/black  letters 
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*/ 

batrib  -  v_setatr (BLUE, WHITE, 0 , 0 ) ; 
watrib  -  v__setatr  (BLUE,  WHITE,  0 ,  D) ; 

/* 

*  Open  title  window 

V 

wl  =  wn_open(0,0,0 , 78,23 .watrib, batrib) ; 
if ( !wl )  exit( 1) ; 

wn_titla(wl,"  Control  Systems  Software  Package  ", batrib); 

/* 

*  Open  credit  window  at  bottom  middle  row 

V 

w2  =  wn_open( 1000 ,24 ,19,39,1 .watrib , batrib ) ; 
if(!w2)  exit(l); 

wn_printf(w2,"  Copyright  (c)  1988-1989  Carl  F.  Adams  "); 


rv  =  popup  (0,3,6,65, 16,  WHITE«4  |BLACK,  BLUE«4  [WHITE,  Sal ,  TRUE)  ; 


switch  (rv) 

{ 

case  1; 

transfunc( ) ; 
wn_close(w2) ; 
wn_close(wl) ; 
break ; 
case  2: 

wn_close(w2) ; 
wn_close(wl) ; 
output_screen( ) ; 
plot(m,n) ; 
printer ( ) ; 

_setvideomode(_TEXTC8Q) ; 
break ; 
case  3; 

bode(m,n) ; 
printert ) ; 

_setvideomode(_TEXTC80) ; 
break ; 
case  4: 

wn_close(wl) ; 
wn_close(w2 ) ; 
exit(0) ; 

} 


/*  Transfer  function  routine  */ 


/*  Initialize  output  screen  */ 
/*  Root  locus  plot  routine  */ 

/*  Set  video  mode  to  text  */ 

/*  Bode  plot  routine  */ 
/*  Set  video  mode  to  text  */ 


} 

} 

/...**************.**„»,  END  OF  MAIN  CSSP  ROUTINE  *»****■***»****»***».**»»/ 


/.»*•*****.*,***«*,***  START  TRANSFER  FUNCTION  ROUTINE  *********»**********/ 
/* 

*/ 

This  routine  prompts  the  user  for  what  type  of  transfer  function  */ 

is  to  be  plotted.  A  popup  menu  of  the  form  types  is  displayed  and 
control  is  sent  to  the  designated  function  once  its'  function  menu 
number  is  pressed  or  else  the  function  is  highlighted  with  the 
highlight  bar  and  the  ENTDl  key  is  pressed. 


/* 

/* 

/* 

/* 

/* 

/ 


int  transfuncO 

{ 

int  rv; 


/*  for  popup  */ 


static  struct  pmenu  m2  •  { 

00,  FALSE,  00, 

14,  16,  { 

00 >  03,  "  Transfer  Function  Forms", 0, 

01,  03,  0, 

°2-  °3'  "  K(s  +  Cl ) (s  +  C2)(s*2  +  C3s  +  C4)  ... 

03,  03,  ”  Factored  Form  G(s)  ”  _ _ _ - _ 
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04,  03,  " 

05,  03,  0, 

06,  03,  ” 

07,  03,  “  Polynomial  Form 

8,  03,  ” 

9,  03,  0, 

10,  00, 

11,  03,  “  Press  the  desired  function  form  number  or  position  the",  0, 

12,  03,  "  highlight  bar  with  the  cursor  keys  and  press  enter.'',  0, 

13,  00, 

- "  _  o  , 

15,  27,  "[1]  Factored  Form",  1, 

17,  27,  "(21  Polynomial  Form",  2, 

19,  27,  "[3]  Return  to  Main  Menu",  3, 

99,  99,  "",99  } 

}; 

rv  *  popupt  0,1,3,72,21,  WHITE«4  j  BLACK,  BLUE«4  [WHITE,  &m2,TRUE); 

switch  (rv) 

{ 

case  1: 

rv-factored( ) ;  /*  Factored  polynomials  */ 

return ( 1) ; 
case  2: 

coeff();  /*  Non-factored  polynomials  */ 

rv=  ROOTC&pcoef [0] ,&poles[01 [0) ,n) ;  /*  Find  the  pole  roots  */ 

rv-  ROOT(&zcoef  [0]  ,Stzeros[01  CO]  ,m) ;  /*  Find  the  zero  roots  */ 

return(2) ; 
case  3: 

return(3);  /*  Return  to  Main  Menu  */ 

> 

} 

/***********************  END  transfer  function  ROUTINE  ********************/ 

/*************************  START  SIGN  ROUTINE  **************»•*****•******«/ 
/*  */ 
/*  This  routine  checks  for  the  sign  of  the  variable  x.  It  returns  a  V 

/*  value  of  0  for  small  values  of  x,  or  it  returns  a  value  of  1  if  x  */ 

/*  is  equall  to  0,  and  if  the  above  is  not  true  it  returns  the  sign  V 

/*  of  x.  */ 

/+***+**  ****************  ********+******  **************************************/ 
int  sign(double  x) 

{ 

if (fabs(x)<0. 00000000001  &&  x!-0.0) 
return(O) ; 
else  if (x— =0 . 0 ) 
return (1 ) ; 
else 

return(x/fabs(x) ) ; 

) 

/**************************  END  OF  SIGN  ROUTINE  ***************************/ 
/*****************************  end  OF  CSSP  FILE  *+++*+++****++++++***+++*+*+*/ 


Is  +  C5)  (s  +  C6Hs  2  +  C7s  +  C8)  ...  ",0, 

K{Cls“m  +  C2s*(m-1)  +  ...  +  C3s  +  C4)  ",0, 

G(s)  -  -  "  ,  0 , 

{C5s*n  +  C6s*(n-1)  +  . . .  +  C7s  +  C8}  ",0, 


APPENDIX  C 


ROOT  LOCUS  SOURCE  CODE 
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/**************.***.**«.*.  START  OF  LOCUS  ROUTINE  **•****»••**»***•»**•****/ 
/*  */ 
/*  This  is  program  plots  the  root  locus  of  a  given  transfer  function.  */ 
/*  */ 

STASI  GLOBAL  VARIABLES  *******•********••*******/ 


♦include  "variables . h" 


/*  Include  variable  list  */ 


double  sigmap[10], 
sigmaz [ 10] , 
omegap] 10] , 
omegazt 10] , 

Xp.Yp,  /* 


/*  Real  part  of  pole  */ 
/*  Real  part  of  zero  */ 

/*  Imaginary  part  of  pole  */ 

/*  Imaginary  part  of  zero  */ 

X  &  Y  components  of  present  point  */ 


/**«********«*****«*****»**  qjd  GLOBAL  VARIABLES  ***»****************»*****/ 
/*.*..*«*»**«**,.***  START  OF  MAIN  LOCUS  PLOT  ROUTINE  **»***•*****«*»•**«***•/ 

void  plot(int  n.int  m) 

{ 


int 


i  ,  j  ,k  .  1 , 
locus , 

converge_ck , 
repeat , 

repeat_counter-l , 
breakout , 
mfinite_red“0, 
inf inite_green=0 , 
green_length“0 , 
red_length=0 ; 


/*  Loop  counters  */ 
/*  Branch  counter  */ 
/*  Convergence  check  l*True  0*False  */ 
/*  Repeated  roots  */ 
/*  Repeated  roots  counter  V 
/*  Breakout  check  l“True  0“False  */ 
/*  Check  for  red  branch  going  to  ®  */ 
/*  Check  for  green  branch  going  to  ®  */ 
/*  Length  of  green  output  string  */ 
/*  Length  of  red  output  string  */ 


double  theta, 
range , 

Xmax , 

Ymax , 

Y_min , 
alpha [9] , 
beta[9] , 
pi-3.1*159265A, 
XI , X2 , X3 , 
Y1,Y2.Y3, 
midx.midy , 
GREEN_max”0 . 0 , 
GREEN_min»0 . 0 , 
RED_max=0 . 0 , 
RED_min“0 . 0 , 
x, 

K, 

K1.K2.K3, 

Ypl,Yp2,Yp3, 

LI , L2 , L3 , 
error“0 . 00001 , 
XX, YY; 


/*  Present  branch  angle  */ 
/*  Allowable  range  for  convergence  to  zero  */ 
/*  Maximum  allowable  X  value  */ 
/*  Maximum  allowable  Y  value  */ 
/*  Minimum  allowable  Y  value  */ 
/*  Array  of  pole  angles  */ 
/*  Array  of  zero  angles  */ 

/*  Sequential  X  values  */ 
/*  Sequential  Y  values  */ 
/*  Middle  index  values  */ 
/*  Maximum  green  branch  value  */ 
/*  Minimum  green  branch  value  */ 
/*  Maximum  red  branch  value  */ 
/*  Minimum  red  branch  value  V 
/*  Present  branch  slope  */ 
/*  Present  gain  value  */ 
/*  Sequential  gam  values  */ 
/*  Sequential  deviation  values  */ 
/*  Sequential  length  modifiers  */ 
/*  Maximum  allowable  branch  deveiation  V 
/*  Temporary  position  values  */ 


char  green_string[20] , 
red_string(20] ; 


/*  Red  branch  string  buffer  V 
/*  Green  branch  string  buffer  */ 


double  slopednt,  int, double  .double) ;  /*  Slope  routine  */ 
double  Y(int  , int  .double  .double  );  /*  Position  equation  routine  */ 
double  gaindnt  ,int  .double  .double  );  /*  Gain  routine  */ 
int  convergednt  ,  double,  double,  double  ) ;  /*  Converged  ?  */ 
void  sortdnt  );  /*  sort  array  values  (lowest  ->  highest)  */ 


delte_s«160 . 0/max/ 75 ; 
range“delta_s ; 

Ll-delta  s*0.0001; 
L2“delta”s*0 . 0002 ; 
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X_max-( vc .numxpixels-4 ) /2/max; 

Y_max-X_max*3/4 ; 

Y_min— (Y_max-32/max*<i/3) ; 

sort(m) , 

for (i-1 ; i<«m ; i++ )  /*  Load  real  &  imaginary  pole  values  */ 

( 

sigmap] 1 1 “poles [ i - 1 ] [15 ; 
if(fabs(poles[i-l][2] )<error ) 
omegap] l  ]-0  ; 
else 

omegap] i ) -poles ( i-ll (2) ; 

} 

for ( i-l ; i<-n ; i++ )  /*  Load  real  &  imaginary  zero  values  */ 

{ 

sigmaz [i]-zeros(i-l)[ll; 
if(fabs(zeros[i-l) [2] )<error) 
omegaz] l ]-0 ; 
else 

omegaz [i]=zeros[i-l)[2); 

} 

for( loeus-1 ; locus<-m; locus++ )  /*  Start  of  main  branch  loop  */ 

{ 

k-0  ; 

BEGIN:  Xl-sigmap [ locus ) ; 

.'l=omegap(  locus]  ; 

_moveto(Xl*max, -Yl*aspect_ratio*max) ; 

theta-0 . 0 ; 
repeat-0 ; 

for(i-l;  i<-m;  i++ )  /*  Compute  pole  departure  angle  component  */ 

{ 

x-sign(sigmap[ 1 ] ) ; 

iff  fabstXl-sigmap] i ] )<-range  &&  fabs(Yl-omegap[i J )<-range  ) 

< 

repeat++ ; 
alpha] 1 J-0 . 0; 

> 

else 

alpha] 1 ]-atan2( (Yl-omegap(i) ) , (Xl-sigmap[ i  ] ) ) ; 

theta-theta-alpha [  i  ] ; 

) 

for] l-l; i<-n; i++)  /*  Compute  zero  departure  angle  component  */ 

{ 

if  (Xl--sigmaz[i]St&Yl--omegaz!i] ) 
beta] i ]-0 . 0  ; 
else 
{ 

beta] 1 ]-atan2( ( Y 1 -omegaz [ i ] ) , (XI -sigmaz [ 1 ] ) ) ; 

) 

theta-theta+beta [ i 1 ; 

) 

if(repeat>l)  /*  Check  for  repeated  roots  */ 

< 

theta“(theta+pi )/repeat  +  2*pi*(repeat_counter++- l)/repeat; 
i f ( repeat_counter>repeat ) 
repeat_counter-l ; 

} 

else 

theta-theta-pi ; 

X2-X1+1 . 5*delta_s*cos( theta ) ; 

Y2-Y1+1. 5*delta  s*sin( theta ) ; 


Yp2-Y(m,n,X2,Y2); 

whila(fabs(Yp2)>error) 

< 


/*  Check  for  branch  convergence  */ 
/*  and  itterate  as  needed.  */ 
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XX-X2; 

YY*Y2; 

X2-XX-Ll*sin(theta); 

Y2“YY+Ll*cos (theta ) ; 

Ypl=Y (id , n ,X2 ,  Y2) ; 

X2“XX-L2*sin(theta) ; 

Y2-YY+L2*cos( theta); 

Yp2«Y(m,n,X2,Y2) ; 

L3=Ll-(Ll-L2)/(Ypl-Yp2)*Ypl; 
if ( f abs(L3 )>delta_s/4 , 0 ) 

L3»sign(L3 )*delta_s/4 . 0 ; 

X2“XX- L3*s in ( theta) ; 

Y2=YY+L3*cos(theta) ; 

Yp2=Y(m,n,X2,Y2); 
theta=atan( ( Y2-Y1 )/ (X2-X1 ) ) ; 

} 

if(X2>0.0)  /*  Plot  in  red  if  branch  is  unstable  */ 

_setcolor(RED) ; 

else  /*  else  it  is  stable  and  plot  in  green  */ 

_setcolor(GREEN ) ; 

_lineto(X2*max, -Y2*aspect_ratio*max ) ; 

K2“-1.0/Xp; 

if (fabs(Y2)>delta_s/2) 
breakout=l; 
else 

break out=0 ; 


l*j“0 ; 

while( j<-800 )  /*  Compute  and  plot  locus  branch  */ 

{ 

START:  j++; 

x«slope(n,m,X2,Y2) ; 

if (fabs(x)<0 . 01 )  /*  Check  for  horizontal  slope  */ 

{ 

X3  =  X2  +  sign(X2-Xl)*delta_s*2; 

Y3=Y2 : 

theta*sign(X2-Xl )*pi ; 

) 

else 

if (fabs(x)>100 . 0&&fabs(Y2)>delta_s)  /‘Check  for  vertical  slope*/ 

{ 

X3-X2; 

Y3  -  Y2  +  sign(Y2-Yl)*delta_s*2 ; 
theta=sign(Y2-Yl)*pi*0. 5; 

) 

else 

if (fabs(x)>100.0&&j>l) 

{ 

X3-X2 ; 

Y3  -  Y2  -  sign(x)*delta_s*2; 
theta-sign(Y2-Yl)*pi*0.5; 

) 

else 

< 

X3  -  X2  +  ( 1 . 0-x*x)/( 1 . 0+x*x)*sign(X2-Xl)*delta_s  + 

(2 . 0*x)/ ( 1 . 0+x*x)*»ign( Y2-Y1 )*delta_s ; 

Y3  -  Y2  +  (2. 0*x)/( 1 . 0+x*x)*aign(X2-Xl)*delta_s  + 

( 1 . 0+x*x)/ ( 1 . 0+x*x)*»ign( Y2-Y1 )*delta_s; 
theta-atan( ( Y3-Y2 ) / (X3-X2 ) ) ; 

) 

/*  Check  for  break-in  point  */ 

if  (aign(Y3)  !*aign(Yl  )M>fabs  ( Y3  )>error&&f  abs(  Y1  )>error&itbreakout— 1) 

{ 

X2-X2-Y2*(X3-X2)/(Y3-Y2) ; 

X3-X2+«ign(Y3)*delta_e*2; 
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Y2=Y3=sign(Yl-Y2) ‘error ; 

Yp2“Y(m,n,X2,Y2) ; 

K2— 1.0/Xp; 

> 

Y3:  Yp3“Y(m,n,X3 , Y3 ) ;  /*  Again,  check  for  branch  convergence  */ 

while(fabs(Yp3)>error)  /*  and  adjust  as  needed,  */ 

{ 

XX=X3 ; 

YY=Y3 ; 

X3-XX-L1* sin (theta) ; 

Y3=YY-H.l*cos(  theta) ; 

Ypl=Y(m,n,X3,Y3); 

X3=XX-L2*sin( theta) ; 

Y3-YY+L2*cos(theta) ; 

Yp2=Y(m,n,X3,Y3); 

L3-Ll-(Ll-L2)/(Ypl-Yp2)*Ypl; 

if(fabs(L3)>delta_s/1.0) 

L3“sign(L3 )*delta_s/4 . 0 ; 

X3“XX~L3*sin( theta) ; 

Y3=YY+L3*cos ( theta ) ; 

Yp3-Y(m,n,X3,Y3); 
j++; 

if ( j>800)  goto  NEXT; 

} 

K3=- 1.0/Xp; 

if (1++>10) 

{ 

1-0: 

midx“X3 ; 
midy=Y3 ; 

) 

if(K3<K2)  /*  Check  to  make  sure  the  branch  is  »/ 

{  /*  going  in  the  correct  direction.  */ 

convergeck-converge (n , range , X2 , Y2 ) j 
if (convergeck  >  0) 
goto  NEXT ; 

Ypl“Y(m,n,Xl,Yl); 

Yp2“Y(m,n,X2,Y2); 

Yp3=Y(m,n,X3,Y3) ; 

if (sign (Ypl) ! -sign (Yp2) ) 

{ 

X2»X2- (X2-X1 )/( Yp2-Ypl >*Yp2 ; 

X1“X2 ; 

Y2“Y3“sign(Yl-Y2)* error ; 
break out=l ; 

} 

else 

{ 

X2-X2 - ( X3 -X2 ) / ( Yp3 - Yp2 ) * Yp2 ; 

X1-X2; 

Y2“Y3“sign(Yl*Y2) ‘error ; 
breakout“l ; 

> 

CHECK:  Yp2-Y(m,n,X2,Y2) ; 

while (fabs(Yp2)>error) 

{ 

XX-X2; 

YY-Y2; 

X2-XX+L1 ; 

Ypl“Y(m,n,X2,Y2>; 

X2-XX-L1; 

Yp3-Y(m,n,X2,Y2); 

L3“Ll*Ypl/(Yp2-Ypl); 


/*  Check  for  non- convergence,  */ 
/*  is  it  going  no  where?  */ 
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if (fabs(L3)>delta_s/4 . 0 ) 

L3“sign(L3)*delta_s/4 .0 ; 

X2=XX+L3 ; 

Yp2=Y(m,n,X2,Y2) ; 
theta*atan( (Y2-YY)/(X2-XX) ) ; 

if ( j>800)  goto  NEXT; 

} 

K2=-1.0/Xp; 

Y2=Y2-sign(X3-Xl)*delta_s; 
goto  START; 

) 

>e 

{ 

if  CX3>0 . Q ) 

{ 

forcolor2=R£D  ; 

_setcolor(forcolor2) ; 

if (RED_max==0 . 0  &&  RED  min“0 . 0  &&  j!-l) 
RED_min«K2-X2*(K3-K2)/(X3-X2) ; 
if  (K3>RED_max )  RED_max”K3 ; 
if (K3<R£D_min)  RED_min=K3; 

) 

else 

{ 

forcolor2=GREEN ; 

_setcolor(forcolor2) ; 

if  (GREEN_max=0 . 0  &&  GREEN_min—0 . 0  &&  j!»l) 

GREEN _min«K2 -X2* (K3 -K2 ) / (X3-X2 ) ; 
if (K3 >GREEN  jnax )  GREEN_max«K3 ; 
i f ( K3 <GREEN_mi n )  GREEN  min-K3; 

) 

converge_ck"converge(n, range, XI ,Y1) ; 
if (convergeck  >  0) 

{ 

_lineto(Xl*max, -Yl*aspect_ratio*max) ; 
goto  NEXT ; 

} 

if (f abs(X3 )>X_max) 

{ 

if Cforcolor2“”GREEN)  inf inite_green-l ; 
if ( forcolor2=*RED )  infinite_red-l ; 
if ( j<10) 

{ 

tnidx-X2 ; 
midy*Y2; 

} 

if  ( (Y3-midy  )“0 . 0  &&  X3<0.0) 
theta=pi ; 

else  if  ( (Y3-midy  )“0 . 0  ) 
theta«0 . 0 ; 

else 

theta*atan2( (Y3-midy ) , (X3-midx) ) ; 
arrow(X3*max , - Y3*aspec t_ratio*max , theta , forcolor2 , aspect_ratio ) 
goto  NEXT ; 

) 

if (Y3>*Y_raax) 

{ 

if  ( forcolor2”<!REEN )  inf inite_green-l ; 
if (forcolor2“RED }  infinite_red“l; 
if ( j<10) 

{ 

midx-Xl ; 
midy«Yl ; 

) 

if((X3-midx)— 0.0) 

th»ta"sign(Y3)*pi/2. 0; 

•Is* 

thata*atan2( (Y2-midy ) , (X2-midx) ) ; 
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arrow(X3*max, -Y3*aspect_ratio*max, theta, forcolor2, aapect_ratio) ; 
goto  NEXT; 

} 

if (Y3<“Y_min) 

{ 

if Cforcolor2“^3REEN)  inf inite_green-l ; 
if  (forcolor2—RED)  infinite_red=l ; 
if ( j<10) 

{ 

midx-Xl; 
midy=Yl ; 

} 

if ( (X3-midx) ! -0 . 0 ) 

theta*atan2( (Y3-midy ) , (X3-midx) ) ; 
else 

theta-signC Y3 )*pi/2 . 0 ; 

arrow(X3*max , -Y3*aspect_ratio*max , theta , f orcolor2 , aspect_ratio ) ; 
goto  NEXT; 

} 

_lineto(X3*max, -Y3*aspect_ratio*max) ; 

K1=K2 ; 

K2*K3 ; 

X1=X2 ; 

X2-X3 ; 

Y1-Y2; 

Y2-Y3 ; 

Ypl=Yp2; 

Yp2”Yp3 ; 

) 

} 

NEXT: 

if ( j>800 ) 

{ 

delta  s-delta  s/2; 
goto  BEGIN; 
k++; 

} 

if (k>5 ) 
break ; 


if  (infinite_red—l  &&  infinite_green!-l  &&  REDjnin ! «0 . 0  ) 

{ 

_settextposition(vc .numtextrows-1 ,2) ; 

greenlength-sprintf Cgreen_string, "Stable  Gain  Range:  0<K<t-.3g" 

, REDmin ) ; 

_outtext(green_string) ; 

_settextposition(vc .numtextrows , 2) ; 

red_length-sprintf Credstring , "Unstable  Gain  Range:  X* . 3g<K<\35V 
,RED_min) ; 

_outtext(red_string) ; 

) 

else  ifCinfinite  red”l  &&  infinite_green!“l  &&  GREEN_min>0.0) 

{ 

green_length"sprintf(green_string , "Stable  Gain  Range:  X- . 3g<K<X- . 3g" 

, GREEN_min , GR£EN_max ) ; 

_settextposition(vc .numtextrows-1 ,2) ; 

_outtext(green_string) ; 

red_length*sprintf (red_string, "Unstable  Gain  Range:  0<K<X-.3g  X- . 3g<K<\354" 
, GREENjnin , GREEN_max ) ; 

_settextposition(vc .numtextrows ,2) ; 

_outtext(red  string); 

) 

else  ifCinfinite  green—1  &&  inf inite_red!"l  &&  RED_min>0.0) 

{ 

green  length~sprintf (green  string , "Stable  Gain  Range:  0<K<X-.3g 

X-.3g<K<\35*" 

,  RZDjnin ,  REDjnax ) ; 
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_settextposition(vc  .numtextrows-1 , 2) ; 

_outtext(green_string) ; 

red_length*sprintf (red_string , "Unstable  Gain  Range:  X- . 3g<K<X- . 3g" 
,RED_min,RED_max) ; 

_settextposition(vc .numtextrows , 2) ; 

_outtext(red_string) ; 

} 

else  i£(infinite_green“l  &&  in£inite_red“l  &&  RED  min>0.0) 

{ 

green_length“sprintf (green_string , "Stable  Gain  Range:  0<K<X-.3g" 

,RED_min) ; 

_settextposition(vc . numtextrows- 1 , 2 ) ; 

_outtext(green_string) ; 

red_length”sprint£(red_string, “Unstable  Gain  Range:  X- . 3g<K<\354” , 

RED_min) ; 

_ settextposi t ion ( vc . numtextrows , 2) ; 

_outtext(red_string) ; 

) 

else  i£(infinite_green=“l  &&  in£inite_red!“l  &&  RED_min“0.0  &&  RED  max'.*0.0) 

{ 

green_length“sprint£(green_string, "Stable  Gain  Range:  X- . 3g<K<\354" 

,RED_max) ; 

_settextposition(vc . numtextrows- 1,2); 

_outtext(green_string) ; 

redlength-sprintf (red_string, "Unstable  Gain  Range:  0<K<X-.3g" 

, REDmax ) ; 

^settextposi t ion ( vc . numtextrows , 2 ) ; 

_outtext(red_string) ; 

T 

else  if (infinite_green“l  &&  inf inite_red! “1  &&  RED  max ! “0 . 0 ) 

{ 

green_length»sprintf(green_string, "Stable  Gain  Range:  X- .  3g<K<\354" 

,  RED  max ) ; 

_settextposition(vc .numtextrows- 1,2) ; 
outtext(greenstring) ; 

red_length“sprintf (red_string, "Unstable  Gain  Range:  0<K<X- . 3g" .REDmax) ; 
_settextposition(vc . numtextrows , 2) ; 
outtext(redstring) ; 

> 

else  if (infinitegreen!*!  &&  inf inite_red!“l  &&  RED  min ! "0 . 0 ) 

{ 

green_length“sprintf (green_string , "Stable  Gain  Range:  0<K<X-.3g" 

, (GREEN_max+RED_min)/2) ; 

_aettextposition(vc . numtextrows-1 ,2) ; 

_outtext(green_string ) ; 

red_length“sprintf(red_string, "Unstable  Gain  Range:  X- . 3g<K<\354" 

, (GR£EN_max+RED_min)/2,RED_max) ; 

_settextposition(vc .numtextrows , 2) ; 

_outtext(red_string) ; 

} 

else  if(infinite_green!“l  &&  in£inite_red!“l  4A  REDjnax !  »0 . 0 ) 

{ 

green_length-sprintf(green_string, "Stable  Gain  Range:  X- . 3g<K<X- . 3g" , 

GREEHjnin, (GREEN_max+RED_min)/2) ; 

_settextposition(vc .numtextrows-1, 2) ; 

_outtext(green_string) ; 

red_length-sprintf(red_string, "Unstable  Gain  Range:  X- . 3g<K<l- . 3g" 

, (GREENjnax+RED  _min7 /2 , RED_max ) ; 

_sattextposition(vc .numtextrows , 2) ; 
outtext(red_string) ; 

} 

else  if  (RED_max~0 . 0) 

< 

green_length“sprintf (green_string , “Stable  Gain  Range:  0<K<\354"); 
_settextposition(vc .numtextrows- 1 ,2) ; 
outtexttgreen  string); 

) 
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else 

{ 

red_length=sprintf ( red_string , "Unstable  Gain  Range:  0<K<I\3S4"); 
_settextposition(vc . numtextrows- 1 ,2) ; 

_outtext(red_string) ; 

} 

_setcolor ( forcolor 1 ) ; 

_setlogorg  (  0,  0  ); 

if (green_length>=red_length ) 

_r ec tangle (_GBORDER, 0 , min, (green_length+2)*8 , (vc .numypixels-1) ) ; 
else 

_rectangle(_GBORDER,  0  ,mm,  (red_length+2)*8,  (vc  .  numypixels-1) ) ; 


> 

/**************.*.**  START  OF  MAIN  LOCUS  PLOT  ROUTINE  ***********************/ 
/***************************  EQUATION  Y  PROGRAM  ***************************/ 


/*  */ 

/*  This  routine  computes  the  equation  for  Y  as  given  in  the  text,  */ 

/*  "Modern  Control  Systems,  A  Manual  of  Design  Methods”,  by  */ 

/*  John  A.  Borrie,  1986  Prentice-Hall  International  (UK).  */ 

/*  The  equation  Y  is  given  in  equations  (2.6)  to  (2.10)  pages  */ 

/*  103  to  105.  */ 

/*  */ 


/***********«**..*.***  START  OF  EQUATION  FOR  Y  ROUTINE  ********************/ 

double  Y(int  m,  int  n,  double  X,  double  Y) 

{ 

double  sigma_diff, 
omega_di f f 
Xpnew ; 

int  i , 
k ; 

sigma_diff=X-sigroapt 1 ) ; 
omega_diff«Y-omegap[ 1) ; 

Xp»sigma_diff/(sigma_diff*sigma_diff  +  omegadif f*omega_dif f ) ; 
Yp=-omega_diff/(sigma_diff*sigma_diff  +  omegadif f*omega_diff ) ; 

for(k-=2;k<^n;k++) 

{ 

sigma_dif f-X-sigmap [k ] ; 
omega_diff=Y-omegap[k] : 

Xp_new“(sigma_diff*Xp  +  omega_diff'‘ fp)/(sigma_diff*sigma_diff 
+  omega_diff*omega_diff ); 

Yp“(sigma_diff*Yp  -  omega_diff*Xp)/(sigma_diff*sigma_diff 
+  otnega_diff*omega_diff ) ; 

Xp*=Xp  new; 

} 

for(i-l; i<”n; i++) 

{ 

sigma_diff“X-sigmaz[i! ; 
omega_diff“Y-omegaz [ il ; 

Xp_new”sigma_diff*Xp  -  omega_dif f*Yp; 

Yp«sigma_diff*Yp  +  omega_dif f*Xp; 

Xp“Xp  new; 

) 

retum(Yp) ; 

} 

/»**.**.******«******«.♦  qjD  OF  EQUATION  Y  PROGRAM  ****•***********•*******/ 
/**••*••****•****•*•••*••••••**  SLOPE  PROGRAM  ****♦»***•»*****»•«***•*•*•*•/ 


/* 

*/ 

/* 

This  routine  computes  the  slope  m  as  give 

the  text, 

*/ 

/* 

"Modern  Control  Systems,  A  Manual  of  Desi* 

-  rhods" ,  by 

*/ 

/* 

John  A.  Borrie,  1986  Prentice-Hall  Interne 

;al  (UK). 

*/ 

/* 

The  slope  equation  is  given  in  equation  (2 

lb)  on  page  108. 

*/ 
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/* 

START  OF  SLOPE  ROUTINE  »*»******»• 
double  slopednt  zeros,  lnt  poles,  double  sigma,  double  omega) 


»/ 


/ 


( 

doub 1  e  m , 

rpk-0.0, 
rzi-0 . 0 , 
omega  p“0  0 , 
omega  z^O . U, 
sigma  p-t)  U, 
sigma  z-(J.D, 

lnt  1  , 

k  , 

for (k- 1 ;k<*poles ;k++ ) 

{ 

rpk*  ( omega-omegapl k 1 )* (omega- omegapf  k ) ) ♦ (sigma- s l gmap 1 k ) )*( sigma- sigmapl k 1 ) 
omega  p'omega  p+ ( omega-omegapl k I )/rpk ; 
sigma  p-sigma  p+ ( sigma- s lgmaplk ) )/rpk ; 

) 

tort  1*1 ; i<*zeros; i++) 

( 

rzi » (omega- omega z ( i ! )* ( omega-omegaz I l ) )* ( sigma-sigmazl l ) )* ( sigma- sigmaz I i ] ) 
omega  z-omega  z+ (omega -omegaz I  1 ] ) /rzl ; 
sigma  z“a lgma  z+ (sigma- sigmaz { i  )  )/rz l ; 

) 

m-lomega  p-omega  z  )/( sigma  p- sigma  z ) ; 

return(m) ; 

} 

Q(D  OF  SLOPE  PROGRAM  **•*»***»*•*•**••»********/ 


CONVERGE  PROGRAM  ••**•*»»******•******•****••/ 


/*  «/ 

/*  This  routine  determines  If  the  (X,Y)  has  converged  to  any  */ 

/*  of  the  n  zeros  It  checks  for  convergence  of  plus  or  minus  */ 

/*  2*range  where  range  Is  given  in  the  main  program  If  it  is  */ 

/*  within  the  limit  the  function  returns  1  for  true,  else  it  */ 

/*  returns  0  for  faluae  */ 

/*  */ 

START  OF  CONVERGE  ROUTINE  »»»************»*****•**/ 
int  ronvergetint.  n,  double  range,  double  X,  double  Y) 

( 

lnt  1  ; 

double  crange*2*range , 


for( 1-1 ; i<-n ;  1++ ) 

{ 

if(  f abs (X- si gmaz 1 1 1 )<“crange  UU  tabs ( Y-omegaz d ] )t_crange ) 
returnt 1 ) . 

> 

return(0 ) ; 

) 

END  0F  CONVERGE  ROUTINE  *•*»•*»»**•«*»«**********/ 


SORT  PROGRAM  eeeeeeeeeeeeeeeeeee 

/* 

J *  This  routine  sorts  the  n  poles  in  lowest  to  highest  order 

/* 

START  OF  SORT  ROUTINE  •**»»****•••*•* 

void  lortdnl  n) 

( 


double  tempi  2 ] ; 

lnt  gap, 

1 , 

J; 


,»/ 

*/ 

«/ 

*/ 

'*/ 


for(gap-n/2;gap>0 ,gap/-2) 
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{ 

for (i-gap; i<n; i++ ) 

{ 

for ( j=i-gap ; j>-0&&poles [j 1 [l]>poles[ j+gap] [1] ; j"“gap) 

{ 

tempi  01 “poles! j ] [ 1 J ; 
tempi  11 -poles! j ] 12] ; 
poles lj  1  [1] -poles! j+gap] 11)  ; 
poles [j ] 12 1 -poles [ j+gap 1 121 ; 
poleslj+gapl tl]-temp[01 ; 
poles! j+gap] [21-temp! 11 ; 

) 

} 

) 

} 

/******rt*ft»r*****wvt*rt-******  enq  QP  SORT  PROGRAM  a**************************/ 

/***********************  START  OF  LOCUS  OUTPUT  SCREEN  *************•*******/ 

/*  */ 

/»  This  is  program  plots  the  root  locus  of  a  given  transfer  function.  */ 

/*  */ 

/******************-*-*-**  START  OF  LOCUS  ROUTINE  *************************/ 
void  output_screen( ) 

{ 

int  l , 
j  . 
x . 

xlow, 
xhigh, 
xcenter , 
ycenter , 
step, 

textlength, 
position! 10) [21 ; 

float  index, 
textx, 
texty , 
textmax. 
fabsl , 
fabs2 ; 

char  indexstring! 10) ; 

if(!set  video_mode( ) ) 

{ 

printf ("NnWarning!  This  program  doesn't  support  this  machine's  video  card"); 
exit(0) ; 

} 

getvideoconf ig(&vc ) ; 

aspect_ratio  -  (float)  (8.0  *  vc .numypixels )/( 6 . 0  *  vc .numxpixels) ; 
xcenter  -  vc .numxpixels  /  2  -  1; 
ycenter  -  vc . numypixels  /  2  -  1; 

_setcolor(GREEN ) ; 
x-vc . numxpixels/ «0 ; 

_moveto( 1,1); 

_lineto(x,x*aspect_ratio) ; 

_moveto( 1 , x*aspect_ratio ) . 

_lineto(x, 1) ; 

greenpole-(char  »)malloc((unsigned  int)  _imagesize(0 , 0 , x, 
x*aapect_ratlo) ) ; 

_getimage(0,0,x,x»aspect_ratio,  greenpole) ; 

xlow-50-x/2 ; 
xhigh-50+x/2; 

_moveto( 50 , 50) ; 

_ellipse(_GB0RDJ31,xlo«,xlow*»spect_ratio,xhigh,xhigh*esp#ct_retio) ; 
greenzero-(char  * Imalloc ( (unsigned  int)  _in>agesize(xlow,xlow+ 
aspect_ratio-l , xhigh, xhlgh»espect_ratio+l ) ) ; 

_getimage(xlow,xlow*aspect_ratlo-l,xhigh,xhigh*aapect_ratio+l , greenzero) ; 


_setcolor(RED) ; 
x“vc .numxpixels/40 ; 

_moveto( 1 ,  1 )  ; 

_lmeto(x , x*aspect_ratio ) ; 

_moveto( 1 , x*aspect_ratio) ; 

_lineto(x, 1) ; 

redpole=(char  *)malloc( (unsigned  int)  _imagesize (0 , 0 , x , x*aspect_ratio) ) ; 
_getimage(0 , 0 , x, x*aspect_ratio ,  redpole) ; 

xlow=50-x/2; 
xhigh=50+x/2 ; 

_moveto(50 ,  50) ; 

_ellipse(_GBORDER, xlow, xlow* aspect_r a tio, xhigh, xhigh*aspect_r at 10) ; 
redzero=(char  * )malloc ( (unsigned  int)  _imagesize (xlow, xlow*aspect_ratio- 1 , 
xhigh , xhigh*aspect_ratio+l ) ) ; 

_get image (xlow, xlow*aspect_ratio- 1 , xhigh , xhigh*aspect_ratio+l , redzero ) ; 

_clearscreen(_GCLEARSCREEN ) ; 

_setcolor (forcolorl ) ; 

_settextcolor(forcolorl) ; 

_settextposition( (vc .numtextrows-1) , 78-maxchar+skipn) ; 

_outtext(nstring) ; 

_settextposition( (vc . numtextrows ) , 78-maxchar+skipd) ; 

_outtext(dstring ) ; 

_moveto(630-maxchar*(vc .numxpixels/vc , numtextcols ) , 460*aspect_ratio) ; 
_lineto(630.460*aspect_ratio) : 
min=“vc  ,numypixels-43*aspect_ratio; 

_rectangle(_GBORDER, 0 ,min, (vc .numxpixels-1) , (vc . numypixels-1 ) ) ; 

_moveto( (vc .numxpixels-2- (maxchar+2)*vc .numxpixels/vc .numtextcols ) ,min) ; 
_lineto( (vc .numxpixels-2- (maxchar+2)*vc . numxpixels/vc .numtextcols) , 

(vc .numypixels-1 ) ) ; 

_settextposition( 1 , vc .numtextcols/2+2) ; 

_outtext( " jw" ) ; 

_settextposition(vc .numtextrows/2, 79) ; 

_outtext(  ”\345'' ) ; 

_rec tangle (  GBQRDER, 0 , 0 , (vc .numxpixels-1) , (vc .numypixels-1) ) ; 
arrow(xcenter ,2,1.5708, forcolorl , aspectratio ) ; 
arrow(vc .numxpixels-2 .ycenter , 0 , forcolorl , aspect_ratio) ; 
_setlinestyle(0x0101) ; 

_moveto(0 .ycenter ) ; 

_lineto(vc . numxpixels .ycenter ) ; 

_moveto(xcenter , 0 ) ; 

_lineto(xcenter ,min) ; 

_setlinestyle(OxFFFF) ; 

_moveto(0,min) ; 

_lineto( vc . numxpixels -2 , min ) ; 
max^O . 0 ; 

for ( i=0 ; i<=n-l ; i++) 

{ 

f abs 1= tabs (poles ( i ] 1 1 ] ) ; 
fabs2“f abs (poles! i ) (2) ) ; 
if (fabsl>fabs(max) )  max“fabsl; 
if ( f abs2>f abs {max ) )  max”£abs2; 

> 

for ( i“0 ; i<^n-l ; i++) 

{ 

f abs 1-f abs ( zeros ( i ] [ 1] ) ; 
fabs2“fabs  (zeros  !  l ] [2! ) ; 
if ( fabsl>fabs(max) )  max~fabsl; 
if (f abs2>fabs(max) )  max*fabs2; 

} 

taxtmax~2*max ; 

/*  Set  up  coord,  tic  marks  */ 
textx*vc . numtextrows/2+2 ; 
for(i-l;i<-9;i++) 

{ 

Step“64*i ; 
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index=(i-5)*textmax/5; 

textlength-sprintf( indexstring , "I- . 3g" , index) ; 
texty*i*vc .nun)textcols/10-textiength/2+l ; 
if  (i“5) 
continue ; 
else 
{ 

_settextposition(textx, texty ) ; 

_outtext(indexstring) ; 

__moveto(  step  .ycenter ) ; 
lineto( step , (ycenter+10*aspec t_ratio) ) ; 

} 

) 

texty»vc . numtextcols/2+3 ; 
for ( i=l ; i<“6 ; i++ ) 

{ 

step“ycenter+( 192-1*6^ )*aspect_ratio; 
index*( i-3 )*textmax/5 ; 
spnntf  ( indexstring .  "X-  .  3g"  ,  index) ; 
textx*step*vc .numtextrows/vc .numypixels+1 ; 
if  ( i”3 ) 
continue; 
else 
{ 

_settextposition(textx, texty ) ; 

_outtext( indexstring) ; 

_moveto(xcenter , step) ; 

_lineto(xcenter+10 , step) ; 

} 

) 

max'=vc  .  nuraxpixels/4/niax; 
x=x/2; 

_setcolor(forcolor2) ; 

_setlogorg  (  xcenter,  ycenter  ); 
for(i”0 ; i<»n-l ; i++) 

{ 

position! i] [O]^nax*poles[i] [l)-x; 
positionti] ( 1 ) “(max*poles f i ) (2) -x)*aspect_ratio ; 
fortj”0; j<»i; j++) 
if C  j ! -i ) 

if  (positionUHOJ  —  position!  j !  [0]  &&  positionti) [1]  —  position!  j)  [  1] ) 
position! i ]  [ 0 ) “position! i ] [0]+3 ; 

if (positionti J (01<“0. 0) 

_put image (posit ion ! i ] [0! .position! i ] ! 1 ! , greenpole ,  _GXQR) ; 
else 

put image (positionti I [ 0 J , position !i ] [ 1 J , redpole,  _GXQR> ; 

} 

for  ( i“0 ;  i<*ea- 1 ;  i++  ) 

{ 

positionti] !0]*max*zeros!i] [l]-x; 
positionti] [l]*(max*zeros[i] (2 ] -x )*aspect_ratio ; 
for ( j”0 ; j<“i ; j++) 
if(j'.-i) 

if  (positionti  ]  1 0  ]  —  position!  j  I  [0]  &&  positionti  ][  II  —  position!.)]  [1] ) 
positionti] !0 Imposition! i] [0)+3 ; 

if (positionti] !0]<-0 . 0 ) 

_putimage(position( i ] [0] .positionli] [l]-l,greenzero,  _GX0R) ; 
else 

_put image (posit ion! i]10],position[i)[l]-l, redzero ,  _GX0R) ; 

} 


) 

/»****.**•*••*****»**•*  END  of  locus  output  program  •«****»***»***»•**•*«*»/ 
/•**..»**«»**•***»**«****»***  END  of  locus  file  **•*****••*••*******•**••*»/ 


APPENDIX  D 


BODE  PLOT  SOURCE  CODE 
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/*.**... 

/* 

/* 

/* 

/*«*•*** 

#mc  lude 

void 

double 


.»»»..***»*»  START  OF  BODE  ROUTINE  »»*»»•**»... 
This  program  computes  and  displays  Bode  plots. 


***********  START  GLOBAL  VARIABLES 
"variables ,h" 
bode_output( ) ; 
magmtudejnax , 

magnitude_min ,  /< 

phasemax , 
phasemin . 
points(6i0] [3] , 
points2 ( 640 ] (31 . 
pi-3. 14159265, 
phase_margin . 
gain_margm, 
gain_cross , 
phase_cross ; 


*/ 

/*  Include  variable  list  */ 
/»  Output  subroutine  */ 
/*  Maximum  magnitude  value  */ 
Minimum  magnitude  value  */ 

/*  Maximum  phase  value  */ 

/*  Minimum  phase  value  */ 

/*  Plotting  points  data  */ 

/*  Plotting  points  data  */ 


/*  Initial  lower  frequency  */ 
/*  Initial  upper  frequency  */ 

/*  Open- loop  gain  value  */ 

/*  Check  for  gain  crossover  */ 
/*  Check  for  phase  crossover  */ 

/..***(,**„*„„„***„„*„„*  END  GLOBAL  VARIABLES  ********•*».«*•**. ***«*„*„«**/ 

/.**»„.**.»***».*„„**  START  OF  MAIN  BODE  ROUTINE  ******•**»•****»**«•»****./ 

void  bodetint  n.int  m) 

{ 


float 


omegal , 
omegau , 
k_mitial ; 


gain_check*0 , 
phase_check=Q ; 


int  1. 
j  . 
rv, 

itterations , 
step , 
count” 1 ; 

float  k_final, 
deltax . 
test; 


/*  Loop  counter  */ 
/*  Loop  counter  */ 

/*  Return  value  */ 

/*  Number  of  frequency  itterations  */ 


/*  Final  gain  value  »/ 
/*  Horizontal  screen  step  size  */ 


double  magnitude”!) .  0  , 
phase«0 . 0 , 
omega , 

error«0 . 00001 , 
start“0 . 0 , 
stop“0 . 0 , 
deltaomega , 
temp; 

char  kbuff [10]  , 
lowerbuff [ 10] . 
upperbuff (10) : 

WINDOWPTR  windl, 
WINDOWFTR  wind2; 
WIFORM  frm; 
int  batrib; 
int  watrib; 
register  i; 


/*  Magnitude  running  total  */ 
/*  Phase  running  total  */ 

/*  Present  frequency  value  */ 

/*  Allowable  error  */ 

/*  Graph  starting  point  */ 

/*  Graph  stopping  point  */ 
/*  Frequency  step  size  */ 

/*  Temporary  value  */ 


/*  input  window  handle  */ 
/*  input  window  handle  V 
/*  input  form  handle  */ 
/*  border  strib  */ 
/*  window  atrib  */ 


/* 

*  Set  window  attributes; 

* 

*  border  -  blue/white  box 

*  window  -  white  beckground/black  letters 

* 

*/ 
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batrib  >  v_setatr (BLACK .WHITE, 0,0); 
watrib  *  v_setatr (BLACK, WHITE , 0 , 0 ) ; 

/* 

*  Open  window  at  0,0  -  78  cells  wide  and  23  cells  high 

V 


batnb  “  v_setatr  (WHITE,  BLACK,  0 , 0) ; 
watrib  =  v_setatr (WHITE , BLACK , 0 , 0 ) ; 

wind2  *  wn_open(Q , 5 , 4  . 70 , 11 , watrib , batrib ) ; 
if(!wind2)  exit(l): 

wn_titla(wind2 , "  Bode  Input  Menu  ".batrib); 

wn_print£(wind2 . "\n  Input  the  transfer  function  information  \n\n"); 

wn_printf (wind2, "  Example : \n" ) ; 

wn_pnntf  (wind2, "  Input  the  open-loop  transfer  function  gain,  K  :  50\n"); 

wn_printf  (wind2 , ''  Input  the  lower  value  of  the  frequency  plotting  range  :  0.01\n"); 

wn_printf (wind2 , "  Input  the  upper  value  of  the  frequency  plotting  range  :  1000\n"); 

wn_printf (wind2 , " - 

"); 

wn  prmtf  (wind2, "  Input  the  open- loop  transfer  function  gain,  K 
wnprintf (wind2 , ”  Input  the  lower  value  of  the  frequency  plotting  range 

wn_prrntf (wind2 , "  Input  the  upper  value  of  the  frequency  plotting  range 


frm*wn_f rmopn( 4 ) ; 
if(!frm)  exit ( 0 ) ; 

kbufffO]  *  0; 
lowerbuff[0]  «  0; 
upperbuff[0]  *  0: 

wn_g£loat(SET, frm, 0 ,wind2 ,8,57 , NSTR, watrib , ' 

’ ,&k_initial , 9 , 2 , 0 . 0 , 100000 .kbuff , NSTR, NSTR) ; 

wngf loat(SET, frm, 1 , wind2 , 9 , 57 , NSTR, watrib , ' 

'  , Aomegal ,11,6,0.0, 100000 , lowerbuf f , NSTR , NSTR ) ; 

wn_g£loat(SET, frm, 2 ,wind2, 10,57 , NSTR .watrib , ' 
' , &omegau ,11,3,0.0,1000000, upperbuf f , NSTR, NSTR) ; 


rvwn_£rmget(frm) ; 
wn_frmcls ( frm) ; 
wn_close(wind2) ; 


l f ( ! set_video_mode( ) )  /*  Set  video  mode  */ 

{ 

pnntf  ( "\nWarnlng!  This  program  doesn’t  support  this  machine's  video  card"); 
exit(O) ; 

> 

_getvideocon£ig(&vc ) ; 

setcolor(YELLCW) ; 

_settextcolor(YELLOW) ; 

_rectangle(_GBCRDER . 0 , 0 , vc . numxpixels- 1 , vc .numypixels- 1 ) ; 

_settextposition( 12,28); 

sprintft string, "Working!  Please  wait."); 

_outtext( string ) ; 
eettextpositiont 14 .28); 

eprintf(  string .  "nniMMHHlMHllHlIMHMilllHiHIHHiMliiH" J ; 

_outtext( string ) ; 
sprintf (string, " (" ) ; 

_eettextposition( 14 , 28) ; 

delta_x”vc ,numxpixels/vc .numtextcols ; 
itteretions“delta_x*(vc . numtextcols- 12 ) ; 
etep-itteretlons/25 ; 


- \n 

\n" ) ; 

\n" ) ; 

"); 
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test-log  10  (omegal)  ; 
start-floor ( test ) ; 
stop-ceil( loglO(omegau) ) ; 

omegal-pow( 10 . 0 , start > ;  /*  Set  omegal  next  lower  power  of  10  */ 

omegau-pow( 10 . 0 , stop ) ;  /*  Set  omegau  next  higher  power  of  10  */ 

delta_omega-(stop-start)/itterationa ;  /*  Set  frequency  step  size  V 


magnitude_max-0 . 0 ; 
magnitude_mm«0 -  0 ; 
phase _roax-0 . 0 ; 
phase_min-0 . 0 ; 


/* 

Initialize  max 

and  min  */ 

/* 

values  to  0 

*/ 

/* 

*/ 

/* 

*/ 

f or ( j-0 ; j<-i tterations ; j++ )  /*  Start  of  frequency  range  loop  V 

( 

k_final-k  initial; 

magnitude-0 . 0 ;  /*  Set  magnitude  and  phase  to  0  */ 

phase-0.0;  /*  for  start  of  each  frequency  */ 

omega-powt 10.0, start+del taomega* j ) ; 
l f ( j  >step* count ) 

{ 

outtext(string) ; 
count++ ; 

) 

for( 1-1 . l<^m ; 1++ ) 

{ 

i f ( f abs( poles [ 1 - 1 ] ( 1 ] )<error  &&  fabslpoles 1 1- 1 ) [2] )<error )  /*  Pole  at  the  Origin 

{ 

magni  tude^nagm  tude-20 . 0* log  10  ( omega  ) ; 
phase-phase -pi/ 2 ; 

) 


else  if ( fabs(poles[ 1- 1 ) 12] )<error )  /*  Real  Axis  Pole  */ 

{ 

k_final-k_final/( -poles 1 1- 1 ] [  1  ]  ) ; 
temp-1 . 0+Tomega/ ( -poles [ 1- 1 J 1 1 ] ) ) ; 
if ( temp! -0 . 0)  ; 

( 

magm tude-snagni  tude-20 . 0* log  10  (  f  aba ( temp) ) ; 
phase-phase -atan2( omega, (-poles[l-l)tl))); 

> 

) 

else  /*  Complex  Pole  */ 

{ 

magni  tude-snagm  tude-20 . 0*  log  10  ( sqrt(  poles  [  l- 1 )  [  1)  ‘poles  ( 1-1  ]  ( 1  ]  +  (omega -poles  [l-l][2])*(om 
ega-poles] 1-1)12)))); 

phase-phase+atan2< ( omega- poles [ 1 - 1 ) (21 ) .poles [ 1- 1 ) [ 1 ) ) ; 

) 


for ( 1-1 ; l<-n ; 1++ ) 

{ 

if ( fabs(zeros ( 1-1 ) ( U )<error  &&  f abs( zeros ( 1 - 1 ) (2) )<error )  /*  Zero  at  the  Origin 

( 

magni tude^nagni tude+20 . 0*log 10 (omega) ; 
phaae-phase+pi/2 ; 

> 

else  if ( fabs(zeros(l -1 J (2) )<error )  /*  Real  Axis  Zero  */ 

{ 

k_finel-k_ final* (-zeros! 1-1)  1 1) ) ; 
temp-1 . 0+ ( omega/ ( -zeros! 1-1)11))); 

1 f ( temp ! -0 . 0 ) ; 

{ 

magni tude^segni tude+20 . 0*logl0( tabs ( temp) ) ; 
phase«phase+atan2( omega , ( -zeros [ 1- 1) ! 1 ) ) ) ; 

> 

) 

else  /*  Complex  Zero  */ 
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( 


magm  tude^nagmtude+20 . 0* log  10 ( aqrt( zeros  [  1-1)11  )*zeros  [  1-1]  [  1 ]+(orosgs-zeros [ 1- 1 ]  (2)  )*(o«n 
ega-zeror [ 1- 1] (2) ) ) ) ; 

phase-phase  atan2( (omega- zeros l 1- 1 ) [2] ) , zeros [ 1-1] [ 1) ) ; 

) 

) 


l f (k_f inal  !-  1.0  &&  kfinal  >  0.0) 
magm  tude-magm tude+20 .0*loglQ(k_ final), ■ 
phase-phase* 180 . 0/pi ; 


points2 ( j )[ 1 ) -omega ;  /*  Store  plotting  data  */ 

points2(.) )  12) -magnitude;  /*  */ 

points2 l j ] [ 3 ] -phase ,  /*  */ 

while(fabs(phase)>180.0)  phase -phase- a ign( phase )* 360 ; 
if  (magm  tude>magni  tude_max )  magm  tudemax-magni  tude ; 
i f (magni tude<magni tudejnin )  magni tudemin-magni tude ; 
l f (phase>phas#_max )  phasejnax— phase ; 
if ( phase<phase_min)  phase_min-phase ; 


points  I j 1 11) -omega ; 
points | j ] [ 2 ) -magnitude ; 
points  t j ) [ 3 ] -phase ; 


/*  Store  plotting  data  */ 
/*  */ 
/*  */ 


for(j-0; j<-itterations; j++)  /*  Start  of  frequency  range  loop  */ 

( 

1 f ( j ! -0  ) 

< 

l  f  ( signCpoints2[  j  ]  [21  )  !-sign(pomts2[  j- 1]  12) ) )  /*  Check  for  phase  margin  V 

( 

phase_margin-180+points2 ( j - 1 1 1 3 ) -points2 [ j - 1 1 [ 2 1  * 

(points2( j )  [3)-pointa2f j-l)  (3)  )/(pomts2[  j]  [2]-polnta2[  j-1)  [2D; 

phase_cross-points2[ j- 1 ) t 1 1 -points2 [ j- 1) [21* 

7points2( j ]  [  1] -points2[ j- 1) [ 1 ) )/ (points2 [ j ) [ 2) -points2 [ j - 1 ) [2] ) ; 
phase_check“l ; 

) 

l  f  (  sign(points2[  j- 1 1  [3) +  180 )  !-  sign(pomts2 [  j  )  [3 1  +  180 ) )  /*  Check  for  gain  margin 

*/ 

( 

gain  raargm--points2[ j-l )  [2 1  -  ( point s2 [  j )  [2]  -pointa2[  j-l)  [2D* 

(-180-polnta2(j-ll[3) )/ (points2 [ j ] [ 3 ) -points2 [ j - 1 ) [ 3] ) ; 

gain_cross*( point s2 [ J - 1 1 [ 1 )+points2[ j ) [ 1 ) )/2 ; 
gain_check-l ; 

) 

} 

) 

outtext(strmg) ; 


bode_output( ) ;  /*  Call  output  routine  */ 

) 

QJ0  0F  Mg!))  BOOE  ROUTINE  *****»•**»**•*****••«•*«***/ 
STgRT  of  OUTPUT  SCREEN  •**•*••••**••*•*•***••***••»»/ 

void  bodeoutputt ) 

( 

lnt  1 , 

j. 

ltterationa; 

float  divisions, 
tlcstsrt, 
tlcstop. 
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rt_upper=0 , 
rt_lower-0 , 
delta_rt=0 , 
angle_upper“0 , 
angle_lower=0 , 
delta_angle-0 , 

x. 

y . 
yi. 
y2, 

tic_y , 
tic_x, 
step_size , 
textstepsize , 
delta_y , 
deltax, 
length; 

double  slope, 
x_slope , 
y_slope, 
x_hi gh , 
x_low. 
yhigh , 
ylow, 

char  string [20); 

aspect_ratio  *  (float)  (8.0  *  vc . numypixels )/ ( 6 . 0  *  vc .numxpixels) ; 
delta_x=vc .numxpixels/vc .numtextcols ; 
delta_y*vc  numypixels/vc . numtextrows ; 

_clearscreen(_GCLEARSCREEN) ; 

_setcolor(forcolorl) ; 
x_low=6*delta_x ; 

x_high-delta_x*(vc .numtextcols-6) ; 
y_high=delta_y ; 

y_low«delta_y*(vc .numtextrows-* )+delta_y/2+3*aspect_ratio; 
i tter at ions-deltax* (vc .numtextcols- 12) ; 

_r ec tangle (  GBQRDER, xlow, y_high, x_high, y_low ) ; 
y_high“delta_y*l . 5; 

/*  Gain  Heading  */ 
settextcolor(MAGENTA) ; 
setcolor (MAGENTA) ; 

_settextposition(0 , 0 ) ; 

sprintf (string, "  dB  Magnitude"); 

_outtext ( string ) ; 

_moveto(delta_x+l,delta_y/2-l) ; 

_lineto(3*delta_x+l ,delta_y/2-l ) ; 
rt_upper”ceil(magnitude_max) ; 
rt_lower”floor(magnitude_min) ; 
for( ; ;  ) 

{ 

if (fmod(rt_upper , 10 . 0) !“0 . 0)  rt_upper*rt_upper+l . 0 ; 

else  break ; 

} 

for( ;  ;  ) 

{ 

if (fmod(rt_lower , 10 . 0) !-0 . 0)  rt_lower-rt_lower-l . 0 ; 

else  break; 

} 

for( ; ; ) 

{ 

if (fmod(rt_upper-rt_lower , (vc .numtextrows- 5) ) !-0 .0) 

{ 

rt_upper-rt_upper+5 . 0 ; 
rt_lower~rt_lower-5 . 0 ; 

) 

else  break ; 

} 

delta_rt”(rt_upper-rt_lower)/(vc .numtextrows- 5) ; 
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tic_stop*6*delta_x; 

tic_start*tic_stop-delta_x/2; 

for (1*0 ; l<*(vc .nurotextrows-5) ; 1++) 

< 

tic_y=delta_y*(l . 5+1)- 1 ; 

_moveto(tic_start,tic_y) ; 

_lineto(tic  stop,tic_y); 

} 

for ( 1*0 ; l<=(vc .numtextrows-5) ; 1++) 

{ 

_  settextposition( 1+2 , 1 ) ; 

sprintf (string , "t5g" , ( rt_upper-delta_rt*l ) ) ; 

_outtext(string) ; 

} 

x_slope=(x_low-x_high)/ (loglO(omegal )-logl0 (omegau) ) ; 
y_slope=(-delta_y*(vc . numtextrows-5) )/(rt_upper-rt_lower ) ; 

_setlinestyle(0x0101) ; 

X*xl0W ; 

y=y_high-rt_upper*y_slope ; 

_moveto(x ,y ) ; 
x*x_high; 

_lineto(x,y ) ; 

_setlinestyle(OxFFFF) ; 

X*x_low-(logl0(omegal)-logl0(points[0] [1) ) )*x_slope; 
y*y_high-(rt_upper-points[0] [2) )*y_slope; 

_moveto(x,y ) ; 

for (1*1 ; l<=itterations; 1++) 

{ 

x*x_low- (loglO (omegal ) -loglQ (points [1] (1] ) )*x_slope; 
y=y_high-(rt_upper-points[l] (2] )*y_slope; 

_lineto(x,y ) ; 

} 

/*  Phase  Heading  */ 
settextcolor (GREEN) ; 

_setcolor (GREEN) ; 

length*sprintf ( string, "Phase  " ) ; 

_aettextposition(0,vc , numtextcols-length) ; 

outtext ( string ) ; 
y»delta_y/2-l; 

for( j*(vc .numtaxtcols-A )*delta_x; j<*(vc . numtextcols~2)*delte_x; j+“+2) 
_setpixel( j ,y ) ; 

tic_start*delta_x+(vc .numtextcols-6) ; 

tic_stop-tic_start+delta_x/2; 

for (1*0 ; <*(vc  .numtextrows-5) ;  1++) 

{ 

tic_y-delta_y*( 1 . 5+1 ) - 1 ; 

_moveto(tic_start , tic_y ) ; 

_lineto(tic_stop,tic_y ) ; 

} 

angle_upper*ceil(phase_max) ; 
angle_lower*f loor(phaaejnin) ; 
for( ; ;  ) 

{ 

if (fmod(angle_upper , 10 . 0) !*0 . 0)  angle_upper*angle_upper+l . 0 ; 
else  break; 

} 

for( ;  ; ) 

( 

if ( fmodtanglelower, 10 . 0) !*0 . 0)  angle_lower*angle_lower- 1 . 0 ; 
elaa  break; 

} 

for( ;  ;  ) 

{ 

if (fmod( (angleupper-anglalowar ) , (vc .numtextrows-5) ) !*0 .0) 

{ 

angle_upper*angla_upper+5 . 0 ; 
angle_lower*angle_lower-5 . 0 ; 
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} 

else  break; 

) 

delta_angle=(angle_upper-angle_lower )/(vc numtextrows-5) ; 

for ( 1=0 ; l<=(vc .numtextrows-5) ; 1++) 

{ 

_settextposition( 1+2 , vc . numtextcols-4 ) ; 

sprintf ( string , "l-g\370” , ( angle_upper-delta_angle*l ) ) ; 

_outtext  (  string ) ; 

} 

x_slope=(x_low-x_high)/(loglO(omegal)-loglO(omegau) ) ; 
y_slope=(-delta_y*(vc .numtextrows-5) )/(angle_upper-angle_lower ) ; 

x=x_low- ( log 10 (omegal ) -log  10 (points! 0] [ 1 ] ) )*x_slope ; 
y=y_high- (angle_upper-points[0] [3] )*y_slope ; 

_moveto(x,y ) ; 

for (1=2 ; l<*itterations ; l+*+2) 

{ 

x=x_low- (logl0( omegal) -log  10 (points! 1] [ 1] ) )*x_slope; 
y2=y_high-(angle_upper-points[l] !3) )*y_slope; 

_setpixel(x.y2) ; 

if (abs(y2-y )>2  &&  sign(  slope)=sign(y2-y ) ) 

{ 

yi=y ; 

fort  j=*l ;  j<=abs(y2-y  )/2 ;  j+=+2) 

{ 

yl=yl+sign(y2-y )*2 ; 
setpixel(x-l ,yl) ; 

T 

for(j=l; j<=abs(y2-y)/2; j+=+2) 

{ 

yl=yl+sign(y2-y )*2; 

_setpixel(x.yl) ; 

} 

} 

slope=y2~y ; 
y=y2 : 

) 

_settextcolor( forcolorl ) ; 

_setcolor(forcolorl ) ; 

di visions=log 10 ( omegau ) -logl0( omegal ) ; 
step_size=(vc .numxpixels-12+delta_x)/divisions ; 
text_step_size=(vc . numtextco Is- 12) /divisions ; 
tic_start=delta_y*(vc  .numtextrows-'i  )+delta_y/2+3*aspect_ratio ; 
tic_stop=tic_start+delta_x/2 ; 
for(l*0 ; l<=divisiona ; 1++) 

{ 

tic_x=6*delta_x+step_size*l ; 
moveto(tic_x,tic_start) ; 

_lineto(tic_x,tic_stop) ; 
if (l!=divisions) 

{ 

for( j”l;j<«8; j++) 

{ 

x=x_low- (logl0( omegal)- log  10 (omegal*pow( 10, 1)*( j+1) ) )*x_slope 
_moveto(x, tic_start) ; 

_lineto(x,tic_stop-delta_x/4) ; 

) 

} 

length-sprintf (string , "Ig" ,omegal*pow(10 , 1) ) ; 

^settextposition(23 , 7+text_step_size*l- length/2) ; 
outtextt string) ; 

T 

string!0]='\0' ; 
atmcpy(string,nstring,25) ; 
length-strlen(string) ; 
for (1=1 ; l<=l«ngth+l ; 1++) 
natrlng(l-l]-atring[l] ; 
atring(0)*'\0' ; 
if (k  initialI-1) 
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{ 

sprintf ( string , "tg" ,k_initial ) ; 
strncat( string .nstring , 25) ; 
strncpy(nstring, string, 25) ; 

} 

if (k_initial-=l  &&  length<2) 

{ 

sprintf (string , "tg" ,k_initial) ; 
strncattstring, nstring, 25) ; 
strncpy (nstring, string  25) ; 

} 

length=strlen(nstring) ; 
i f ( length>maxchar ) 

{ 

skipd=skipd+(length-maxchar )/2; 
maxchar=length ; 

} 

else 

skipn=(maxchar-length)/2+2 ; 

_settextposition( (vc . nuratextrows- 1 ) . 78-maxchar+skipn) ; 

_outtext(nstring) ; 

_settextposition( (vc . numtextrows ) , 78-maxchar+skipd ) ; 

_outtext(dstring) ; 
y=460*aspect_ratio; 

_moveto(630-maxchar*(vc .numxpixels/vc .numtextcols) ,y ) ; 

_lineto(630 ,y ) ; 

min=vc ,numypixels-43*aspect_ratio; 

_rec tangle (  GBORDER, 0 ,min+2 , (vc.numxpixels-1) , (vc .nuraypixels-1) ) ; 

movetot (vc .numxpixels-2- (maxchar+2)*vc . numxpixels/vc .numtextcols) ,min+2) ; 

_lineto( (vc .numxpixels-2- (maxchar+2)*vc .numxpixels/vc .numtextcols) , (vc . nuraypixels- 1 ) ) 

length-sprintf (string,"  Bode  Plot"); 

_settextposition(0, (vc .numtextcols- length) /2) ; 
outtext(string) ; 

x-(vc  numxpixels-( length* l)*delta_x)/2 ; 

_moveto(x, 0) ; 

_lineto(x, delta_y ) ; 

x*(vc .numxpixels* (length* l)*delta_x) /2-del ta_x-2; 

moveto ( x , 0 ) ; 

_lineto(x,delta_y ) ; 

if (phase_check--l) 

{ 

length-sprintf (string , "Phase  Margin  -  t.lg\370  ",  phasemargin) ; 

_settextposition(24 ,2) ; 

_outtext( string) ; 

length-sprintf ( string , Freq.  *  t.lg".  phase_cross ) ; 

_settextposition(24 , 24 ) ; 

_outtext( string) ; 

} 

else 

{ 

length-sprintf (string , "Phase  margin  is  not  defined."); 

_settextposition(24,2) ; 

_outtext( string) ; 

) 

if ( gain_check--l ) 

{ 

length-sprintf (string , "Gain  Margin  “  t.lg  dB",  gain_margin) ; 

_settextposition(25 , 2) ; 

_outtext( string) ; 

length-sprintf (string , "8  Freq.  -  t.lg",  gaincross); 

_settextpoaition(25 , 24 ) ; 

_outtext(string) ; 

} 

else 

{ 

length-sprintf (string, "Gain  margin  is  not  defined."); 

_settextpoaition(25,2) ; 
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_outtext ( string ) ; 

> 

_rectangle(_GBQRDER, 0 , 0 ,vc . numxpixels- 1 , vc .numypixels-1 ) ; 

> 

/************************  END  OF  OUTPUT  SCREEN  **********************★*****/ 


/**************************  fjij])  OF  BODE  ROUTINE  ******************★********/ 


APPENDIX  E 


MODIFIED  WINDOWS  SOURCE  CODE 
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/*****************»  START  OF  MODIFIED  WINDOW  ROUTINES  *«****************»**/ 


/*  */ 

/*  These  two  routines  are  modified  versions  of  the  popup  and  formget  V 
/*  routines  supplied  in  the  windowing  package  */ 
/*  */ 
/*  The  Window  BOSS  */ 
/*  by  Phil  Mongelluzzo  */ 
/*  '  */ 
/*  Revision:  08.15.88  */ 
/*  */ 
/*  Star  Guidance  Consulting,  Inc.  */ 
/*  273  Windy  Drive  */ 
/*  Waterbury,  Connecticut  06705  */ 
/*  */ 
/*  (203)  571-2449  */ 
/*  */ 


/A***************************************************************************/ 

/* 

**  POPUP  -  Popup  Menu/Window  Driver 
** 

**  Copyright  (c)  1985-1988  Philip  A.  Mongelluzzo 
**  All  rights  reserved. 

*/ 


^include  "windows. h" 

tthindef  UARROW 
#undef  DARROW 
#undef  LARROW 
#undef  RARROW 
#undef  BS 
#undef  DEL 
#undef  RET 
#undef  ESC 


♦define 

UARROW 

0x48 

♦define 

DARROW 

0x50 

♦define 

LARROW 

0x4b 

♦define 

RARROW 

0x4d 

♦define 

BS 

OxOe 

♦define 

DEL 

0x53 

♦define 

RET 

Oxlc 

♦define 

ESC 

0x01 

♦define 

SPACE 

0x39 

♦define 

ZERO 

0x00 

♦define 

ONE 

0x02 

♦define 

TWO 

0x03 

♦define 

THREE 

0x04 

♦define 

FOUR 

0x05 

♦define 

FIVE 

0x06 

♦define 

SIX 

0x07 

♦define 

SEVEN 

0x08 

♦define 

EIGHT 

0x09 

♦define  v_setrev( atrib )  ( ( atrib&0*88 )  ] 

/* 

*»  Popup  structure  templates 

*/ 


/*  window  header  */ 

/*  dont  argue  with  windows. h  */ 


/*  define  key  codes  */ 


(<atrib»4)&0x07)  |  ( ( atrib«4  )&0x70 )  ) 


struct  m 
int  r; 
int  c ; 
char  *t; 
int  rv; 

); 


/*  menu  item  template  */ 
/*  row  */ 

/*  col  */ 

/*  text  */ 

/*  return  value  */ 
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struct  pmenu  { 

WINDCWPTR  wpsave; 
int  winopn; 
int  lndx; 
int  fm; 
int  lm; 

struct  mitem  scrn[25] ; 
k 


/*  popup  menu  structure  *7 
/*  place  to  hold  window  id  */ 
/*  leave  window  open  flag  */ 
/*  last  index  */ 

/*  first  menu  item  index  */ 

/*  last  menu  item  index  */ 

/*  a  bunch  of  menu  items  *7 


/* 

★  ****■**** 

*  popup  * 

********* 

V 

popup (page , row, col .width , height, atrib 
int  page; 
int  row,  col; 
int  width,  height; 
struct  pmenu  *mx; 
int  cflag; 
int  atrib,  batrib; 


batrib ,mx , cflag ) 

/*  video  page  */ 

/*  window  -  upper  row/col  */ 

/*  window  -  height  &  width  */ 

/*  pointer  to  popup  menu  struct  */ 
/*  close  on  return  strike  flag  */ 
/*  attributes  -  window  &  border  */ 


int  i ; 

WINDOWPTR  w; 
int  c ; 


/*  scratch  integer  */ 

/*  window  pointer  */ 

/*  key  scan  code,,char  */ 


if (mx->winopn)  { 
goto  dO; 

) 

mx->lndx  =  -1; 

w  «  wn_open(page ,  row,  col.  width, 

wn_sync (w , FALSE ) ; 

mx->wpsave  “  w; 

if  (w  ”=  NULL)  return(99); 

mx->winopn  =  TRUE; 

i  =  0; 


/*  window  is  still  open  */ 

/*  enter  processing  loop  */ 

/*  set  index  out  of  range  */ 
height,  atrib,  batrib); 

/*  sync  cursor  */ 

/*  save  pointer  */ 

/*  out  of  mem  or  just  fubar  */ 

/*  say  window  is  open  */ 

/*  init  index  */ 

for  as  long  as  we  have  data  *7 
c,  mx->scrn[i] .t,  atrib); 


while(mx->scrn(i) ,r  !=  99)  {  /* 

wn_putsa(w,  mx->scrn [ i ] , r ,  mx->scrnti) 
i++; 

) 


w  ”  mx->wpsave; 
i  *  mx->lndx; 

if(i  <  mx->fm)  i  *  mx->fm; 
if(i  >  mx->lm)  i  ■  mx->fm; 
while (TRUE)  { 

wn_putsa(w,  mx->scrn[i J . r , 
c  *  v_getch( )  >>8; 


/*  restore  pointer  */ 

/*  BLINDLY  assume  that  we're  back  */ 
/*  and  reset  if  "i"  is  now  out  of  */ 
/*  range  -  (dumb,  but  it  works)  */ 

/*  till  we  exit  */ 

mx->scrn[i)  .c,  mx->scrn[ i ) . t ,  v_setrev( atrib ) ) 


if ( c  --  ESC)  { 
if (cflag)  { 
wn_close(w) ; 
mx->winopn  -  FALSE; 

} 

mx->lndx  “  i; 
return(i ) ; 

) 

if(c  “  RET)  { 
if (cflag )  { 
wn_close(w) ; 
mx->winopn  ■  FALSE; 

) 

mx->lndx  “  1 ; 
return(mx->scm[  i )  .  rv) ; 


/*  ESC  (user  wants  out)  -  swk  */ 

/*  close  window  on  return  strike  ?  */ 
/*  close  the  window  */ 

/*  say  window  is  closed  */ 

/*  remember  last  indx  V 
/*  return  with  rv  •/ 

/•  swk  RETURN  */ 

/*  close  window  on  return  strike  ?  */ 
/*  close  the  window  */ 

/*  say  window  is  closed  */ 

/*  remember  last  indx  */ 

/*  return  with  rv  */ 
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) 


if (c  =  ONE)  { 
if(cflag)  { 
wn_close(w) ; 
mx->winopn  =  FALSE; 
> 


/*  swk  RETURN  */ 

/*  close  window  on  return  strike  ? 
/*  close  the  window  */ 

/*  say  window  is  closed  */ 


mx->lndx  =  i ; 
return(l) ; 

) 

if(c  ==  TWO)  { 
if ( cf lag )  { 
wn_close(w) ; 
mx->winopn  =  FALSE; 

} 

mx->lndx  =  1 ; 
return(2) ; 

} 

if (c  —  THREE)  { 
if(cflag)  { 
wn_close(w) ; 
mx->winopn  =  FALSE; 
} 

mx->lndx  =  i; 
return ( 3 ) ; 

) 

if (c  =  FOUR)  { 
if (cflag)  { 
wn_close(w) ; 
mx->winopn  *  FALSE; 
> 

mx->lndx  “  i; 
return(4) ; 


/*  remember  last  indx  */ 

/*  return  with  rv  */ 

/*  swk  RETURN  */ 

/*  close  window  on  return  strike 
/*  close  the  window  */ 

/*  say  window  is  closed  */ 

/*  remember  last  indx  */ 

/*  return  with  rv  */ 

/*  swk  RETURN  */ 

/*  close  window  on  return  strike 
/*  close  the  window  */ 

/*  say  window  is  closed  */ 

/*  remember  last  indx  */ 

/*  return  with  rv  */ 

/*  swk  RETURN  */ 

/*  close  window  on  return  strike 
/*  close  the  window  */ 

/*  say  window  is  closed  */ 

/*  remember  last  indx  */ 

/*  return  with  rv  V 


1 f ( c  —  DARROW)  , 
if (c  -«  RARRCW) 
if (c  —  LARROW)  i 
if (c  =  U ARROW)  , 
if (c  —  SPACE  ij 

wn_putsaCw,  mx 

if (c  =  SPACE) 
i++; 


=  SPACE;  /* 
-  SPACE;  /» 
=  BS;  /* 
=  BS;  /• 


c  -=  DEL  ] |  c  BS)  ■ 
>scrn[i].r,  mx->scrn[ 

/* 

/* 


swk  conversion  */ 
swk  conversion  */ 
swk  conversion  */ 
swk  conversion  V 

. ) . c ,  mx->scrn[ i ) . t , 
foreward  ??  */ 
bump  index  */ 


atrib) ; 


else 
i-- ; 

if(i  <  mx->fm)  i  “  mx->lm; 
if(i  >  mx->lm)  i  “  mx->fm; 

} 

) 

wn_close(w) ; 
mx->winopn  -  FALSE; 
return ( 99 ) ; 


/*  decrement  */ 

/*  wrap  on  either  */ 

I*  end  */ 

/*  end if  */ 

/*  end  while  */ 

/*  bye  bye  */ 

/*  say  window  is  closed  */ 
/*  nothing  selected  */ 


? 


*/ 


*/ 


*/ 


*/ 


#ifdef  COMMENTS 

/* 

********** 

*  qpcpup  * 
********** 

*/ 

#endif 


/*  quick  popup. . .  */ 

/*  this  is  a  hack  of  popup  */ 
/*  but  the  code  fragment  */ 

/*  could  be  of  value  */ 


WINDCWPTR  qpopuptpage , row, col .width, 

int  page; 

int  row,  col; 

int  width,  height; 

int  atrib,  batrib; 

struct  pmanu  *mx; 

{ 

int  i ; 

WINDCWPTR  w; 

w  "  wn_open(page,  row,  col,  width, 
i  -  0; 


height , atrib .batrib ,mx) 

/*  video  page  V 

/*  window  -  upper  row/col  */ 

/*  window  -  height  &  width  */ 

/*  atributes  -  window  &  border  V 
/*  pointer  to  text  struct  */ 

/*  scratch  integer  */ 

/*  window  pointer  */ 

height,  atrib,  batrib); 

/*  init  index  */ 
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while(mx->scrn[i) ,r  !*  99)  {  /*  for  as  long  as  we  have  data  */ 

wn_puts(w,  mx->scrn[i]  .r,  mx->scrn[  i  ]  .  c ,  nuc->scrn[i  ]  .  t) ; 
i++; 

} 

return(w) ; 


/*  End  V 


/* 

**  The  Window  BOSS'S  Data  Clerk 

**  Copyright  (c)  1988  -  Philip  A.  Mongelluzzo 

*"  All  rights  reserved. 

*  A 

**  wn_frmget  -  Get  (read)  data  entry  form 
** 

**  Copyright  (c)  1988  -  Philip  A.  Mongelluzzo 
**  All  rights  reserved. 

** 

V 

/*  standard  stuff  */ 

/* 

*  ★•Hr********** 

*  wn_f rmget  * 

************* 

*/ 


/* 

**  wn_frmget(frm) 

** 

**  int  nfields  -  number  of  fields  in  form. 


**  RETURNS: 

** 

**  TRUE  -  indicating  all  fields  of  the  form  in  question  have  been 
**  fetched  and  verified  (where  required). 

** 


** 

** 


** 


or 

Never  Returns! ! 


**  NOTES: 

aa 


** 

aa 

aa 

** 

aa 

aa 

aa 

aa 

aa 

aa 


This  code  can  be  used  as  is  or  can  be  customized  to  suite  each 
applications  need.  The  section  to  customize  is  at  the  tail  end 
of  this  file.  As  distributed,  logic  diplays  all  prompts  and 
display  fields,  positions  to  the  first  field,  performs  data 
entry  on  a  field  by  field  basis  from  the  first  to  the  last 
(allowing  editing  along  the  way),  asks  for  a  confirmation  to 
accept  the  fields  on  the  form  after  the  last  field  is  entered, 
either  accepts  the  form  or  drops  into  edit  mode  for  all  the 
fields  on  the  form  starting  at  the  first  field. 


wn_f rmget  will  not  return  until  all  data  has  been  entered  and 
verified  (where  required) . 


aa 

aa 


This  routine  must  be  called  after  wn_frmopn,  and  before  wn_frmcls. 


V 


/* 


*  wn_f rmget  * 


V 


♦define  FLD  frm[indx] 


/*  some  shorthand  */ 


wn_frmget2(  frrn) 


/*  input  form  processor  V 
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WIFORM  frm; 

{ 

int  mdx; 
int  rv; 
int  c : 

WINDOWPTR  wn; 

static  char  *msg  =  "[Press  ENTER  to 


/*  array  of  field  pointers  */ 

/*  input  field  index  */ 

/*  wn_g???'s  return  value  */ 
/*  scratch  */ 

/*  ray  CWN  window  */ 

Accept,  any  other  key  to  Edit]"; 


indx  =0;  /*  set  index  */ 

while(TRUE)  {  /*  display  output  fields  */ 

if (FLD->pself  !  =  (char  *)FLD)  /*  corrupted  memory  check  */ 

wns_ierr ( "wm_frmget@wn_puts" ) ;  /*  die  if  memory  is  mangled  */ 

if( !FLD->fcode)  break;  /*  all  done  with  output  fields  V 

wn_puts(FLD->wn,  FLD->row,  FLD->col,  FLD->prmpt); 
indx++;  /*  bump  index  */ 


begin ; 

indx  =0;  /*  reset  index  */ 

while(TRUE)  {  /*  fetch  input  fields  */ 

if ( ! FLD->f code )  break;  /*  all  done  */ 

switch  (FLD->fcode)  {  /*  based  on  function  code  */ 

case  GDATE: 

if (FLD->pself  !=  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr( "wra_frmget@GDATE" ) ; 
rv  ■  wn_gdate(XEQ,NFRM,NFLD,FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vip, 
FLD->v3.vip,  FLD->v4.vcp, 
FLD->v5.vcp,FLD->v6.vcp) ; 

break  ; 
case  GTIME: 

if (FLD->pself  !=  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr("wm_frmget@GTIME") ; 
rv  «  wn_gtime(XEQ,NFRM,NFLD,FLD->wn,  FL0->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vip, 
FLD->v3. vip.  FID->v4 . vcp, 
FLD->v5.vcp,FLD->v6.vcp) ; 

break  ; 
case  GINT: 

if (FLD->pself  !“  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr ( "wm_f rmget@GINT" ) ; 
rv  “  wn_gint(XEQ.NFRM,NFLD , FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vi, 
FLD->v3.vi,  FLD->v4.vi,  FLD->vS.vcp, 
FLD->v6,vcp,FLD->v7.vcp) ; 

break ; 
case  GUINT: 

if (FLD->pseif  !  =  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr ( "wm_frmgetSGUINT" ) ; 
rv  -  wn_guint(XEQ , NFRM ,NFLD , FLD->wn ,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD-?fill,  FLD->vl . vuip ,  FLD->v2.vi, 
FLD->v3.vui,  FLD->v4.vui,  FLD->v5.vcp, 

FLD->v6 . vcp , FLD->v7 . vcp) ; 

break : 
case  GLONG: 

if (FLD->pself  !-  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr ("wro_frmget6GLONG" ) ; 
rv  “  wn_g long (XEQ, NFRM, NFLD , FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib ,  FLD->fill,  FLD->vl .  vlp ,  FLD->v2.vi, 
FLD->v3.vl,  FLD-xv4.vl,  FLD->v5.vcp, 

FLD->v6 , vcp, FL0->v7 .vcp) ; 

break ; 

case  GFLOAT ; 

if (FLD->pself  !-  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr("wm_frmget8GFLOAT") ; 
rv  -  wn_gfloat(XEQ, NFRM, NFLD, FLD->*m,  FLD->row, 
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} 


FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD~>£ill,  FLD->vl.vfp,  FLD->v2.vi, 
FLD->v3.vi,  FLD->v4.vf,  FLD->v5.v£,  FLD->v6.vcp, 


FLD->v7 . vcp, FLD->v8 . vcp) ; 

break ; 

case  GFHONE: 

if (FLD->pself  !=  (char  *)FLD)  /*  corrupted  memory  check  */ 
vms_ierr("vm_frmget@GPHONE'' ) ; 
rv  «  vm_gphone(XEQ , NFRM,NFLD , FLD->wn ,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->£ill,  FLD->vl.vip,  FLD->v2.vip, 
FLD->v3.vip,  FLD->v4.vcp, 

FLD->v5.vcp,FLD->v6.vcp) ; 

break  ; 
case  GTEXT : 

if (FLD->pself  !-  (char  *)FLD)  /*  corrupted  memory  check  */ 
vms_ierr( "wm_f rmgetSGTEXT" ) ; 
rv  -  wn_gtext (XEQ , NFRM, NFLD , FLD->wn ,  FLD->row, 

FLD->col+strlen(FLD->pnnpt ) .  NSTR, 

FLD->atrib,  FLD->£ill,  FLD->vl.vi,  FLD->v2.vcp, 

FLD->v3 , vcp, FLD->v4 .vcp) ; 

break  ; 
case  GBOOL: 

if (FLD->pself  !=  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr(''wm_frmget@GBOOL") ; 
rv  =  wn_gbool(XEQ, NFRM, NFLD, FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt ) ,  NSTR, 

FLD->atrib,  FLD->£ill,  FLD->vl.vip,  FLD->v2.vcp, 
FLD~>v3 . vcp,FLD->v4 . vcp) ; 

break ; 

case  DTEXT :  /*  display  text  */ 

if (FLD->psel£  !  =  (char  * )FLD )  /*  corrupted  memory  check  */ 

wns_ierr("wm  f rmgetSDTEXT" ) ; 

rv  ■  wn_dtext(XEQ, NFRM, NFLD , FLD->wn ,  FLD->row,  FLD->col,  FLD->prmpt); 
break ; 


} 

if (rv  !«  BKTAB  &&  rv  !-  UARROW) 
indx++; 

else  iftrv  “  ESC) 
return(O) ; 
else  { 
indx- - ; 

if (indx  <-  0)  indx  •  0; 


/*  end  switch(frm.code)  */ 
/*  previous  field  ??  */ 

/*  no.  just  bump  index  */ 


/*  previous  field  request  */ 
/*  decrement  index  */ 

/*  but  dont  be  stupid!  */ 


> 


/ftftft* ft* ****************************** ************************************* 
ft  * 

*  CUSTOMIZE  THE  SECTION  BELOW  ONLY  * 


**************************** 


) 

/*  End  */ 
andif ; 


/* 

**  The  Window  BOSS'S  Data  Clerk 

**  Copyright  (c)  1986  -  Philip  A.  Mongelluzzo 

**  All  rights  reserved. 

** 

**  wn_frmget  -  Get  (read)  data  entry  form 
** 

**  Copyright  (c)  1988  -  Philip  A.  Mongelluzzo 
**  All  rights  reserved. 

** 

*/ 

/* 
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************* 

*  wn_f  rmget  * 
************* 


**  wn_frmget ( frm) 


nfields  -  number  of  fields  in  form. 


TRUE  -  indicating  all  fields  of  the  form  in  question  have  been 
fetched  and  verified  (where  required). 


**  Never  Returns!! 

** 

**  NOTES: 

** 

**  This  code  can  be  used  as  is  or  can  be  customized  to  suite  each 
**  applications  need.  The  section  to  customize  is  at  the  tail  end 

**  of  this  file.  As  distributed,  logic  diplays  all  prompts  and 

**  display  fields,  positions  to  the  first  field,  performs  data 
**  entry  on  a  field  by  field  basis  from  the  first  to  the  last 
**  (allowing  editing  along  the  way),  asks  for  a  confirmation  to 

**  accept  the  fields  on  the  form  after  the  last  field  is  entered, 

**  either  accepts  the  form  or  drops  into  edit  mode  for  all  the 

**  fields  on  the  form  starting  at  the  first  field. 

*  * 

**  wn_frmget  will  not  return  until  all  data  has  been  entered  and 
**  verified  (where  required). 

** 

**  This  routine  must  be  called  after  wn_frmopn,  and  before  wnfrmcls. 


************* 

*  wnfrmget  * 
************* 


♦define  FLD  frm[indx] 


/*  some  shorthand  */ 


wn_frmget(frm) 
WIFORM  frm; 


/*  input  form  processor  V 
/*  array  of  field  pointers  */ 


int  mdx; 
int  rv; 
int  c ; 

WINDOWPTR  wn; 


/*  input  field  index  V 
/*  wn_g???'s  return  value  */ 
/*  scratch  */ 

/*  my  CWN  window  */ 


static  char  *msg  ■  "[Press  ENTER  to  Accept,  any  other  key  to  Edit)' 


indx  “  0; 
while (TRUE)  { 

if (FLD->pself  !-  (char  *)FLD) 
wns_ierr ( "wm_frmget8wn_puts" ) ; 
if ( !FLD->fcode)  break; 
wn_puts(FLD->wn,  FLD->row,  FLD->col 
indx++ ; 


/*  set  index  V 

/*  display  output  fields  */ 

/*  corrupted  memory  check  */ 

/*  die  if  memory  is  mangled  */ 

/*  all  done  with  output  fields  */ 
FLD->pnnpt); 

/*  bump  index  */ 


begin: 
indx  ■  0; 
while(TRUE)  { 

if ( !FLD->fcode)  break; 
switch  (FlD->feode)  { 
case  GDATE: 

if(FLD->pself  !-  (char  *)FLD> 


/*  reset  index  */ 

/*  fetch  input  fields  */ 

/*  all  done  */ 

/*  based  on  function  code  */ 


/•  corrupted  memory  check 
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wns_ierr("wm_frmget8GDATE" ) ; 
rv  -  »m_gdate(XEQ, NFRM, NFLD, FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vip, 
FLD->v3.vip,  FLD->v4.vcp, 
FLD->v5.vcp,FLD->v6.vcp) ; 

break ; 
case  GTIME: 

if (FLD->pself  !*  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr( "wm_f rmgetSGTIME" ) ; 
rv  =  wn_gtime(XEQ, NFRM, NFLD, FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vip, 
FLD~>v3.vip,  FLD->v4.vcp, 

FLD->v5 . vcp , FLD->v6 . vcp) ; 

break  ; 
case  GINT: 

if (FLD->pself  !=  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr ( "wm_£rmget8GINT" ) ; 
rv  -  wn_gmt (XEQ , NFRM , NFLD  ,FLD->wn ,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vi, 
FLD->v3.vi,  FLD~>v4.vi,  FLD->v5.vcp, 

FLD->v6 . vcp, FLD->v7 . vcp) ; 

break ; 
case  GUINT: 

if (FLD->pself  !  =  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr("wm_frmget8GUINT" ) ; 
rv  =  wn_guint(XEQ, NFRM, NFLD, FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl . vuip ,  FLD->v2.vi, 
FLD->v3.vui,  FLD->v4.vui,  FLD->v5.vcp, 

FLD->v6 . vcp, FLD->v7 . vcp) ; 

break ; 
case  GLONG: 

if (FLD->pself  !«  (char  *)F1D)  /*  corrupted  memory  check  */ 

wns_ierr("wm_frmget8GL0NG" ) ; 
rv  ■  wnglongt XEQ, NFRM, NFLD, FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD~>atrib ,  FLD->fill ,  FLD->vl . vlp,  FLD->v2.vi, 
FLD->v3.vl,  FLD->v4,vl.  FLD->v5.vcp, 
FLD->v6.vcp,FLD->v7,vcp) ; 

break ; 

case  GFLOAT: 

if (FLD->pself  !“  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr("wm_£nnget@GFLOAT" ) ; 
rv  -  wn_gfloat(XEQ, NFRM, NFLD, FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->prmpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vfp,  FLD->v2,vi, 
FLD->v3,vi,  FLD->v4.vf,  FLD->v5.vf,  FLD->v6.vcp, 
FLD->v7 . vcp, FLD->v8 . vcp) ; 

break ; 

case  GFRONE: 

if (FLD->pself  !■  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr ( "wm_f rmgetBGPHONE" ) ; 
i  -  wn_gphone(XEQ, NFRM, NFLD , FLD->wn,  FLD->row, 

FLD->col+strlen(FLD->pnnpt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vip, 
FLD->v3.vip,  FLD->v4,vcp, 

FLD->v5.vcp,FLD->vE  vcp); 

break ; 
case  GTEXT : 

if (FLD->pself  !■  (char  *)FLD)  /*  corrupted  memory  check  */ 
*ms_ierr ( "wm_£rmget#GTEXT" ) ; 
rv  -  wn_gtext( XEQ, NFRM, NFLD, FLD->wn,  FLD->row, 

FLD->eol+atrlen(FLD->prnipt> ,  NSTR, 

FLD->atrib ,  FLD->f ill ,  FLD->vl,vi,  FLD->v2.vcp, 
FLD->v3 . vcp, f LD->v4 , vcp) ; 

break; 
caae  GBOOL: 


if(FLD->paelf  !-  (char  *}FLD)  /*  corrupted  memory  check  */ 
wna_ierr(  "**m_£rmget#GBOOL'‘ ) ; 
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rv  «  wn_gboo 1 ( XEQ , NFRM , NFLD , FLD- >wn ,  FLD->row, 

FLD->col+strlen(FLD->prnipt) ,  NSTR, 

FLD->atrib,  FLD->fill,  FLD->vl.vip,  FLD->v2.vcp, 
FLD->v3 . vcp, FLD->v4 . vcp) ; 

break : 

case  DTEXT:  /*  display  text  */ 

if (FLD->pself  !=  (char  *)FLD)  /*  corrupted  memory  check  */ 
wns_ierr ( "wm_f rmgetSDTEXT" ) ; 

rv  =  wn_dtext(XEO , NFRM , NFLD , FLD->wn ,  FLD->row,  FLD->col,  FLD->prmpt) 
break  ; 


> 

/* 

end  switch(frm.code)  */ 

if (rv  ! =  BKTAB  &&  rv  !=  UARROW) 

/* 

previous  field  ??  * / 

indx++ ; 

/* 

no,  just  bump  index  */ 

else  if(rv  =*  ESC) 

return( 0) ; 

else  { 

/* 

previous  field  request  */ 

indx-- ; 

/* 

decrement  index  */ 

if(indx  <=  0)  indx  =  0; 

/* 

but  dont  be  stupid!  */ 

} 

} 


/★A*********************************  **************************  ************ 
*  * 


* 

* 


CUSTOMIZE  THE  SECTION  BELOW  ONLY 


********* A****************************************************************/ 


wn,!,wn_open  ( 1000 ,24,0 
if(!wn)  exit(l); 
wnputs  (wn,  0 , 0  ,msg)  ; 
c  *  v_getch( )  &  Oxf f 
wn_close(wn) ; 
if (c  “  CR) 

re turn (TRUE) ; 
else 

goto  begin; 


/*  PRESS  Return  Stuff  */ 
(int)strlen(msg) , 1. (RVIDEO) .NVIDEO) ; 

/*  display  message  */ 

/*  fetch  response  */ 

/*  make  message  go  away  */ 
/*  ENTER  ??  */ 

/*  return  if  so  */ 

/*  else  edit  the  */ 

/*  form  !  */ 


/*  End  */ 


APPENDIX  F 


FACTORED  POLYNOMIAL  SOURCE  CODE 
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/***.*******.***********  FACTORED  polynomial  PROGRAM  ************************/ 


/*  */ 

/*  This  program  prompts  the  user  for  the  transfer  function  information.  */ 

/*  This  transfer  function  is  in  the  factored  polynomial  form.  */ 

/*  */ 


/a***********************  START  GLOBAL  VARIABLES  **************************/ 
^include  "variables . h"  /*  Include  variable  list  */ 

/**************************  FND  GLOBAL  VARIABLES  a**************************/ 
/**.*************»****,»*»  START  OF  FACTORED  ROUTINE  **********»*****»»*****/ 
int  factored() 

{ 


int 

j , 

/*  Loop  counter  */ 

k , 

/*  Loop  counter  */ 

1, 

/*  Loop  counter  V 

rv; 

/*  Return  value  */ 

float 

coef [9] [3] , 

/*  Temporary  coefficient  values  */ 

a.b  ,  c  ,  d; 

/*  Temporary  quadratic  values  */ 

char 

input : 

/*  User  input  variable  V 

WINDOWPTR  w2 ; 

/*  Factored  input  window  handle 

V 

WINDOWPTR  w3 : 

/*  Warning  window  handle 

V 

WINDOWPTR  w 4; 

/*  Final  equation  window  handle 

V 

WINDOWPTR  m; 

/*  Lower  prompt  window 

*/ 

WIFORM 

frm; 

/*  Factored  input  menu  form  handle 

*/ 

WIFORM 

f  rm2 ; 

/*  Coefficient  input  menu  form  handle 

V 

int 

batrib ; 

/*  border  atrib 

V 

int 

watrib ; 

/*  window  atrib 

*/ 

register  i; 

char 

nlbuff [ 10) , 

/****,*********************  / 

nqbuff [ 10] , 

/*  */ 

dlbufftlO) , 

/*  Input  buffer  arrays  */ 

dqbuff [ 10 ] , 

/*  V 

fltbuff [15) [30] ; 

/********** *************** 

int 

nnlinear , 

/ *  Number  of  linear  factors  in  numerator 

V 

nnquadratic , 

/*  Number  of  quadratic  factors  in  numerator  */ 

dnlinear , 

/»  Number  of  linear  factors  in  denominator  V 

dnquadratic , 

/*  Number  of  quadratic  factors  in  denominator  */ 

total , 

/*  Number  of  input  window  fields  needed  */ 

field-0, 

/*  Current  input  field  number  */ 

col-0 , 

/*  Current  input  column  number  */ 

jl.j2; 

/*  Input  loop  counters  */ 

static  char  *msg  -  "[Press  ENTER  to  Accept,  any  other  key  to  Edit]1'; 


nstring[0]=' \0 ' ; 
dstring[0)-'\0' ; 

batrib  -  v_setatr (WHITE, BLACK , 0 , 0 ) ; 
watrib  -  v_setatr (WHITE , BLACK , 0 , 0) ; 

w2  -  wn_open(0, 2, 12, 52, 10 .watrib .batrib) ; 
if(!w2)  exit(l); 

wn_title(w2, "  Factored  Input  Menu 


wn_printf(w2,"\n  Examples  of  a  linear  factor:  (s  +  0),(s  +  5)  \n"); 

wn_printf (w2,"\n  Example  of  a  quadratic  factor:  (s*2  +  5s  +  10)  \n\n" ) ; 

wn_printf (w2, " - - - \n" ) ; 

wn_printf (w2,"  Number  of  linear  factors  in  the  numerator  :  \n"); 

wn_printf(w2,"  Number  of  quadratic  factors  in  the  numerator  :  \n” ) ; 

wn_printf  (w2, "  Number  of  linear  factors  in  the  denominator  :  \n'*); 

wn_printf(w2,”  Number  of  quadratic  factors  in  the  denominator  :  "); 


f rm-wn_f raopn(  5 ) ; 
iff !frm)  exit(0); 
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START : 

nlbufftO]  -  0; 
nqbuff [0]  =  0; 
dlbuff [0]  «  0; 
dqbuff [ 0 !  -  0; 

wn_gint(SET, frm,0 ,w2, 6, 50 , NSTR, watrib , ' O' .Stnnlinear ,1,0, 9 , nlbuff , NSTR, NSTR) ; 
wn_gint(SET, f rm, 1 ,w2 , 7 , 50 , NSTR, watrib , ’ O' ,&nnquadratic ,1,0, 9 , nqbuff , NSTR , NSTR ) ; 
wn_gint(SET , f rm, 2 ,w2 , 8 , 50 , NSTR , watrib , ’ 0  * ,&dnlinear ,1,0 , 9, dlbuff , NSTR, NSTR) ; 
wn_gint(SET, frm, 3 ,w2, 9, 50 , NSTR , watrib , ' O' ,&dnquadratic ,1,0 , 9 , dqbuff , NSTR, NSTR) ; 


rv-wn_frmget(frm) ; 
i£(rv=0)  return(O); 

total=nnlinear+dnlinear+2*(nnquadratic+dnquadratic ) ; 
nlength=l+13*nnlinear+25*nnquadratic ; 
dlength=13*dnlinear+25*dnquadratic ; 

if (nlength>dlength) 

{ 

maxchar*>nlength; 
skipn=2 ; 

skipd=(nlength-dlength)/2+3 ; 

) 

else 

{ 

maxchar=dlength; 

skipn=(dlength-nlength)/2+2; 

skipd=2; 

} 

batrib  ■  v_setatr (RED .WHITE , 0 , 0 ) ; 
watrib  “  v_setatr (RED .WHITE , 0 , 0 ) ; 

if (maxchar>78) 

{ 

if (nlength>78) 

{ 

w3  *  wnopent 0 , 15 , 5, 68, 1 .watrib .batrib ) ; 
if(!w3)  exit(l); 
wn_title(w3 , "  Warning  "); 

wnprintf (w3 , "  The  numerator  is  too  large,  try  again  with  a  lower  order 
equation. " ) ; 

> 

else 

{ 

w3  “  wn_open( 0 , 15 , 4 , 70 , 1 .watrib , batrib) ; 
i£(!w3)  exit(l); 
wn_title(w3 , "  Warning  "); 

wn_print£(w3 , "  The  denominator  is  too  large,  try  again  with  a  lower  order 
equation. ; 

} 

wn-wn_open( 1000 , 2* , 0 . 27 , 1 , (RVIDEO) .NVIDEO) ; 
if(!wn)  exit(l); 

wnjprintf (wn, " [Press  any  key  to  continue] \a" ) ;  /*  display  message  */ 

input=NUL; 

v_kf lushC ) ; 

input“getch( ) ; 

wn_cloae(wn) ; 

wn_close(w3) ;  /*  make  message  go  away  */ 

batrib  ■  v_setatr(WHITE, BLACK, 0,0) ; 
watrib  -  v_aetatr(WHITE,BLACX,0,D); 
goto  START; 

> 

m“nnlinear+2*nnquadratic ; 
n»dnlinear+2*dnquadratic ; 

batrib  -  v_aetatr (MAGENTA , WHITE. 0, 0) ; 
watrib  "  vie t a tr (MAGENTA, WHITE, 0,0) ; 
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w3  -  wn_open(0 , 9,  ( 37-maxchar*0 . 5) ,  ( 4+maxchar ) ,  5, watrib  .batrib) ; 
if(!w3)  exit(l); 

wn_title(w3 , "  Input  Data  "); 

wn_locate(w3 , 1 ,  skipn) ; 
wn_printf (w3 , "K" ) ; 

for( j=l ; j<=nnl inear ;j++) 

wn_printf(w3, "(s+  )"); 

for tj=l ; j<=nnquadratic ;j++ ) 

wn_printf (w3 , " ( s”2+  s+  )“); 

wn_locate(w3 ,2,2); 
fort j=l ; j<=tmaxchar+l ) ; j++) 
wn_printf (w3 , ; 
wn_locate(w3 , 3 , skipd) ; 

for(j=l; j<=dnlinear; j++) 

wn_printf (w3 , "(s+  )"); 

fort j=l;  j <=dnquadra tic; j++) 

wn_printf (w3 , " t s~2+  s+  )"); 

frm2=wn_frmopn(total+l) ; 
if(!frm2)  exittO); 

1=0; 

for''  j=0  ; j<=tnnlinear-l) ;j++) 

{ 

f ltbuff [fie ’ d ] ( 0 ]=0 ; 
col=skipn+4+j*13 ; 

2coef [ j ]=0 . 0 ; 

wngf loattSET , frm2 , f ield ,w3 , 1 , col , NSTR, watrib , 1  1 ,&zcoef [ 1++ ] , 9 , A , 
-10000. 0, 10000. O.fltbuff [field) , NSTR, NSTR) ; 
field++; 

) 

f  ;k<=tnnquadratic-l)  ;k++) 

{ 

fltbuff [field] [0]=0; 
col=skipn+13*nnlinear+6+26*k ; 
zcoef [nnlinear+k 1=0.0; 
tcoef [nnlinear+l+k]=0 . 0 ; 

wn_gf Itat [SET, f rm2 ,  f  ield  ,w3 , 1 , col , NSTR ,  watri1' .  '  ‘  , &zcoef  [1++] , 9, 4 , 
-10000.0, 10000. C,f ltbuff [fi,:d] .NSTR, NSTR); 
f ield++ ; 

fltbuff [f  eld) l 0 ] -0 ; 
col=col+l.l  ; 

wn_gf loattSET,  f  i2 .  neld,w3 ,  l.col , NSTR,  watrib,  '  '  ,&zcoef  [1++]  ,9,4, 
- 10000 . 0 , 10000 . 0 , fltbuff! field] , NSTR , NSTR) ; 
f ield++ ; 

) 

1=0; 

fort j=0 ; j<=(dnli„oar-l) ; j++) 

< 

fltbuff [field] [ 0 ] =0 ; 
col=skipd+3+j*13 ; 
pcoef [ j ]=0 . 0 ; 

wn_gf loattSET, frm2 , field, w3 ,3, col  NfTR .watrib , '  ' .fcpcoef [ 1++) . 9 , 4 , 
-10000. 0,10000. O.fltbuff [field] , NSTR, NSTR), 
f ield++ ; 

> 

for (k=0 ;k<=(dnquadratic-l ) ;kt+) 

{ 

fltbuff [field] [ 0 ] =0 ; 
eol=skipd+13*dnlinear+5+26*k ; 
pcoef [ dnlinear+k ] =0 . 0 ; 
pcoef [dnlinear+k+l]=0 . 0 ; 

wn_gf loattSET. frm2 , f ield,w3 , 3 , col , NSTR, watrib, '  ' .ipcoef [ 1++] ,9,4 , 
-10000 . 0 , 10000 . 0 , fltbuff [field] ,NSTR, NSTR) ; 
f ield++; 

fltbuff [field] (0]=0; 
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col=col+ll; 

wn_gf loattSET, fnn2, f ield,w3 , 3 , col ,NSTR,watrib, '  ' ,&pcoef [1++] , 9, A , 

- 10000.0, 10000, O.fltbuff [ field], NSTR.NSTR); 
field++; 

} 

input=NUL ; 
while (input  !=  CR) 

{ 

wn_£rmget2 ( f rm2 ) ; 

nlength=dlength=0 ; 

nlength  =  sprintf (nstring , "K" ) ; 

1-0; 
k=0 ; 

fort jl=0;jl<=(nnlinear-l) ; jl++) 

{ 

if  (zcoef  [11=0 , 0 ) 

nlength  +=sprintf (nstring+nlength, " (s)" ) ; 
else 

nlength  +*sprintf (nstring+nlength, " (st+g)" , zcoef [1] 1 ; 
zerosfk] [ l]=-zcoe£[l++] ; 
zeroslk] [21=0.0; 
k++; 

} 

for (j 2=0 ; j2<=(nnquadratic-l) ;  j2++) 

{ 

nlength  +=sprintf (nstring+nlength, " ( s' 2X+gsX+gl" , zcoef [1] , zcoef [1+11 ) ; 

a=1.0; 

b=zcoef [1++1 ; 
c=zcoef [ 1++1 ; 
d=b*b-4*a*c ; 
if (d>0 . 0 ) 

{ 

d=sqrt(d) ; 

zeros [k] [ 1 1 ”( -b+d)/ (2*a) ; 
zeros [k++] [21=0.0; 
zeroslk]  [  1}*°( -b-d)/ (2*a) ; 
zeros [k++] [21=0.0; 

} 

else 

{ 

d=sqrt(-d) ; 
zeros [k] [ 1] =-b/ (2*a) ; 
zeros[k++l [2]=d/(2*al ; 
zeros [k]  [  11— b/(2*a) ; 
zeros[k++]  [21—d/ (2*a) ; 

} 

} 

1=0; 

k=0; 

f or ( j 1*0 ; j l<=(dnlinear- 1 ) ; jl++) 

{ 

if (pcoef [1]==0 . 0) 

dlength  +=sprint£(dstring+dlength, " (s)" ) ; 
else 

dlength  +=sprintf (dstring+dlength, "(sl+g)" , pcoef [1] ) ; 
poles[k) [l]=-pcoef [1++1 ; 
poles[k++) [21=0.0; 

) 

for (j 2=0 ; J2<“(dnquadratic-1) ; j2++) 

{ 

dlength  +=sprintf(dstring+dlength,"(s*2Z+gsX+g)", pcoef [1] .pcoef [1+11 ) ; 
a=l . 0 ; 

b-pcoef [1++] ; 
c=pcoef [1++1 ; 
d=b*b-A*a+c ; 
if(d>0.0) 

{ 

d=aqrt(d) ; 

poles  tkJ(U=(-b+d)/(2*a); 
poles[k++) [21=0.0; 
poles [kj (l]=(-b-d)/(2*a) ; 
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poles  [k++]  [2 J=0 . 0 ; 

> 

else 

{ 

d-sqrt(-d) ; 

poles (k  H 1 1— b/  (2*a)  ; 

poles[k++][21«d/(2*a); 
poles [k] [l]=-b/(2*a) ; 
poles [k++]  [21—  d/(2*a); 

} 

) 

if (nlength>=dlength) 

{ 

raaxchar=nlength ; 
skipn=2 ; 

skipd=(nlength-dlength)/2+2; 

} 

else 

{ 

maxchar=dlength ; 
skipn=(dlength-nlength)/2+2; 
skipd=2 ; 

> 

batrib  -  v_setatr (RED , WHITE, 0 , 0 ) ; 
watrib  =  v_setatr (RED , WHITE, 0, 0 ) ; 

w4  «  wnopentO, 16, ( 37-maxchar*0 . 5) , (4+maxchar) , 5, watrib, batrib) ; 
if(!w4)  exit(l); 

wn_locate(w4 , 1 , skipn) ; 
wnprintf (w4 , "Zs" .nstring ) ; 
wnlocateCwA .2,2) ; 
for (j”l ; j<^naxchar : j++) 
wn_printf  (w4 , ) ; 
wn_locate(w4 , 3, skipd) ; 
wn_printf(w4 , “Zs" .dstring) ; 

vm-wn_open( 1000 ,24 ,0 , (int)strlen(msg) , 1 , CRVIDEO) .NVIDEO) ; 
if ( !wn)  exit(l); 

wn_puts (wn ,0,0 ,  msg) :  /*  display  message  */ 

input=NUL; 

vkflushO; 

input=*getch( ) ; 

wnclose(wn) ;  /*  ®ake  message  go  away  */ 

) 

wn_frmcls(frm) ; 
wn_frmcls(frm2) ; 
wn_close(w2) ; 
wn_close(w3) ; 
wn_close(w4) ; 

/************************  END  of  factored  routine  *************************/ 
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/**********„*,*»*****,  NON-FACTORED  POLYMONIAL  PROGRAM 


This  program  prompts  the  user  for  the  transfer  function  information. 
This  transfer  function  is  in  the  non-factored  polynomial  form. 


**/ 

*/ 

*/ 

V 

*/ 


/*..***»**»***«**,**»**.*  START  GLOBAL  VARIABLES  *****»****************»*»*/ 


#include  "variables .h"  /*  Include  variable  list  */ 

/***»**»,**.»**,****.»*.*»,  jjjjj  qloBAL  VARIABLES  ***************»**********/ 
/****«*.*.*.«**.»,*»„  START  OF  NON-FACTORED  ROUTINE  ***********»*****»*****/ 
int  coeff() 

{ 

int  j,k, 

rv, 

start; 


float 

at  12]  , 

b 

12]; 

WINDOWFTR 

w2 ; 

/*  Order  of  polynomial  window  handle 

*/ 

WINDOWPTR 

w3  ; 

/*  Numerator  window  handle 

*/ 

windowptr 

w4  ; 

/*  Denominator  window  handle 

*/ 

WINDOWPTR 

wn; 

/*  Message  window  handle 

*/ 

WIFORM 

frm2; 

/*  Order  of  polynomial  form  handle 

*/ 

WIFORM 

frm3 ; 

/*  Numerator  form  handle 

*/ 

WIFORM 

frraA  ; 

/*  Denominator  form  handle 

*/ 

int 

batrib ; 

/*  Border  atrib 

*/ 

int 

watrib ; 

/*  Window  atrib 

*/ 

register 

i; 

char 

mbuff [10] , 

/*  Denominator  variable  buffer 

*/ 

nbuff [  10  ]  ,  /*  Numerator  variable  buffer  */ 

fltbuff [30] [30]  ;  /*  Float  variable  buffer  array  */ 


static  char  *msg  =  "[Press  ENTER  to  Accept,  any  other  key  to  Edit]"; 
static  char  *errmsg  -  "[Order  of  polyinomial  must  be  <»  9]"; 

nstring[0]=*\0' ; 
dstring [0] \0 1 ; 


batrib  -  v_setatr (CYAN , BLACK , 0 , 0 ) ; 
watrib  =  vaetatr (CYAN , BLACK, 0 , 0) ; 


/*  Set  border  atrib  */ 
/*  Set  window  atrib  */ 


w2  -  wn_open(0 , 2, 15, 46 , 5, watrib .batrib) ; 
if(!w2)  exit(l); 


/*  Open  order  window  */ 
/*  Exit  if  unable  */ 


wn_title(w2 , "  Polynomial  Input  Menu  "); 


/*  Order  window  title  */ 


wn_printf(w2,"\n  Cls'm  +  C2s'(m-1)  +  C3s*(m-2)  ...  C4s  +  C5\n"); 

wn_jprintf  (w2 , " - - - \n” ) ; 

wn_printf(w2, ''  Highest  order  of  the  numerator  :  \n"); 

wnjprintf (w2,"  Highest  order  of  the  denominator  : 

mbuff [0]  ■  0; 
nbuff [0]  -  0; 


frm2*wn_frmopn(3) ; 
if(!frm2)  exit(0); 


wn_gint(SET, frm2,0,w2, 3, 42.NSTR, watrib, ’0’ ,Sm, 1,0, 9, mbuff ,NSTR, errmsg) ; 
wn_gint(SET , frm2 , 1 ,w2, * , 42, NSTR .watrib ,  ’  0  ’  ,ta ,  1 , 0 , 9 , nbuff ,NSTR, errmsg ) ; 

rv=wn_irmget(frm2) ; 
if(rv—0)  retum(O); 

i£(m>-l) 

{ 

batrib  -  v_setatr (GREEN, BLACK. 0,0) ; 
watrib  -  v~setatr(GREEN, BUCK. 0,0); 


*•3  -  wn_open(0, 10, 10, 23, m+1, watrib, batrib); 
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if(!w3)  exit(l); 
wn_title(w3 , "  Numerator  "); 

fnE3^*n_£nnopn(n>+2) ; 
if ( ! f rm3 )  exit(O); 

for(  j*0;  j<=^n;  j++) 

{ 

wn_locate(w3, j ,3) ; 
if ( j”=(m-l ) ) 

wn_printf (w3 , " (  )s"); 

else  if(j“=m) 

wn_pnntf  (w3 ,  "  (  }"); 

else 

wn_printf(w3  , ''  (  )s*td",m~j); 

fltbuff [j] [01=0; 
a[j)«0.0; 

wn_gf loat(SET , f rm3 , j ,w3 , j , A , NSTR , watrib , *  ',&a[j],ll,A, 
-1000000.0, 1000000. 0, fltbuff ( j] , NSTR, NSTR) ; 

} 

) 

else  if(m“0) 
a(0)-1.0; 

batrib  «  v_setatr (MAGENTA , BLACK, 0 , 0 ) ; 
watrib  =  v_setatr (MAGENTA, BLACK , 0 , 0 ) ; 

wA  *  wn_open(0 , 10 , AS , 23 ,n+l .watrib , batrib ) ; 

if(!wA)  exitd); 

wn_title(wA , "  Denominator  "); 

£rmA-wn_frmopn(n+2) ; 
if(!frmA)  exit(0); 

for(k-0;k<«n;k++) 

{ 

wn_locate(wA ,k , 3 ) ; 
if(k—(n-l)) 
wn_printf(wA , "( 
else  if(k-n) 
wnjprintf (wA , "( 
else 

wn_printf (wA , "( 

fltbuff [j+k] [0)-0; 
b[k]-0.0; 

wngf loat(SET, frmA , k , wA ,k , A , NSTR, watrib, '  ',&bIk),ll,A, 

- 1000000 .0,1000000.0, fltbuff! j+k ) , NSTR , NSTR) ; 

} 

i£(m>0) 

< 

wn_activate(w3) ; 

batrib  ■  v_setatr (GREEN , BLACK, 0 , 0 ) ; 
watrib  •  v_»etatr (GREEN, BLACK, 0 , 0) ; 
rv-wn_frmget(£rm3) ; 
if(rv—0)  retum(O); 

) 

wn_aetivate(wA) ; 

batrib  “  v_eetatr (MAGENTA , BLACK , 0 , 0 ) ; 
watrib  ■  v_eetatr (MAGENTA .BLACK, 0 ,0) ; 
rv*wn_fnnget(frmA ) ; 
if(rv-“0)  retum(O); 

nlength“0 ; 

nlength  +-eprintf (netring+nlength, "K(") ; 
for(  j-0;  J<^n;  j++) 

{ 

lf(a[j]—  0.0) 
continue; 

elee  if(a[j]— 1.0  BA  J<(m-1)  BA  J!-0) 


)  5  " ) ; 

)"); 

)s*Id" ,n-k ) ; 
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nlength  +=sprintf  (nstring+nlength,  "+s*Xd" ,m-j ) ; 
else  if(a[j)“1.0  &&  j“0) 

nlength  +*sprintf (nstring+nlength , "s*Xd" ,m- j ) ; 
else  if(j=^n) 

nlength  +-sprintf (nstring+nlength, "X+g" , a[j ] ) ; 
else  if(j=*ro-l  &&  a[j]!“1.0) 

nlength  +»sprintf (nstring+nlength, "X+gs" , a[ j ] ) ; 
else  if(j=^n-l  &&  a(j)“1.0) 
nlength  +=sprintf (nstring+nlength, "+s" ) ; 
else 

nlength  +=sprintf  (nstring+nlength.  "X+gs" Xd’1 ,  a!  j  ]  ,m-j  ) ; 

} 

nlength  +=sprint£( nstring+nlength, " )" ) ; 
dlength=0 ; 


for ( j=0 ; j<=n ; J++ ) 

{ 

if(b[j]—0.0) 
continue ; 

else  if (b [ j ] “1 . 0  &&  j“0) 

dlength  +-sprintf (dstring+dlength, "s*Xd" ,n-j ) ; 
else  if (bt j  J— — 1  -  0  &&  j<(n-l)  &&  j!-0) 
dlength  +=sprintf (dstring+dlength , "+s"Xd" ,n-j ) ; 
else  if(j=“n) 

dlength  +-sprintf (dstring+dlength, "X+g" , b [ j ] > ; 
else  if(j—n-l  &&  btj)!-1.0) 

dlength  +-sprintf (dstring+dlength, "X+gs" ,b( j ] ) ; 
else  if ( j=n-l  &&  bfjl”1.0) 
dlength  +-sprintf (dstring+dlength, "+s” ) ; 
else 

dlength  +“sprintf (dstring+dlength, "X+gs*Xd" ,b[ j ) ,n-j ) ; 

> 

i f  ( n  1 eng  th>“d 1 eng th ) 

{ 

maxchar-nlength; 
skipn*2 ; 

skipd=(nlength-dlength)/2+2 ; 

) 

else 

{ 

maxchar-dlength; 

sk ipn“( dlength -nlength ) /2+2 ; 

skipd”2; 

) 


fortj-O; j<-m; j++) 
zcoef [ j+lj=a[ j] ; 

£or(j“0; j<”n; j++) 
pcoef [ j+1 ] “b  t j ] ; 

return; 

} 

/*********»*******»*****  qjxi  OF  NON-FACTORED  ROUTINE  ********•*************/ 
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/***•***************  START  OF  VIDEO  SUPPORT  FUNCTIONS  ***»*****************/ 


/*  */ 

/*  This  files  includes  the  source  for  the  following  video  funcions:  */ 

/*  */ 

/*  int  set_video_tnode( ) :  set  video  mode  to  EGA  or  CGA  */ 

/*  void  arrow(double,  double,  double,  int,  double):  plot  arrow  */ 

/*  */ 

j  *****  ********  *  ********  *********  *  *  *  ****************  **************************/ 

/it************************  START  GLOBAL  VARIABLES  ********»**********»*****/ 


♦include  "variables .h"  /*  Include  variable  list  */ 

/***«**»***«*»*************  end  GLOBAL  VARIABLES  **************************/ 


/a****************),*******  SET  VIDEO  MODE  PROGRAM  *************************/ 
/*  */ 


/*  This  routine  sets  the  video  screen  to  the  EGA  or  CGA  video  mode.  */ 
/*  It  first  attempts  to  set  the  video  mode  to  EGA,  if  this  adapter  */ 
/*  is  not  supported  it  then  trys  to  set  it  to  CGA.  If  CGA  is  not  */ 
/*  supported  either,  the  program  exits  with  an  error  message.  */ 
/»  */ 


/***.*«****.*..***«•»**»  START  OF  SET  VIDEO  ROUTINE  ***********************/ 
int  set_video_mode( ) 


{ 

iff  setvideomode (  ERESCOLOR))  /*  Try  to  set  EGA  */ 


{ 


forcolorl^ll ; 
forcolor2*2 : 
return (  ERESCOLOR ) ; 
} 


if (_setvideomode(_HRESBW) )  /*Try  to  set  CGA  BLACK  &  WHITE*/ 

T 

forcolorl*3 ; 
forcolor2*2 ; 
retum(_HRESBW) ; 

) 

else 

return(O) ; 

) 

/no**********************  end  OF  SET  VIDEO  ROUTINE  ****»*******************/ 


/***********************  PLOT  INFINITY  ARROW  PROGRAM  ************************/ 
/* 

/» 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/* 

/..***..*.****.*.*. 

♦include  <math.h> 

♦include  <graph.h> 


This  routine  plots  an  arrow  on  the  screen  to  indicate  infinity. 
A  description  of  the  passing  variables  follows. 

xl  *  x  starting  point  of  arrow  tip. 
yl  “  y  starting  point  of  arrow  tip. 
slope  “  exit  angle  of  arrow. 

color  m  color  that  the  arrow  is  to  be  plotted  with. 
aspect_ratio  ■  axspect  ratio  of  the  video  adapter . 


*/ 

*/ 

*/ 

*/ 

*/ 
*/ 
*  / 

V 

*/ 

V 


START  OF  MAIN  ARROW  ROUTINE  **********************/ 


/oid 


arrowtdouble  xl, double  yl, double  slope, int  color, double  aspect_ratio) 


double  x2,x3,y2,y3; 

_moveto(xl,yl) ; 

x2-xl+15*eos(3 . l*159265+slope+0 . 2618 ) ; 

y2-y 1- ( 15*sin(3 . l*159Z65+slope+0 . 2618 ) )*aspect_ratio ; 

_lineto(x2,y2) ; 

x3“xl+15*cos( 3 . 1* 159265+slope-0 . 2618 ) ; 

y3-yl-(  l'i*sln(3.  lA159265+slope-0.2618))*aspect_ratio; 
_llneto(x3,y3) ; 

_lineto(xl,yl) ; 

x2-xl+10*cos(3.  H159265+*lope+0. 1309) ; 
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y2=y l-( 10*sin(3 . 14159255+slope+0 . 1309) )*aspect_ratio ; 

_f loodf ill ( x2 ,y2, color ) ; 
x2-xl+10*cos(3  .  1 **  1 5926  5+s  1  ope -0 . 1309) ; 
y2-yl-(10*sin(3 . 14159265+slope-0 . 1309) )*aspect_ratio ; 

_£loodfill(x2,y2, color) ; 

x2-xl+12.5*cos(3.  U159265+slope+0. 1309) ; 

y2=y  1-  ( 12. 5*sin(3 .  14159265+slope+O . 1309 ) )*aspect_ratio ; 

_floodfill(x2,y2, color) ; 

x2-xl+12. 5*cos(3  .  14 15926 5+s lope -0 . 1309) ; 

y2=yl- ( 12. 5*sin(3 . 14159265+slope-0 .1309) )*aspect_ratio ; 

_f loodf Ill(x2,y2 , color) ; 

) 

/**********★************  pjjo  Qp  pix)T  ARROW  ROUTINE  ************************/ 
/****■»«,.•«■»*******»< ,**  ejjq  of  VIDEO  SUPPORT  FUNCTIONS  »»»***•**»»******»****/ 


APPENDIX  I 


ROOT  SOLVING  SOURCE  CODE 
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/*************************  KOOT  FINDING  PROGRAM  **************************/ 


/*  V 
/*  This  program  uses  the  Bairstow's  Method  to  extract  the  */ 
/*  the  real  and/or  complex  roots  of  a  polynomial  equation  */ 
/*  of  the  form  */ 
/*  */ 
/*  A(l)X'tn)  +  A(2)X'(n-l)  +  ...  +A(n+1)  =0  */ 
/*  */ 


/it************************  START  GLOBAL  VARIABLES  it************************/ 

^include  "variables .h"  /*  Include  variable  list  */ 

/**************************  Ejjd  GLOBAL  VARIABLES  **************************/ 
/**************************  START  OF  ROOT  ROUTINE  *************************/ 
int  ROOTtfloat  *A, double  *RT,int  N) 

{ 

double  B[l«], 

criA] . 

D, 

E, 

F, 

P-0.0, 

PI, 

q=o.o, 

01, 

DELTP , 

DELTQ, 

EPS=0. 0000001, 

DEN, 

CBARL, 

SUM, 

SUM1, 

INITIAL[6] ; 

int  I , 

J, 

L, 

M, 

K, 

IITAG-1000, 

ITAG=0; 


fort  1=0 ;I<=13;I++) 
B[I]-C[IJ-0.0; 

INITIAL! 11-0.0; 
INITIAL[2]-1.0; 
INITIAL! 3] -10 . 0 ; 
INITIAL [4  J -A [2 J/2 . 0 ; 
INITIAL [5] -A [3]/2.0; 

D-A [ 1 ] ; 


for ( I-l ; I<-N ;  I++ ) 
A[I)-A[I+1)/D; 
for (K-l ;K<-1 0 ;K++) 

{ 

LAO : 

if(N>l)  goto  LA3; 

ITAG-ITAG+1; 
L-ITAG+ITAG; 
*(RT+L)-0 . 0 ; 
*(RT+L-1 )— At  1 1  ; 
return (A ) ; 

LA3 : 

if (N>2)  goto  LA6 ; 


P-Alll; 
0-A(2] ; 
goto  L8 ; 
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U6: 

P=INITIAL[K] ; 

Q=0 . 0 ; 

M=l; 

L51: 

B [ 1 ] “A [ 1 ] -  P ; 
B[2]=A[2]-P*B[1]-Q; 

L=N - 1 ; 

C[l]-B[l]-P; 

C[2]-B[2] -P*C[ X] “0; 

if(L=2)  L— 3  ; 

for  ( J“3  ;  J<=L ;  J++ ) 

{ 

B[J]-A[J]-P*B[J-l]-Q*BfJ-2] : 
C[J]=B[J]-P*C[J-l]-Q*C[J-2] ; 
} 


CBARL=C[L] -B[L] ; 
DEN=-CBARL ; 

if (N ! =3 )  DEN=DEN*C [N-3] ; 

DEN=DEN+C(N-21*C[N-2] : 

if (DEN ! *0 . 0 )  goto  LI; 
P*P+1.0; 

0=0+1. 0; 
goto  L5 1 ; 


B[N]-A[N]-P+B[N-1)-Q*B[N-21 ; 
DELTP=-B[N] ; 

if (N!*3)  DELTP=DELTP*C1;N-31  ; 

DELTP=(BtN-l]*C(N-2]+  DELTP)/DEN ; 
DELTQ-(BtN]*C[N-2]-B[N-l]*CBARL)/DEN; 
P=P+DELTP ; 

Q-Q+DELTQ; 

SUM»f abs ( DELTF ) +f abs ( DELTQ ) ; 

if CM==1 )  SUM1  -  SUM; 
if(M!“5| |SUM<“SUM1)  goto  Lll; 

LI  1 : 

if (SUM<=EPS)  goto  L8 ; 

if (M<IITAG)  goto  L2; 

goto  L100; 

L2: 

M-M+l; 
goto  L51; 


D”-P*0 . 5; 
ITAG-ITAG+2 ; 
F-Q-D+D; 

E”pow(£abs (F) ,0.5); 
L-ITAG+ITAG; 

if CF>0 , 0)  goto  L31; 

*(RT+L)-0.0; 
*(RT+L-1)-D-E; 
*(RT+L-2)-0.0; 
*(RT+L-3)“D+E; 
goto  L32; 


L31 : 


*(RT+L)— E; 
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*  CRT+L- 1)-D ; 

*(RT+L-2)=E; 

*  CRT+L-3 )-D ; 

L32: 

N-N-2; 

if (N>0)  goto  LSI; 
return  ('ll ; 

L81 : 

for ( 1*1 ; I<=N ; I++ ) 

A [ 1 1 =B [ I ]  ; 

goto  L**0 ; 

} 

L100: 

K-K; 

)***.*«*»**<.**»****»««**»*  end  OF  ROOT  ROUTINE  ** ************************** 
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/****.*«*•»•***«**  START  OF  PRINTER  SUPPORT  FUNCTIONS  *********************/ 


/*  */ 

/*  This  files  includes  the  source  for  the  following  printer  funcions:  */ 

/*  */ 

/*  int  set_print_screen( ) :  loads  print  screen  program  */ 

/*  int  printer!):  prompt  user  for  print  screen  */ 

/*  */ 


/*********************  ******  *************************************************  / 
/**•»*******•»•******»»*»*  START  GLOBAL  VARIABLES  *************»*****»*»»**/ 


#include  "variables . h"  /*  Include  variable  list  */ 
/**************************  GLOBAL  VARIABLES  **************************/ 
/*****«*******«*********  LOAD  PRINT  SCREEN  PROGRAM  *****«******************/ 


/*  */ 

/*  This  routine  loads  the  proper  print  screen  program  into  memory.  */ 
/*  It  first  attempts  to  load  the  EGA  to  EPSON  print  screen  program,  */ 
/*  if  this  video  adapter  is  supported  it  trys  to  load  the  CGA  print  */ 
/*  screen  program.  */ 
/*  */ 


/*******.»***„««*.**  START  OF  SET  PRINT  SCREEN  ROUTINE  ******»************»/ 
int  set_print_screen( ) 


{ 


int  rv;  /*  Return  value  */ 

if (_setvideomode(_ERESCOLOR) )  /*  Try  to  load  EGA  */ 

{ 

rv=sy stem! "egaepson" ) ; 
return! 0) ; 

} 


else  if!  setvideomode!_HRESBW) )  /*  Try  to  load  CGA  BLACK  &  WHITE*/ 

{ 

rv=system! "graphics" ) ; 
return! 0) ; 

) 

else 

{ 

_clearscreen!_GCLEARSCR£EN) ;  /*  Video  is  not  supported  */ 

printf ("\n\n\n  This  video  adapter  is  not  supported"); 

exit!0) ; 

} 

} 

/*********************  Qjp  QP  SET  print  SCREEN  ROUTINE  a*******************/ 


/*************»**»*„**  START  PRINTER  OUTPUT  ROUTINE  a*********************/ 
/*  */ 
/*  This  routine  prompts  the  user  once  the  program  is  finished  plotting  */ 
/*  to  the  screen.  The  program  sends  the  screen  image  to  the  printer  */ 
/*  and  returns  to  the  main  menu  if  the  user  presses  the  P  key,  else  */ 
/*  any  other  key  returns  control  back  to  the  main  menu  */ 
/*  */ 


/****************************************************************************/ 

int  printer!) 

{ 


int  rv; 


/*  Return  value  */ 


fflush(stdin) ;  /*  Flush  keyboard  buffer  */ 

rv-getch();  /*  Prompt  user  for  key  */ 

if!  rv*"“ '  P '  ||  rv” '  p ' ) 

{ 

print_screen! ) ;  /*  Screen  dump  to  printer  */ 

return! 1) ; 

) 

else 

return! 0) ; 

> 

/***•««***»***•*»***,.***  END  PRINTER  OUTPUT  ROUTINE  •*•*****••••••****•***/ 
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/**«********•*********  START  PRINTER  INTERRUPT  ROUTINE  «*****»*******«*****/ 
/*  */ 
/*  This  routine  sends  the  contends  needed  to  do  s  screen  dump  to  the  */ 
/*  printer.  */ 


/*  V 

/**************************************************************************** / 
void  print_screen( ) 

{ 

union  REGS  registers; 

mt86(5, Reregisters  .Reregisters) ;  /*  Send  interrupt  */ 

fprintf(stdprn,''\f") ;  /*  Form  feed  »/ 


> 


/**•******************«  end  PRINTER  INTERRUPT  ROUTINE  **************«******/ 


/ 


END  OF  PRINTER  SUPPORT  FUNCTIONS  **************«****/ 
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/*********************  START  OF  GLOBAL  VARIABLE  LIST  *********•************/ 
/*  */ 

/*  This  is  the  list  of  global  variables  used  throughout  the  CSS P  V 

/*  routines.  */ 

/*  V 

/****.»*»********«...***  START  GLOBAL  INCLUDE  FILES  ***************»«*«****/ 


♦include  <stdio.h> 

♦include  <graph.h> 

♦include  <math.h> 

♦include  <malloc.h> 

♦include  <dos.h> 

♦include  "windows .h" 

/***********************»  end  GLOBAL  INCLUDE  FILES  **************»*********/ 


/**************************  START  GLOBAL  ROUTINES  *************************/ 


int  factoredC), 

wn_frmget2( struct  wi_scb  *  *frm), 
input ( ) , 
coeff ( ) . 

set_video  mode(void), 

ROOT (float  *, double  *,int), 
printer ( ) ; 


/*  Factored  polynomial  routine  */ 
/*  Modified  formget  routine  */ 
/*  Input  routine  */ 
/*  Coefficients  input  routine  */ 
/*  Set  video  mode  routine  */ 
/*  Root  solving  routine  */ 
/*  Printer  output  routine  */ 


void  output_screen( ) .  /*  Output  screen  routine  */ 

plotCint,  int),  /*  Root  locus  routine  V 

arrowCdouble,  double,  double,  int,  double),  /*  Arrow  routine  V 

print_screen( ) ;  /*  Print  screen  routine  */ 

/******•*********»**********  END  GLOBAL  ROUTINES  *******»*****»»*»**»»*»**»/ 


/**************************  START  GLOBAL  VARIABLES  ****•**************»****/ 

struct  videoconfig  vc ;  /*  Video  configuration  structer  V 


char  *greenzero, 
*redzero, 
•grrenpole, 
*redpole , 
nstring[90) , 
dstring[90] , 
string [90]  , 
atringg[3] ; 


/*  Character  pointers  */ 


/* 

V 

/* 

V 

/* 

•/ 

/*  Numerator  string  array  */ 
/*  Denominator  string  array  */ 
/»  Temporary  string  array  */ 
/*  Temporary  string  array  */ 


int  n, 

m, 

maxchar , 

skipn, 

skipd. 

forcolorl , 

forcolor2, 

dlength, 

nlength; 


/*  Number  of  numerator  roots  */ 
/*  Number  of  denominator  roots  */ 

/*  Maximum  number  of  charaters  */ 

/*  Skip  n  number  of  spaces  */ 

/*  Skip  d  number  of  spaces  */ 

/*  Forground  color  ^1  V 
/*  Forground  color  */ 

/*  Length  of  denominator  string  */ 

/*  Length  of  numerator  string  */ 


double  poles[12J [2] , 
zeros! 12] [2] , 
delta  s, 

K; 


/*  Complex  pole  roots  array  */ 
/*  Complex  zero  roots  array  */ 

/»  S  step  size  */ 

/*  Gain  value  */ 


float  pcoef[20], 
zcoef (20 J , 
aspect_ratio, 
max, 
min; 

/***.****•••**.*•.< 


/*  Pole  coefficients  array  */ 
/*  Zero  coefficients  array  V 
/*  Screen  aspect  ratio  */ 

/*  Maximum  plot  value  */ 

/*  Minimum  plot  value  */ 

END  GLOBAL  VARIABLES  ***********»*»»**•*»»**»»/ 
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/**..***.*•****».«*  STAST  OF  MODIFIED  WINDOW  ROUTINES  «********************/ 

/*  V 
/*  This  file  is  a  modified  version  of  the  windows. h  file  supplied  */ 
/*  in  the  windowing  package  */ 
/*  */ 
/*  The  Window  BOSS  */ 
/*  by  Phil  Mongelluzzo  */ 
/*  */ 
/*  Revision:  08.15.88  */ 
/*  */ 
/*  Star  Guidance  Consulting.  Inc.  */ 
/*  273  Windy  Drive  */ 
/*  Waterbury,  Connecticut  06705  */ 
/*  V 
/*  (203)  574-2449  */ 
/*  */ 


/ft*************************************************************************** / 

/* 

/* 

**  WINDOWS  -  Simple  but  Elegant  Window  Functions 
**  (Lattice,  Computer  Innovations,  Microsoft, 

**  Datalight,  Aztec,  Watcoro,  Mix  Power  C) 

** 

**  Copyright  (c)  1984,  1985,  1986  -  Philip  A.  Mongelluzzo 
**  All  rights  reserved. 

** 

*/ 


(Small  Model) 
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Hr* 

**  LCS  -dLC3  sourcef ile_name 

** 

**  Datalight 

** 

**  DLC1  -dDLC  -p  -w  -mS  sourcef ile_name 

**  DLC2  sourcef ilename  -S 

** 

**  Computer  Innovations  CI86 

** 

**  ccl  -cm  -dC86  sourcef ile_name 

**  cc2 

*  * 

*/ 

/* 

*«  Computer  Innovations  comes  first.... 

*/ 


#ifdef  C86 

#def ine 

BORLAND 

0 

#def  me 

MSCV3 

0 

♦define 

MSCV4 

0 

♦define 

MSC 

0 

♦define 

MSC3 

0 

♦define 

MSC  A 

0 

♦define 

DLC 

0 

♦define 

CI86 

1 

♦define 

LC2 

0 

♦define 

LC3 

0 

♦define 

MIXPC 

0 

♦define 

AZTEC 

0 

♦define 

WATCOM 

0 

♦ifdef  . 

C86  BIG 

♦define 

LPTR 

1 

♦define 

SPTR 

0 

♦else 

♦define 

LPTR 

0 

♦define 

SPTR 

1 

♦endi f 

♦define 

UTTICE 

0 

♦define 

void  int 

struct  WGRDREGS 

{ 

unsigned  int  ax; 
unsigned  int  bx; 
unsigned  int  cx; 
unsigned  int  dx; 
unsigned  int  si; 
unsigned  int  di ; 
unsigned  int  ds; 
unsigned  int  et; 
unsigned  int  fiags; 

); 

struct  BYTERBGS  { 

unsigned  char  al,  ah; 
unsigned  char  bl,  bh; 
unsignad  char  cl,  ch; 
unsigned  char  dl,  dh; 
); 

union  REGS  { 

struct  WCRDREGS  x; 
struct  BYTEREGS  h; 

); 

struct  SRBGS  { 

unsigned  int  cs; 
unsigned  int  as; 
unsigned  int  da; 
unsigned  int  ea; 

); 

extern  unsigned  «ns  atypeO; 

#endif 


/*  define  void  as  int  */ 

/*  register  layout  is  */ 

/*  different  from  the  rest  ! 


/*  <»  KB  */ 
/*  <-  NB  */ 


/*  make  everyone  happy  */ 
/*  end  C86  Stuff  V 


!  */ 
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♦  lf  WATCCM 

♦pragma  aux  v_stksp  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  (ax  bx  cx  dx  es] ; 

♦pragma  aux  _putca  parm  caller  []  \ 

value  struct  float  struct  routine  (ax]  \ 
modify  (ax  bx  cx  dx  es] ; 

♦pragma  aux  _getca  parm  caller  (]  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

♦pragma  aux  _absloc  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

♦pragma  aux  v_wca  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es); 

♦pragma  aux  v_wtty  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es] ; 

♦pragma  aux  v_cls  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

♦pragma  aux  v_spage  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

♦pragma  aux  v_smode  parm  caller  [ )  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es); 

♦pragma  aux  vlocate  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es] ; 

♦pragma  aux  v_hidec  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es}; 

♦pragma  aux  v_aapu  parm  caller  []  \ 

value  atruct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

♦pragma  aux  v_sapd  parm  caller  [ )  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es); 

♦pragma  aux  v_rcpos  parm  caller  [ )  \ 

value  atruct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es); 

♦pragma  aux  v_rcvs  parm  caller  (]  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

♦prags  aux  v_getch  parm  caller  [ ]  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es] ; 

♦pragma  aux  v_kflush  parm  caller  (]  \ 

value  struct  float  struct  routine  (ax)  \ 
modify  [ax  bx  cx  dx  as]; 

♦pragma  aux  v_kstat  parm  caller  ()  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  (ax  bx  cx  dx  as] ; 
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♦pragma  aux  v_sctype  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  esj; 

♦pragma  aux  xferdata  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

#pragma  aux  vborder  parm  caller  [)  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es] ; 

♦pragma  aux  _vidblt  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es]; 

♦pragma  aux  _putca  parm  caller  []  \ 

value  struct  float  struct  routine  [ax]  \ 
modify  [ax  bx  cx  dx  es); 

#endif 


/* 

**  Microsoft  4.0,  5.X,  QuickC,  PowerC 

*/ 

♦if  MSCV4 
♦define  MSC  1 

♦define  MSCV3  0 

♦define  MSC 4  1 

♦define  DLC  0 

♦define  Cl  86  0 

♦define  LC2  0 

♦define  LC3  0 

♦define  BORLAND  0 

♦ifndef  MIXPC 
♦define  MIXPC  0 
♦end if 

♦define  AZTEC  0 
♦ifdef  MI86SM 
♦define  SPTR  1 

♦define  LPTR  0 

♦end if 

♦ifdef  M_I86LM 
♦define  SPTR  0 

♦define  LPTR  1 

♦end if 

♦ifdef  M_I86CM 
♦define  SPTR  0 

♦define  LPTR  1 

♦end if 

♦ifdef  M_I8»W 
♦define  SPTR  1 

♦define  LPTR  0 

♦end if 

♦define  LATTICE  1 
♦end if 

/* 

**  Microsoft  3.00 

•/ 


♦if  MSCV3 
♦define  MSC  1 

♦define  MSC 4  o 

♦define  DLC  0 

♦define  CI86  0 

♦define  LC2  0 

♦define  LC3  0 

♦define  BORLAND  0 

♦define  MIXPC  0 

♦ifndef  AZTEC 
♦define  AZTEC  0 


/*  small  code,  small  data  */ 


/*  large  code,  large  data  */ 


/*  small  code,  large  data  */ 


/*  large  code,  small  data  */ 
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#endif 

#1 fdef  M_I86SM 
♦define  SPTR  1 

♦define  LPTR  0 

♦end if 

♦ifdef  M_I86LM 
♦define  SPTR  0 

♦define  LPTR  1 

♦endif 

♦ifdef  M_I86CM 
♦define  SPTR  0 

♦define  LPTR  1 

♦endif 

♦ifdef  M_I86t*l 
♦define  SPTR  1 

♦define  LPTR  0 

♦endif 

♦define  LATTICE  1 
♦endif 

/* 

**  BORLAND 
*/ 

♦ifdef  _ TURBOC _ 

♦ifndef  BORLAND 
♦define  BORLAND  1 
♦endif 
♦endif 

♦if  BORLAND 
♦define  MSC  1 

♦define  MSC*.  1 

♦define  DLC  0 

♦define  CI86  0 

♦define  LC2  0 

♦define  LC3  0 

♦define  MIXPC  0 

♦define  AZTEC  0 

♦ifdef  _ SMALL _ 

♦define  SPTR  1 

♦define  LPTR  0 

♦endif 

♦ifdef  _ LARGE _ 

♦define  SPTR  0 

♦define  LPTR  1 

♦endif 

♦ifdef  _ COMPACT _ 

♦define  SPTR  0 

♦define  LPTR  1 

♦endif 

♦ifdef  _MEDIUM_ 
♦define  SPTR  1 

♦define  LPTR  0 

♦endif 

♦define  LATTICE  1 
♦endif 

♦define  TRUE  1 

♦define  FALSE  0 

♦define  The  BOSS  TRUE 


♦if  MSCV3  i  MSCV*.  j  BORLAND 
♦define  LINT  ARCS 
♦if  BORLAND  ~ 

♦include  "elloc.h" 

♦else 

♦if  BORLAND  |  MSC  j  DLC  j  LC2  |  LC3 
♦if  AZTEC 

♦cite 

♦include  "Belloc. h" 


/*  small  code,  small  data  */ 

/*  large  code,  large  data  */ 

/*  small  code,  large  data  */ 

/*  large  code,  small  data  */ 


/*  small  code,  small  data  */ 

/*  large  code,  large  data  V 

/*  small  code,  large  data  */ 

/*  large  code,  small  data  */ 

/*  truth  */ 

/*  lies  */ 

/*  convienent  equate  * / 

/*  Microsoft  C  or  BORLAND  */ 
/*  enforce  type  checking  */ 

/*  Borland  */ 

|  MIXPC  J  HATCOM  |  Cl 86 

/*  for  aalloc. . .  MSC  */ 
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♦endif 

#endif 

♦endif 

♦else 

char  *malloc(),  ‘callocO; 
#endi£ 


/*  if  not  Microsoft  C  */ 

/*  keep  everybody  happy  */ 


#include  "stdio.h" 

♦if  AZTEC 

♦include  "stdlib.h" 

#endif 

#if  BORLAND  !  MSC  |  DLC  |  LC2  J  LC3 

♦if  AZTEC 

♦else 

♦include  "dos.h" 

♦end if 
♦endi f 

♦include  "ctype.h" 

♦if  MSC 4 

♦include  "stdarg.h" 

♦endif 

♦if  AZTEC 


/*  standard  header  »/ 

|  MIXPC  !  WATCOM  !  Cl  66 

/*  Lattice  stuff  */ 

/*  character  conversion  stuff  */ 
/*  variable  arg  list  marcos  */ 

/*  AZTEC  DOS.H  */ 


struct  WQRDREGS  { 
unsigned  int  ax; 
unsigned  int  bx; 
unsigned  int  cx; 
unsigned  int  dx; 
unsigned  int  si ; 
unsigned  int  di ; 
unsigned  int  cflag; 

); 

struct  BYTEREGS  { 

unsigned  char  al,  ah; 
unsigned  char  bl,  bh; 
unsigned  char  cl,  ch; 
unsigned  char  dl,  dh; 

); 

union  REGS  { 

struct  WQRDREGS  x; 
struct  BYTEREGS  h; 

}; 

struct  SREGS  { 

unsigned  int  cs; 
unsigned  int  ss; 
unsigned  int  da; 
unsigned  int  es; 

); 

♦define  FP_SEG(fp)  (‘((unsigned  *)&(fp)  +  D) 

♦define  FP_OFF(fp)  (‘((unsigned  *)&(fp))) 

struct  RS  { 

int  ax,  bx,  cx,  dx,  ai,  di ,  ds,  es; 

); 


♦endif 

/* 

End  AZTEC  DOS.H  */ 

♦define 

SAVE 

TRUE 

/* 

similar  truth  */ 

♦define 

RESTORE 

FALSE 

/* 

fibs  •/ 

♦define 

PAINT 

TRUE 

/* 

screen  update  modes  */ 

♦define 

FLASH 

FALSE 

/* 

ditto  */ 

♦define 

REPLACE 

1 

/* 

for  flicker  free  */ 

♦define 

ERASE 

0 

/* 

scroll  w_sapd  A  wsapu  */ 

♦define 

FAST 

0x01 

/* 

fast  retrace  ‘/ 

♦define 

F'O W 

0x06 

/* 

slow  retrace  •/ 

♦define 

NULPTR 

(char  *)  0 

/* 

null  pointer  */ 

♦define 

NUL 

•\0' 

/• 

NUL  char  */ 
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♦define 

NVAL 

0 

/* 

null  value  -  INTEGER  */ 

♦define 

BEL 

0x07 

/* 

beep  */ 

♦define 

BS 

0x08 

/* 

backspace  */ 

♦define 

ESC 

Oxlb 

/• 

Escape  */ 

♦define 

CR 

OxOd 

/* 

carriage  return  */ 

♦define 

LF 

0x0a 

/* 

linefeed  */ 

♦define 

RUB 

0x7f 

/* 

delete  */ 

♦define 

NAK 

0x15 

/* 

'U  */ 

♦define 

ETX 

0x03 

/* 

*C  */ 

♦define 

CAN 

0x18 

/* 

*X  */ 

♦define 

Del 

0x53 

/* 

Del  key  scan  code  */ 

♦define 

ECHO 

0x8000 

/* 

echo  disable  bit  */ 

♦define 

BIOS 

0x01 

/* 

BIOS  Scrolling  */ 

♦define 

DMAS 

0x02 

/* 

The  BOSS’S  DMA  Scrolling  */ 

/* 

**  Externals 
*/ 


extern  int  wn_dma£lg; 
extern  char  wn_sbit; 
extern  int  wnblank ; 
extern  int  wns_bchars[ ] ; 
extern  unsigned  int  vmsmtflg; 
extern  int  wns_cflag; 

extern  struct  SREGS  wns_srp; 

extern  unsigned  wni  segM; 
extern  unsigned  wni_off[]; 
extern  unsigned  vmi_ptr[); 


#de£ine  BCUL 
♦define  BCUR 
♦define  BCTB 
♦define  BCSD 
♦define  BCLL 
♦define  BCLR 

/* 

**  Misc  Stuff 

*/ 


wns_bchars[0] 
wns_bchars[l] 
wns_bchars [2] 
wns_bchars[3] 
wns_bchars[4) 
wns_bchars[5] 


♦if  LC2 

extern  unsigned  wnsmtype ( ) ; 
♦endif 


/*  dma  flag  */ 

/*  retrace  test  bit  8  slow,  1  fast  */ 
/*  vidon  &  vidoff  control  flag  */ 

/*  box  chars  */ 

/*  monitor  type  flag  */ 

/*  close  in  progress  flag  */ 

/*  for  segread  */ 

/*  for  wns_push/pop  */ 

/*  ditto  */ 

/*  ditto  */ 

/*  some  shorthand  for  later  V 


/*  make  everyone  happy  */ 


♦define  WMR  wn->bsize 


/*  shorthand  */ 


typedef  struct  web 
{ 

int  ulx, 
uly, 
xsize, 
ysize, 

ccx, 

ccy , 
style, 
bstyle, 
bsize; 

char  ‘semsave; 
int  page, 

oldx, 

oldy , 
wrpf lg , 
aynflg; 

char  *handle; 

char  *pravptr; 
char  *nextptr ; 
unsigned  tmpsag; 
unsigned  tmpoff; 


/*  Window  control  block  */ 

/*  upper  left  corner  x  coordinate  */ 

/*  upper  left  corner  y  coordinate  */ 

/*  width  of  window  -  INSIDE  dimension  */ 
/*  height  of  window  -INSIDE  dimension  */ 
/*  virtual  cursor  offset  in  window  */ 

/*  attribute  to  be  used  in  window  */ 

/*  border  attribute  */ 

/*  total  border  size  0  or  2  only  */ 

/*  pointer  to  screen  save  buffer  */ 

/*  current  video  page  being  used  •/ 

/*  cursor  position  when  window  was  */ 

/*  opened  (used  for  screen  restore)  */ 

/*  wrap  flag  */ 

/*  cursor  sync  flag  */ 

/*  my  own  id  */ 

/*  linked  list  -  previous  */ 

/*  linked  list  -  next  */ 

/*  for  activate  */ 

/*  ditto  */ 
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int  smeth;  /*  scroll  method  to  use  */ 

}  WINDOW,  *WINDOWPTR; 

extern  WINDCMPTR  wns_last;  /*  last  window  opened  */ 

♦if  MSCV3  j  MSCV4  j  BORLAND  j  DLC  j  LC3  /*  allow  for  LINT_ARGS  */ 
♦ifndef  GENFNS 

♦include  "windows . fns"  /*  enforce  type  checking  */ 

♦end if 

♦else  /*  and  almost  lint  args  */ 

struct  web  *wn_open(); 

struct  web  *wn_move(); 

struct  web  *wn_save(); 

char  *wn_gets ( ) ; 

struct  web  *wn_save(); 

struct  wi_scb  *  *wn_frmopnO ; 

unsigned  int  wns_mtype(); 

♦endif 


♦define 

BLACK 

0x00 

/* 

foreground  */ 

♦define 

RED 

0x04 

/* 

background  */ 

♦define 

GREEN 

0x02 

/* 

colors  */ 

♦define 

YELLOW 

0x14 

/* 

bg  «  4  |  fg  */ 

♦define 

BLUE 

0x01 

♦define 

MAGENTA  0x05 

♦define 

CYAN 

0x03 

♦define 

WHITE 

0x07 

♦define 

BROWN  0x2e 

♦define 

BLINK 

0x80 

♦define 

BOLD 

0-08 

♦define 

NDISPB 

0x00 

/* 

non  display  black  */ 

♦define 

NDISFW 

0x77 

/* 

non  display  white  */ 

♦define 

RVIDEO 

0x70 

/* 

reverse  video  */ 

♦define 

UNLINE 

0x01 

/* 

under  line  (BLUE)  */ 

♦define 

NVIDEO 

0x07 

/* 

normal  video  */ 

♦define 

NORMAL 

0x03 

/• 

cyan  is  normal  for  me  */ 

/* 

**  Display  Mode 
V 

;  Atributes 

♦define 

B4025 

0 

/* 

black  &  white  40  x  25  */ 

♦define 

C4025 

1 

/* 

color  40  x  25  */ 

♦define 

B8025 

2 

/* 

black  &  white  80  x  25  */ 

♦define 

C8025 

3 

/* 

color  80  x  25  */ 

♦define 

C320 

4 

/* 

color  graphics  320  x  200 

♦define 

B320 

5 

/* 

black  &  white  graphics  * 

♦define 

HIRES 

6 

/* 

B&W  hi  res  640  *  200  V 

♦define  MONO 

7 

/* 

monocrocne  80  x  25  */ 

/* 

**  Macro  to  sat  attribute  byte 

*/ 


♦define  v_setatr(bg,fg,blink,bold)  ( (blink  |  (bg«4 ))  |  (fg 'bold) ) 

/* 

**  Key  Scan  Scodes  &  Window  Input  Stuff 
V 


♦define 

BELL 

0x0007 

/* 

ring  e  ding  */ 

♦define 

ENTER 

OxOd 

/* 

enter/return  key  */ 

♦define 

LARROW 

0x4b00 

/* 

left  arrow  */ 

♦define 

RARROW 

0x4d00 

/* 

right  arrow  •/ 

♦define 

DARROW 

0x5000 

/* 

down  arrow  */ 

♦define 

U  ARROW 

0x4800 

/• 

up  arrow  */ 

♦define 

HOME 

0x4700 

/• 

home  key  */ 

♦define 

END 

0x4 fOO 

/* 

end  key  */ 

♦define 

PAGEUP 

0x4900 

f* 

page  up  key  */ 

♦define 

PAGEDN 

0x5100 

/* 

page  down  key  */ 

♦define 

INS 

0x5200 

/• 

inaert  key  */ 
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♦define  DEL  0x5300 

♦define  FI  Ox3bOO 

♦define  HELP  FI 

♦define  TAB  0x0f09 

♦define  BKTAB  OxOfOO 


♦define  GDONE  0 
♦define  GDATE  10 
♦define  GTIME  11 
♦define  GINT  12 
♦define  GUINT  13 
♦define  GLONG  14 
♦define  GFLOAT  15 
♦define  GPHONE  16 
♦define  GTEXT  17 
♦define  GBOOL  18 
♦define  DTEXT  19 


♦define  NSTR  "" 

♦define  NFRM  (WIFQRM)(0) 
♦define  NFLD  0 
♦define  SET  1 
♦define  XEQ  2 
♦define  MAXSTR  80 

union  wiargs  { 
int  vi ; 
int  *vip; 
unsigned  int  vui ; 
unsigned  int  *vuip; 
char  vc ; 
char  *vcp ; 
long  vl: 
long  *vlp; 
float  vf; 
float  *vfp: 
double  vd; 
double  *vdp ; 

}  ; 

typedef  struct  wi_scb  { 
char  *pself; 
int  fcode; 

WINDOWPTR  wn; 
int  row; 
int  col; 
char  *prmpt; 
unsigned  int  atrib; 
char  fill; 
union  wi_aigs  vl; 
union  wi_args  v2; 
union  wi_args  v3; 
union  wi_args  v4; 
union  wi_args  v5; 
union  wi_args  v6; 
union  wi_args  v7 ; 
union  wi_args  v8; 

}  WIFLD,  *WIFLDPTR,  **WIF0RM; 

♦define  WNLPTR  (WINDOWPTR)  0 


/*  delete  key  */ 

/*  FI  aka  HELP  */ 

/*  same  as  FI  */ 

/*  tab  */ 

/*  back  (shift)  tab  */ 

/*  0  to  100  are  reserved!!  */ 

/*  end  of  list  */ 

/*  wn_gdate  */ 

/*  wn_gtime  */ 

/*  wn_ gint  */ 

/*  wn_guint  */ 

/*  wn_glong  */ 

/*  wn_gfloat  */ 

/*  wn_gphone  */ 

/*  wn_gtext  */ 

/*  wn_gbool  */ 

/*  Display  text  only  */ 

/*  from  above-100  are  reserved! !  */ 

/*  null  string  */ 

/*  null  form  pointer  */ 

/*  must  be  int  0  */ 

/*  form  setup  */ 

/*  immediate  execution  */ 

/*  max  size  -  vm_gtext,  wn_input  */ 

/*  variable  arg  type  union  */ 

/*  (int)  */ 

/*  (int  *)  •/ 

/*  (unsigned  int)  */ 

/*  (unsigned  int  *)  */ 

/*  (char)  */ 

/*  (char  *)  */ 

/*  (long)  */ 

/*  (long  *)  V 
/»  (float)  */ 

/*  (float  *)  */ 

/*  (double)  »/ 

/*  (double  *)  */ 


/*  screen  control  block  */ 

/*  pointer  to  myself  */ 

/*  input  funtion  code  */ 

/*  the  window  */ 

/*  window  (wn)  location  -  row  */ 
/*  window  (wn)  location  -  col  */ 
/*  prompt  string  for  field  */ 

/*  input  field  attribute  */ 

/*  input  field  fill  character  */ 
/*  whatever  */ 

/*  whatever  */ 

/*  whatever  */ 

/*  whatever  */ 

/*  whatever  */ 

/*  whatever  */ 

/*  whatever  •/ 

/*  whatever  */ 


/*  A  TRUE  NULL  WINDOW  POINTER  */ 


/*  End  */ 


